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the  vertical  structure  and  temporal  behavior  of  Internal  waves,  par- 
ticularly those  with  periods  near  the  local  inertial  period.  An  im- 
portant feature  of  the  observed  structure  is  the  polarizatldrof  hor- 
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with  temporal  changes  of  the  vertical  wave  structure  seen  In  a time 
series  of  profiles  made  at  one  location,  has  been  related  to  the  di- 
rection of  vertical  energy  flux  due  to  the  observed  waves.  Whereas 
the  observed  vertical  phase  propagation  can  be  affected  by  horizontal 
advection  of  waves  past  the  point  of  observation,  the  use  of  wave  po- 
larization to  Infer  the  direction  of  vertical  energy  propagation  has 
the  advantage  that  it  is  not  influenced  by  horizontal  advection.  The 
result  shows  that  at  a location  where  profiles  were  obtained  over 
smooth  topography,  the  net  energy  flux  was  downward,  indicating  that 
the  energy  sources  for  these  waves  were  located  at  or  near  the  sea 
surface.  An  estimate  of  the  net,  downv/ard  energy  flux  (~  .2  - .3 
erg/cmVsec)  has  been  obtained.  Calculations  have  been  made  which 
show  that  a frictional  bottom  boundary  layer  can  be  an  important  en- 
ergy sink  for  near-inertial  waves.  A rough  estimate  suggests  that 
the  observed,  net,  downward  energy  flux  could  be  accounted  for  by 
energy  losses  In  this  frictional  boundary  layer.  A reflection 
coefficient  for  the  observed  waves  as  they  reflect  off  the  bottom 
has  been  estimated. 

In  contrast,  some  profiles  made  over  a region  of  rough  topog- 
raphy indicate  that  the  rough  bottom  may  also  be  acting  to  generate 
near-inertial  waves  which  propagate  energy  upv/ard. 

Calculations  of  vertical  flux  of  horizontal  kinetic  energy, 
using  an  empirical  form  for  the  energy  spectrum  of  internal  waves, 
show  that  this  vertical  flux  reaches  a maximum  for  frequencies  10X  - 

20%  greater  than  the  local  inertial  frequency.  Comparison  with  pro-  ■ 

filer  velocity  data  and  frequency  spectra  supports  the  conclusion 

that  the  dominant  waves  had  frequencies  10%  - 20%  greater  than  the 

inertial  frequency.  The  fact  that  the  waves  were  propagating  energy 

in  the  vertical  is  proposed  as  the  reason  for  the  observed  frequency 

shift. 

Finally,  energy  spectra  in  vertical  wave  number  have  been 
calculated  from  the  profiles  in  order  to  compare  the  data  with  an 
empirical  model  of  the  energy  density  spectrum  for  internal  waves 
proposed  by  C.  Garrett  and  W.  Hunk  (1975).  The  result  shows  that 
although  the  general  shape  and  magnitude  of  the  observed  spectrum 
compares  well  with  the  empirical  model,  the  two-sided  spectrum  is 
not  symmetric  in  vertical  wave  number.  This  asymmetry  has  been 
used  to  infer  that  more  energy  was  propagating  downward  than  upward. 

These  calculations  have  also  been  used  to  obtain  the  coherence  be- 
tween profiles  made  at  the  same  location,  but  separated  in  time 
(the  so-called  dropped,  lagged,  rotary  coherence).  This  coherence 
is  compared  with  the  aforementioned  empirical  model.  The  coherence 
results  show  that  the  contribution  of  the  semidiurnal  tide  to  the 
energy  of  the  profiles  is  restricted  to  long  vertical  wave  lengths. 
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ABSTRACT 

A set  of  vertical  profiles  of  horizontal  ocean  currents,  ob- 
tained by  electro-magnetic  profilers  in  the  Atlantic  Ocean  southwest 
of  Bermuda  in  the  spring  of  1973,  has  been  analyzed  In  order  to  study 
the  vertical  structure  and  temporal  behavior  of  Internal  waves,  par- 
ticularly those  with  periods  near  the  local  Inertial  period.  An  im- 
portant feature  of  the  observed  structure  is  the  polarization  of  hor- 
izontal velocity  components  in  the  vertical.  This  polarization,  along 
with  temporal  changes  of  the  vertical  wave  structure  seen  in  a time 
series  of  profiles  made  at  one  location,  has  been  related  to  the  di- 
rection of  vertical  energy  flux  due  to  the  observed  waves.  Whereas 
the  observed  vertical  phase  propagation  can  be  affected  by  horizontal 
advectlon  of  waves  past  the  point  of  observation,  the  use  of  wave  po- 
larization to  Infer  the  direction  of  vertical  energy  propagation  has 
the  advantage  that  it  Is  not  influenced  by  horizontal  advection.  The 
result  shows  that  at  a location  where  profiles  were  obtained  over 
smooth  topography,  the  net  energy  flux  was  downward,  indicating  that 
the  energy  sources  for  these  waves  were  located  at  or  near  the  sea 
surface.  An  estimate  of  the  net,  downward  energy  flux  (~  .2  - .3 
erg/cm2/sec)  has  been  obtained.  Calculations  have  been  made  which 
show  that  a frictional  bottom  boundary  layer  can  be  an  important  en- 
ergy sink  for  near-inertial  waves.  A rough  estimate  suggests  that 
the  observed,  net,  downward  energy  flux  could  be  accounted  for  by 
energy  losses  In  this  frictional  boundary  layer.  A reflection 
coefficient  for  the  observed  waves  as  they  reflect  off  the  bottom 
has  been  estimated. 

In  contrast,  some  profiles  made  over  a region  of  rough  topog- 
raphy Indicate  that  the  rough  bottom  may  also  be  acting  to  generate 
near-inertial  waves  which  propagate  energy  upward. 

Calculations  of  vertical  flux  of  horizontal  kinetic  energy, 
using  an  empirical  form  for  the  energy  spectrum  of  internal  waves, 
show  that  this  vertical  flux  reaches  a maximum  for  frequencies  10%  - 
20%  greater  than  the  local  Inertial  frequency.  Comparison  with  pro- 
filer velocity  data  and  frequency  spectra  supports  the  conclusion 
that  the  dominant  waves  had  frequencies  10%  - 20%  greater  than  the 
Inertial  frequency.  The  fact  that  the  waves  were  propagating  energy 
In  the  vertical  Is  proposed  as  the  reason  for  the  observed  frequency 
shift. 


Finally,  energy  spectra  In  vertical  wave  number  have  been 
calculated  from  the  profiles  in  order  to  compare  the  data  with  an 
empirical  model  of  the  energy  density  spectrum  for  internal  waves 
proposed  by  C.  Garrett  and  W.  Munk  (1975).  The  result  shows  that 
although  the  general  shape  and  magnitude  of  the  observed  spectrum 
compares  well  with  the  empirical  model,  the  two-sided  spectrum  is 
not  symnetric  In  vertical  wave  number.  This  asymmetry  has  been 
used  to  infer  that  more  energy  was  propagating  downward  than  upward. 
These  calculations  have  also  been  used  to  obtain  the  coherence  be- 
tween profiles  made  at  the  same  location,  but  separated  in  time 
(the  so-called  dropped,  lagged,  rotary  coherence).  This  coherence 
is  compared  with  the  aforementioned  empirical  model.  The  coherence 
results  show  that  the  contribution  of  the  semidiurnal  tide  to  the 
energy  of  the  profiles  is  restricted  to  long  vertical  wave  lengths. 
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Chapter  I Introduction  and  Historical  Review 

For  many  years  the  subject  of  internal  waves  has  been  of  in- 
tense interest  to  meteorologists,  aeronomists,  and  oceanographers. 
These  internal  waves  are  oscillations  in  the  atmosphere  or  ocean 
that  are  supported  by  buoyancy  and  pressure  forces,  and  also  by  the 
fictitious  Coriolis  force  introduced  by  the  rotation  of  the  earth. 
Oscillations  having  time  and  length  scales  appropriate  to  internal 
waves  have  been  observed  at  many  locations  in  the  ocean.  In  the 
atmosphere,  motions  having  time  and  space  scales  appropriate  to 
internal  waves  have  also  been  observed.  However,  until  recently, 
meteorologists  and  oceanographers  have  been  studying  internal  wave 
phenomena  from  decidedly  different  points  of  view.  This  has  been 
brought  about  primarily  by  the  different  types  of  instrumentation 
used  in  the  two  disciplines.  Observations  of  atmospheric  internal 
waves  are  discussed  first. 

In  the  years  since  World  War  II,  several  methods  have  been 
developed  by  meteorologists  for  observing  variables  such  as  wind 
and  temperature  in  the  atmosphere.  Wind  measurements,  which  are 
of  primary  interest  in  the  present  study,  have  been  performed  using 
the  following  methods:  tracking  the  distortion  of  smoke  (actually 

sodium)  trails  left  by  high  altitude  rockets;  tracking  of  radar  re- 
flecting balloons  with  ground-based  radar;  visual  tracking  of  meteor 
trails  by  telescopes;  tracking  chaff  (strips  of  aluminized  reflect- 
ing material)  with  ground-based  radar;  and  radio-echo  tracking  of 


W 


meteor  trails.  All  of  these  methods  have  their  disadvantages. 
Rockets,  for  example,  are  expensive  to  use,  and  they  have  not 
usually  been  applied  to  rapidly  repeated  measurements  in  the  at- 
mosphere. Visual  and  radio-echo  tracking  of  meteors  depend,  of 
course,  on  a sufficient  number  of  meteors  entering  the  upper  at- 
mosphere. Except  in  rare  cases  when  meteor  showers  occur,  the 
number  of  meteors  in  the  atmosphere  has  usually  been  insufficient 
to  study  high-frequency  oscillations  in  the  wind  field,  such  as 
internal  waves.  More  recently,  tracking  of  ascending,  radar  re- 
flecting balloons  has  provided  better  measurements  of  atmospheric 
winds . 

All  of  these  methods,  although  diverse  in  concept,  have  in 
common  the  fact  that  they  provide  profiles  of  horizontal  atmos- 
pheric winds  as  a function  of  altitude.  Many  examples  of  these 
profiles  have  been  presented  in  the  literature  (see,  for  example, 
the  paper  by  L.  Broglio  in  "Proceedings  of  the  First  International 
Symposium  on  Rocket  and  Satellite  Meteorology,”  in  which  results 
from  rocket  profiles  are  presented,  and  the  radar  balloon  observa- 
tions of  Endlich,  Singleton  and  Kaufman  (1969)).  Thus,  the  atten- 
tion of  meteorologists  and  aeronomists  has  historically  been  di- 
rected toward  observing  the  vertical  structure  of  the  atmospheric 
wind  field.  Attempts  at  using  repeated  profiles  to  determine  the 
amount  of  energy  contributed  by  hiqh-frequency  (internal  wave  and 
turbulent)  oscillations  at  different  time  scales  or  frequencies 
have  been  notably  less  successful  than  the  measurements  of  the 
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profiles  themselves.  These  attempts  have  usually  consisted  of  try- 
ing to  fit  diurnal  and  semidiurnal  period  sinusoids  to  the  winds 
observed  at  a given  altitude  (El ford  and  Robertson,  1953;  Smith, 
1960). 

The  situation  In  the  oceans  has  been  substantially  different. 
Until  recently,  the  primary  tools  which  oceanographers  have  used  to 
observe  oscillations  with  frequencies  in  the  internal  wave  range 
have  been  fixed-location  current  meters,  temperature  sensors,  and 
neutrally  buoyant  floats.  The  current  meter  allows  one  to  measure, 
at  a fixed  point  in  the  ocean,  the  horizontal  current  vector  as  a 
function  of  time.  If  current  meters  are  placed  at  many  locations 
in  a part  of  the  ocean  (by  being  put  on  moorings),  then  time  meas- 
urements In  a spatial  array  are  obtained.  These  arrays  allow  one 
to  obtain  good  temporal  resolution  of  velocity,  but  only  at  a 
limited  number  of  points  in  space. 

Current  meters  have  provided  a large  volume  of  data  on  the 
temporal  behavior  of  the  ocean.  Fofonoff  (1969)  showed  that  the 
behavior  of  horizontal  velocity  components  and  temperature  in  the 
internal  wave  frequency  band  (those  frequencies  corresponding  to 
periods  between  one  half  of  the  local  pendulum  day,  l/(2sin(Lat. )) 
days,  and  the  local  stable  oscillation  period  of  a particle  of 
water,  or  the  Brunt-ValsSlS  period)  could  be  reasonably  well  re- 
lated to  the  expected  behavior  of  internal  waves  in  a linear  model. 
Webster  (1968)  summarizes  many  of  the  current  meter  observations 
that  have  been  used  to  study  internal  waves  and  other  high-frequency 
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oceanic  phenomena.  Most  of  these  current  meter  records  have  shown 
that  the  dominant  contributions  to  horizontal  kinetic  enerqy  spec- 
tra come  from  frequencies  near  the  local  inertial  frequency,  cor- 
responding to  the  aforementioned  period  of  one  half  of  a pendulum 
day,  and  from  the  tides. 

Temperature  observations  have  been  made  using  several  dif- 
ferent methods.  Hauritz,  Storrmel  and  Munk  (1959)  observed  temp- 
erature oscillations  as  a function  of  time  at  two  fixed  depths, 

50  and  550  meters,  near  Bermuda.  These  data  yielded  frequency 
spectra  for  the  temperature  oscillations. 

Another  method,  employed  by  Charnock  (1965)  and  LaFond  and 
LaFond  (1971),  is  to  tow  a weighted  thermistor  chain  behind  a 
vessel.  If  the  speed  of  the  vessel  is  large  compared  to  the  hori- 
zontal phase  propagation  speed  of  the  waves  beinq  observed,  spec- 
tral analysis  of  the  temperature  trace  obtained  at  a oiven  depth 
gives  a horizontal  wave  number  spectrum  of  the  observed  isotherm 
displacements.  Similar  observations  can  be  obtained  by  towing  an 
isotherm- following  "fish"  behind  a research  vessel  (Katz,  1973, 

1975). 

The  use  of  neutrally  buoyant  floats  has  also  provided  data 
on  internal  waves.  Fochapsky  (1972)  used  such  floats  to  observe 
oscillations  of  temperature  and  horizontal  velocity  with  periods 
in  the  internal  wave  range.  Voorhis  (1968)  obtained  measurements 
of  vertical  motion  with  a neutrally  buoyant  float  which  rotates  in 
the  presence  of  vertical  flow  past  the  float.  The  rate  of  rota- 
tion of  the  float  gives  an  estimate  of  this  vertical  flow. 
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Many  of  the  observations  described  above  (as  well  as  others 
which  have  not  been  mentioned)  have  been  gathered  together  by 
Garrett  and  Munk  (1972,  1975)  In  an  attempt  to  model  the  energy 
spectrum  of  internal  waves  in  frequency  and  horizontal  wave  number. 
The  model  depends  not  only  on  Individual  spectra  of  temperature  or 
velocity,  but  also  on  the  correlation,  or  coherence,  between  obser- 
vations made  at  different  points  of  a spatial  array.  Coherence 
observations  have  also  been  gathered  together  and  described  by 
those  authors  (Garrett  and  Munk,  1972,  1975). 

A major  difficulty  with  this  type  of  model  construction  has 
been  lack  of  knowledge  concerning  the  vertical  structure  of  the 
wave  field.  This  vertical  structure  Is  difficult  to  determine 
using  velocity  (or  temperature)  measurements  from  moorings,  since 
only  a limited  number  of  instruments  can  be  placed  on  any  one 
mooring.  Thus,  vertical  profiles  of  horizontal  velocity,  which 
have  been  the  most  common  method  of  observation  in  the  atmosphere, 
have  been  the  least  common  type  of  observation  in  oceanographic 
work. 

Recently,  however,  several  methods  have  been  developed  that 
allow  the  oceanographer  to  obtain  vertical  profiles  of  horizontal 
current  in  the  ocean  that  are  analogous  to  profiles  obtained  in  the 
atmosphere.  One  such  method  involves  the  use  of  acoustic  trans- 
ponders located  on  the  ocean  bottom  to  track  a float  as  it  descends 
through  the  water  column  (Rossby,  1969;  Pochapsky  and  Malone,  1972). 

Another  method  uses  a moored  vertical  cable  with  inclinom- 
eters located  at  various  points  along  the  cable  (McNary,  1968). 
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Horizontal  currents  cause  the  cable  to  distort.  This  distortion, 
as  measured  by  the  inclinometers,  can  then  be  converted  to  values 
of  the  horizontal  velocity  component  along  the  cable.  A third 
method  uses  horizontal  electrical  currents  induced  by  the  flow  of 
ocean  water  through  the  terrestrial  magnetic  field  (Sanford,  1971). 
The  device  which  measures  these  electrical  currents  is  called  an 
electro-magnetic  velocity  profiler,  or  EMVP.  A picture  of  this 
device,  developed  by  T.  Sanford  and  R.  Drever  of  the  Woods  Hole 
Oceanographic  Institution,  is  shown  in  Figure  1.  The  EMVP  is  re- 
leased from  a research  vessel,  and  as  it  falls  to  the  ocean  bottom 
it  continuously  records  the  horizontal  component  of  electrical 
current,  as  well  as  temperature,  conductivity,  and  pressure.  When 
the  EMVP  reaches  the  bottom,  it  reverses  and  ascends  through  the 
water  column,  again  recording  the  above  variables.  There  are  no 
hard  connections  (cables,  etc.)  between  the  research  vessel  and 
the  EMVP:  the  only  communication  from  EMVP  to  ship  is  through 

acoustic  telemetry.  All  data  are  recorded  on  7-track  magnetic 
tape  within  the  EMVP.  Upon  return  to  the  ship,  the  data  on  the 
tape  is  reformatted  onto  9-track  tape,  and  several  computer  programs 
are  used  to  convert  the  internally  recorded  variables  to  their 
physical  counterparts,  such  as  temperature.  Another  routine,  using 
a conversion  equation,  converts  horizontal  electrical  current  den- 
sity to  horizontal  water  flow  relative  to  an  unknown,  but  depth 
independent,  horizontal  velocity.  A further  description  of  the 
EMVP  and  its  operation  can  be  found  in  Sanford,  et.  al . (1974). 


Figure  1.  Electro-magnetic  velocity  profiler 


The  horizontal  velocity  component  Is  not  the  only  quantity 
that  has  been  measured  by  vertical  profiling.  Vertical  profiles 
of  temperature,  conductivity  and  salinity  have  been  obtained  by 
raising  and  lowering  a CTD  (conductivity-temperature-depth  sensor). 
Such  data  have  been  presented  by  Hayes  (1975)  and  Hayes,  et.  al. 
(1975). 

We  turn  now  to  a description  of  the  experiment  which  is 
discussed  in  this  report.  In  the  spring  of  1973  a large-scale 
oceanographic  experiment  called  MODE  (Mid-Ocean  Dynamics  Experi- 
ment) was  carried  out  in  the  Atlantic  Ocean,  southwest  of  Bermuda. 

A map  showing  the  bathymetry  in  the  region  of  this  experiment  is 
shown  in  Figure  2.  The  EMVP  was  used  to  take  a large  number  of 
vertical  profiles  during  the  experiment.  The  locations,  dates  and 
times  of  all  profiles  are  given  in  Appendix  A.  We  are  going  to 
concentrate  particularly  on  two  subsets  of  all  the  profiles.  The 
first  subset  consists  of  a time  series  of  repeated  profiles  made 
at  the  center  of  the  experiment  (69°  40'  W,  27°  59'  N).  This  will 
be  called  the  "central  mooring  time  series,"  since  the  central 
mooring  of  the  MODE  experiment  was  located  near  these  coordinates. 
This  time  series,  made  up  of  20  profiles,  lasted  from  June  11-15, 
1973.  (These  profiles  are  designated  by  an  asterisk  in  Appendix  A) 
As  can  be  seen  from  Figure  2,  the  location  of  the  series  is  over  a 
region  of  smooth  bottom  topography.  The  second  subset  is  made  up 
of  some  profiles  taken  around  a ridge  east  of  the  time  series,  in  a 
region  of  rough  topography.  These  locations  are  shown  in  Figure  2. 


Figure  2.  Bathymetry  of  the  MODE  region.  Location  of  central 

mooring  time  series  is  shown  by  an  X.  The  ridge  over 
which  rough  topography  profiles  were  made  is  circled. 
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Any  profile  that  has  ever  been  made  by  the  EMVP  has  shown 
two  dominant  contributions  to  the  velocity  structure  in  the  verti- 
cal. The  first  contribution  is  a low-frequency  part,  which  does 
not  change  appreciably  from  one  profile  to  the  next  at  the  same 
location,  and  which  corresponds  well  with  current  shear  derived  by 
geostrophic  shear  calculations  in  that  region.  The  second  part  is 
a high-frequency  contribution  which  is  evidently  due  to  internal 
waves,  particularly  those  with  frequencies  near  the  local  inertial 
frequency.  It  is  this  high-frequency  part  that  will  be  of  main 
interest  in  what  follows. 

Chapter  II  gives  a description  and  analysis  of  the  velocity 
data  obtained  from  the  five-day  time  series  of  profiles.  Plots  of 
mean  flow  and  average  perturbation  horizontal  kinetic  energy  versus 
depth  are  given,  as  are  contour  plots  showing  the  evolution  of 
features  in  the  velocity  field  in  time,  and  spectra  of  the  high- 
frequency  profile  structure  in  vertical  wave  number. 

One  of  the  most  interesting  results  obtained  from  the  spec- 
tral analysis  in  Chapter  II  is  that  there  is  a definite  elliptical 
polarization  of  the  high-frequency  waves  in  the  profiles.  That  is, 
the  horizontal  current  vector  in  a profile  (after  the  mean  has  been 
removed)  clearly  tends  to  rotate  with  depth.  This  polarization  and 
its  relation  to  the  vertical  propagation  of  internal  wave  energy  is 
described  in  Chapter  III.  In  Chapter  III,  comparisons  will  be  made 
between  our  observations  and  similar  observations  made  in  the  atmos- 
phere. 


Chapter  IV  describes  the  reflection  of  Internal  waves  from 
a smooth  bottom.  The  spectra  presented  in  Chapter  II  are  used  to 
attempt  to  calculate  reflection  coefficients  for  these  waves. 

Chapter  V presents  a comparison  of  the  vertical  wave  number 
spectra  In  Chapter  II  to  the  theoretical  models  of  Garrett  and 
Munk  (1972,  1975).  The  coherence  of  velocity  components  between 
pairs  of  profiles  separated  by  a time  lag  is  also  presented.  Ver- 
tical energy  flux  of  the  observed  waves  is  discussed  in  more  detail. 

Chapter  VI  discusses  the  profiles  obtained  over  rough  topog- 
raphy. Chapter  VII  provides  some  conclusions,  a further  discussion 
of  some  of  the  data,  and  recommendations  for  future  experiments. 
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Chapter  II  Description  and  Analysis  of  Data  from  the  Five-day 
Time  Series  of  Profiles 

It  was  mentioned  In  the  Introduction  that  EMVP  data  is  his- 
torically more  akin  to  meteorological  measurements  in  the  upper  at- 
mosphere than  to  oceanographic  data.  In  particular,  profiler  data 
emphasizes  the  vertical  structure  of  the  ocean  or  atmosphere,  while 
until  recently  most  oceanographic  work  has  emphasized  the  temporal 
behavior  of  oceanic  motions.  This  chapter  will  present.  In  various 
forms,  horizontal  velocity  data  obtained  during  the  five-day  series 
of  profiles,  and  will  describe  some  of  the  methods  used  to  reduce 
and  analyze  this  data.  Succeeding  chapters  are  devoted  to  the 
Interpretation  of  results  obtained  from  this  analysis.  The  reader 
may  refer  to  Appendix  B for  a brief  discussion  of  internal  wave 
measurements  made  by  the  EMVP  and  to  Appendix  C for  a discussion  of 
the  accuracy  and  precision  of  data  presented  in  the  following  sec- 
tions. 

a.)  Dominance  of  inertial -period  motions  in  EMVP  profiles--the 
"mirror-imaging"  of  profiles 

Figure  3 shows  the  east  and  north  components  of  two  profiles, 
219D  and  221 D . (The  letter  D signifies  that  the  down  portion  of  a 
drop  was  used.  The  letter  U will  be  used  to  denote  the  up  portion 
of  a drop.)  The  time  separation  between  the  two  profiles  is  approx- 
imately one  half  of  an  inertial  period.  The  vertical  sampling  rate 
is  one  point  every  10  dbar.  It  is  evident  from  the  figure  that, 
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superimposed  on  a low-frequency  shear  in  the  profiles,  there  is  a 
considerable  amount  of  high-frequency  activity.  This  high-frequen- 
cy activity  is  evident  in  the  almost  total  reversal  of  the  velocity 
structure  over  the  time  interval.  If  we  remove  the  low-frequency 
shear  of  the  profiles  (by  averaging  all  drops  in  the  time  series 
and  subtracting  the  low-frequency  shear  profile  from  each  drop  sep- 
arately), we  find  that  the  high-frequency  part  that  is  left  over 
is  negatively  correlated  in  profiles  made  one  half  of  an  inertial 
period  apart.  This  negative  correlation,  or  "mirror-imaging,"  is 
an  indication  that  the  EMVP  profiles  are  strongly  dominated  in  the 
high  frequencies  by  motions  with  periods  close  to  the  local  iner- 
tial period.  It  will  be  shown  later  that  some  of  the  observed 
variability  is  also  contributed  by  the  semidiurnal  tide. 

b.)  Low-frequency  "geostrophic"  profiles 

For  ease  of  discussion  we  will  call  oscillations  with  fre- 
quencies less  than  the  local  inertial  frequency  "low-frequency" 
oscillations,  or  the  "geostrophic"  part  of  the  flow.  Similarly, 
oscillations  with  frequencies  greater  than  or  equal  to  the  inertial 
frequency  will  be  called  "high-frequency,"  or  "internal  wave" 
oscillations.  It  should  be  clear  that  these  terms  are  used  solely 
to  simplify  the  description  of  the  data.  We  do  not  intend  to 
imply  that  all  motions  with  frequencies  less  than  inertial  are  in 
geostrophic  balance.  Also,  high-frequency  motions  are  not  necessar- 
ily due  to  internal  waves.  A certain  amount  of  high-frequency  en- 
ergy may  be  contributed  by  turbulence,  for  example.  It  is  known, 
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however,  that  there  is  a "well,"  or  depression  in  frequency  spectra 
obtained  in  the  r’ODF  reqion  (see  Fioure  29).  This  well  exists  in 
the  ranqe  of  periods  from  the  inertial  period  to  periods  of  several 
days.  Therefore,  averaqino  velocity  data  over  several  days  will 
partition  the  profile  data  into  hioh-frequency  and  low-frequency 
parts.  In  order  to  characterize  more  fully  this  partition  of  en- 
ergy between  the  low-frequency  "qeostrophi c'  and  hinh-frequency 
'internal  wave"  parts  of  the  profiles,  we  have  averaoed  all  20  pro- 
files in  the  time  series.  (The  20  profile  time  series  occupies  a 
time  interval  of  about  103  hours.)  This  averaging  was  done  simply 
by  removina  the  "barotropic,"  or  depth-averaged,  part  of  each  pro- 
file, and  then  taking  an  average  over  the  20  drops  at  each  level. 

The  depth-averaaed  part  of  each  profile  was  removed  because  it  is 
unknown  and  is  not  the  true  barotropic  part  of  the  velocity  profile 
(see  Sanford  (1971)).  The  resulting  low-frequency  profiles  are 
shown  in  Figures  4 and  5 for  the  east  and  north  components. 

There  is  some  indication  in  these  mean  profiles  of  a decrease 
in  the  mean  shear  in  the  18°  water  (around  400  dbar).  There  is  also 
an  indication  that  the  mean  shear  increases  again  near  the  bottom, 
in  the  Antarctic  Bottom  Water.  Comparison  of  the  profile  of  average 
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horizontal  velocity  (Fioures  4 and  5)  and  the  average  profile  of 
Brunt-Vaisala  frequency  (Figure  8)  appears  to  show  that  there  is  a 
relation  between  them.  We  note  that  there  is  a relative  minimum  in 
the  average  Brunt-Vaisala  profile  at  about  400  dbar.  This  depth 
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Figure  4.  Average  profile  of  the  east  velocity  component 
during  the  five-day  time  series. 
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corresponds  to  the  depth  region  in  which  the  vertical  shear  in  the 
profile  of  mean  horizontal  velocity  is  seen  to  decrease.  Also,  with- 
in several  hundred  decibars  of  the  bottom  the  average  Brunt-VaisSIS 
frequency  starts  to  increase,  and  it  continues  to  increase  as  the 
bottom  is  approached.  There  is  a corresponding  increase  in  the 
vertical  shear  of  the  average  profile  of  horizontal  velocity  as  the 
bottom  is  approached.  Thus,  relatively  small  values  of  the  mean 
Brunt-Vaisala  frequency  appear  to  correspond  to  a decrease  of  the 
vertical  shear  in  the  profile  of  average  horizontal  current,  and 
vice  versa. 

c.)  High-frequency  velocity  components  in  the  upper  2.5  km  as  a 
function  of  time 

From  the  20  high-frequency  profiles  in  the  five-day  time 
series,  (that  is,  the  profiles  obtained  after  the  mean  velocities 
of  Figures  4 and  5 have  been  removed),  contour  plots  of  east  and 
north  velocity  components  have  been  made  (Figures  6 and  7).  These 
plots  show  how  various  structures  in  the  horizontal  velocity  field 
evolve  as  time  progresses.  This  evolution  appears  as  changes  in 
the  depth  and  size  of  individual  features.  Before  discussing  the 
significance  of  these  plots,  a few  comments  on  their  construction 
should  be  made. 

First,  the  distribution  of  profiles  in  time  is  not  uniform. 

In  particular,  there  are  relatively  large  qaps  in  time  between  2240 
and  2260,  and  between  2280  and  230D.  Between  these  two  pairs  of 
drops  the  profiler  had  to  be  used  for  inter-comparison  drops  with 
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other  vertical  profilers,  and  these  inter-comparison  drops  were  not 
done  at  the  central  mooring.  Hence,  gaps  in  the  time  series  re- 
sulted. Also,  the  time  interval  between  profiles  in  the  second 
half  of  the  series  (after  230D)  is  generally  less  than  it  is  for 
orofiles  in  the  first  half.  This  uneven  time  distribution  allows 
some  regions  of  the  plots  to  be  contoured  more  easily  than  others. 

Second,  only  the  0 cm/sec  contour  line  is  shown.  Since  the 
main  intent  in  drawing  the  contour  plots  was  to  attempt  to  see 
vertical  phase  propagation  of  the  near-inertial  waves,  a more  de- 
tailed contouring  was  not  required.  Regions  with  negative  velocity 
components  are  indicated  by  oblique  hatching  in  Figures  6 and  7. 

Finally,  features  for  which  contouring  was  ambiguous  are 
denoted  by  dashed-line  contours.  There  are  two  main  reasons  for 
this  ambiguity.  One  is  that  as  the  depth  increases,  the  overall 
energy  level  of  the  profiles  decreases  (see  Figure  8).  Thus, 
at  depths  of  2000  dbar  to  2500  dbar,  the  true  behavior  of  the 
zero-crossinq  lines  becomes  unclear.  Below  2500  dbar  it  is  very 
difficult  to  draw  contours  with  any  confidence.  The  other  reason 
is  that,  even  at  shallow  depths,  certain  profiles  have  regions 
where  a velocity  component  is  close  to  zero  without  being  clearly 
positive  or  negative.  In  cases  where  it  appeared  that  structures 
on  either  side  of  this  ambiguous  region  were  clearly  negative, 
dashed-line  contours  were  used  to  join  these  negative-velocity- 
component  features  together.  In  other  cases,  there  are  regions  on 
one  drop  which  are  clearly  positive,  while  adjacent  drops  have 
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negative  structures  at  the  same  depth.  Such  a case  occurs  at  500 
dbar,  around  drop  2420  in  the  contour  plot  for  the  north  component. 
In  these  cases,  no  attempt  was  made  to  join  adjacent  negatlve-vel- 
ocity-component  features.  Taking  note  of  these  difficulties,  we 
can  now  point  out  some  of  the  significant  features  in  these  plots. 

Perhaps  the  most  striking  characteristic  is  the  tendency  for 
regions  of  negative  velocity  to  move  upward  in  time.  If  we  assume 
that  we  can  connect  contours  through  regions  where  a velocity  com- 
ponent is  not  clearly  positive  or  negative,  then  most  of  the  struc- 
tures with  negative  velocity  can  be  followed  across  much  of  the  103 
hours  over  which  the  time  series  lasted.  If,  on  the  other  hand,  we 
assume  that  this  connecting  of  structures  is  not  valid,  these  nega- 
tive-velocity regions  will  be  divided  into  a series  of  "hot-dogs," 


with  gaps  of  weakly  positive  velocity  in  between.  In  either  case, 
the  upward  motion  is  evident. 

In  a few  cases,  features  do  appear  to  move  downward  in  time. 
An  example  occurs  at  about  1000  dbar  in  the  north  plot,  around  drop 
228D.  It  should  be  noted  that  this  occurs  in  a region  where  the 
profiles  are  widely  spaced  in  time,  and  the  contouring  there  is  not 
clear. 

Several  other  points  can  be  made  concerning  these  two  fig- 


ures. First,  if  we  follow  the  negative-velocity-component  struc- 
tures as  they  move  upward  in  time,  we  see  that,  at  a given  depth. 
It  usually  takes  about  a day  for  two  successive  negative-velocity 
zones  to  pass  through  a given  depth.  This  indicates  that  these 
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structures  have  a period  of  about  a day,  which  is  near  the  local 
inertial  period  at  this  latitude. 

Second,  if  Figures  6 and  7 are  overlaid,  it  is  found  that, 
in  general,  negative-velocity-component  zones  in  the  east  plots 
occur  at  depths  which  are  somewhat  greater  than  the  depths  for  the 
corresponding  zones  in  the  north  plots.  In  other  words,  the  east 
and  north  velocity  components  are  not  in-phase  with  respect  to 
depth. 

Third,  the  speed  with  which  the  negati ve-velocity-component 
zones  move  upward  is  about  2 mm/sec  at  500  dbar,  and  this  speed  in- 
creases to  about  6 mm/sec  at  2000  dbar.  We  will  consider  these  ob- 
servations in  more  detail  in  later  chapters. 

d.)  The  average  perturbation  kinetic  energy  profile,  and  a smoothed 
profile  of  Brunt-Vaisala  frequency 

The  20  high-frequency  profiles  (with  vertical  averages  and 
the  low-frequency  or  "geostrophic"  shear  flow  removed)  were  then 
used  to  calculate  profiles  of  high-frequency,  or  "perturbation," 
kinetic  energy.  This  was  done  for  each  profile  simply  by  squaring 
the  horizontal  velocity  values  at  each  depth.  These  horizontal 
kinetic  energy  profiles  were  then  averaged,  at  each  depth,  over  the 
20  drops  and  a mean  kinetic  energy  profile  resulted.  This  profile 
is  shown  in  Fiqure  8.  To  simplify  plotting,  this  mean,  high-fre- 
quency, horizontal  kinetic  energy  profile  was  averaged  in  the  ver- 
tical by  taking  centered  averages  50  dbar  long  around  the  points 
50  dbar,  100  dbar,  150  dbar,  and  so  on.  Plotted  along  with  the 
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Figure  8.  Average  profiles  of  high-frequency  horizontal 
kinetic  energy  and  Brunt-Vaisala  frequency. 
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mean  horizontal  perturbation  kinetic  energy  profile  of  Figure  8 is 
a smoothed  Brunt-Vaisala  frequency  profile,  obtained  from  CTD  (con- 
ductlvity-temperature-depth)  stations  made  at  roughly  the  same  time 
and  location  as  the  time  series.  Because  of  difficulties  with 
measuring  velocities  by  CMVP  near  the  surface  (to  depths  of  about 
70  dbar),  the  points  at  0 dbar  and  50  dbar,  showing  rather  high 
kinetic  energies,  are  probably  not  accurate. 

It  is  possible  to  calculate  an  error  bar  for  Figure  8.  We 
assume:  the  variance  of  the  noise  signal  is  0.25  cm2/sec2  for 

either  velocity  component  (Appendix  C);  noise  at  a given  depth  is 
uncorrelated  with  the  signal  and  is  uncorrelated  between  profiles, 
but  the  noise  of  points  in  the  50  dbar  averages  may  be  correlated 
as  may  the  noise  in  the  east  and  north  velocity  components.  The 
contribution  of  the  noise  (if  it  is  Gaussian)  will  then,  for  each 
point  in  Figure  8,  be  distributed  as  a Chi-square  random  variable 
with  (at  least)  19  degrees  of  freedom  (Jenkins  and  Watts,  1969). 

The  result  of  applying  the  Chi-square  distribution  law  is  that  the 
difference  between  the  measured  and  actual  kinetic  energy  may  be 
expected  to  be  between  .18  and  .43  erg/cm3  95%  of  the  time  (the 
measured  kinetic  energy  being  greater  than  the  actual  kinetic  en- 
ergy). This  error  estimate  has  not  been  plotted  in  Figure  8,  since 
it  is  rather  difficult  to  see  (it  is  about  three  times  the  width  of 

4 

the  line  depicting  the  kinetic  energy  profile).  We  emphasize  that 

«• 

this  error  estimate  only  takes  into  account  instrumental  noise.  It 

- " does  not  take  into  account  other  sources  of  error,  such  as 
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contamination  of  the  high-frequency  profiles  due  to  leakage  of 
high-frequency  energy  into  the  low-frequency  profiles  (Figures  4 
and  5) . 

We  see  from  Figure  8 that  the  profiles  of  mean  horizontal 
perturbation  kinetic  energy  and  mean  Brunt-Vais51S  frequency  are 
very  similar.  The  horizontal  kinetic  energy  decreases  from  the 
surface  down  to  about  400  dbar,  increases  in  the  main  thermocline 
to  about  1000  dbar,  and  then  steadily  decreases  down  to  about  4500 
dbar.  This  parallels  the  behavior  of  the  Brunt-Vaisala  frequency. 

It  is  interesting  to  note  that  while  the  Brunt-Vaisala  frequency 
slowly  increases  from  about  4500  dbar  to  the  bottom,  the  mean  ki- 
netic energy  profile  Increases  similarly  to  about  100-150  dbar  off 
the  bottom,  but  then  decreases  again.  This  near-bottom  decrease  of 
energy  may  be  indicative  of  the  presence  of  a bottom  boundary  layer. 
The  possible  existence  of  a bottom  boundary  layer  will  be  considered 
in  more  detail  in  Chapter  IV. 

We  recall  that  the  time  series  lasted  for  only  103  hours. 

This  means  that  the  mean  kinetic  energy  profile  was  obtained  by 
averaging  over  approximately  four  wave  cycles  for  the  inertial 
waves.  It  is  not  surprising,  then,  that  there  is  considerable 
structure  in  the  mean  kinetic  energy  profile.  Presumably,  if  the 
time  series  of  profiles  were  of  longer  duration,  the  average  per- 
turbation horizontal  kinetic  energy  profile  would  follow  the  mean 
Brunt-Vaisala  profile  more  closely,  particularly  in  the  main 
thermocl ine. 
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e.)  The  stretching  of  the  vertical  coordinate  and  the  normaliza- 
tion of  velocity  by  the  mean  Brunt-Vaisala  profile 

Figures  6 and  7 suggest  that  the  high-frequency  motions  ob- 
served are  dominated  by  internal  waves  of  near-inertial  period 
which  are  moving  vertically  through  the  water  column.  An  important 
quantity  which  could  be  calculated  from  the  high-frequency  profiles 
is  the  horizontal  kinetic  energy  spectrum  as  a function  of  vertical 
wave  number.  However,  the  data  on  which  we  would  like  to  perform 
these  calculations  (that  data  being  the  ensemble  of  profiles  in  the 
time  series)  is  clearly  a nonhomogeneous  process.  Figure  8 shows 
that  the  overall  variance  (or  horizontal  perturbation  kinetic  en- 
ergy) of  the  high-frequency  profiles  is  a function  of  depth.  Also, 
an  examination  of  Figure  3 shows  that  the  most  energetic  wave s have 
vertical  length  scales  which  appear  to  be  shorter  in  the  thermo- 
cline  than  in  the  deep  water.  This  means  that  the  statistical 
properties  of  the  wave  field  are  functions  of  the  Brunt-Vaisala 
frequency,  which  in  turn  is  a function  of  pressure.  Some  method 
must  be  found  which  will  at  least  partially  convert  the  process 
represented  by  the  ensemble  of  profiles  into  a homogeneous  one. 
Figure  8 gives  an  idea  of  the  method  to  be  used.  Since  the  hori- 
zontal kinetic  energy  tends  to  follow  the  mean  Brunt-Vaisala  pro- 
file, vie  could  make  the  variance  more  uniform  with  depth  by  normal- 
izing the  horizontal  velocity  at  a given  depth  with  the  square  root 
of  the  Brunt-Vaisala  frequency  at  that  depth.  But  this  type  of 
normalization  is  what  would  be  expected  if  the  waves  were  obeying  a 
WKB  approximation  (Phillips,  1966),  since  in  this  approximation  the 
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horizontal  velocity  components  of  an  internal  wave  are  proportional 

1/2 

to  (N)  , where  N is  the  local  Brunt-Vaisala  frequency.  There- 

fore, we  chose  to  normalize  the  hiqh-frequency  velocity  profiles 


usino: 


u*(z)  = u(z)/(N(z)/No)1/2 


where  u is  an  original  velocity  component,  u is  a normalized  ve- 
locity component  (botli  functions  of  pressure,  z),  N(z)  is  the 
Brunt-VSisSlff  frequency  at  pressure  z and  N is  a reference  Brunt- 
V<fisffl£  frequency,  11  = 3 cph.  N has  been  given  this  value  to 
facilitate  comparison  with  theoretical  vertical  wave  number  spec- 
tra presented  later  (see  Chapter  V). 

We  have  also  pointed  out  that  the  scale  of  the  dominant 
waves  appears  to  change,  being  shorter  in  the  main  thermocline 
than  in  the  deep  water.  This  is  also  in  accord  with  the  WKB  ap- 
proximation, which  predicts  that  the  vertical  wave  length  of  an 
internal  wave  should  be  shorter  in  regions  of  greater  N.  Thus, 
it  is  also  possible  to  use  the  WKB  approximation  to  normalize  wave 
lengths  in  the  vertical.  This  has  been  done  by  stretching  the 
pressure  coordinate  according  to  the  following  differential  law: 


- ( Mi)-  > 

"o 


where  z is  the  original  pressure  coordinate,  and  z is  the  stretched 


pressure  coordinate. 


A 


38 


Fiqure  9 qives  z as  a function  of  z,  where  z is  in  decibars 
and  z*  is  in  "stretched  decibars."  Below  about  2000  dbar,  z(z*)  is 
approximately  a straioht  line.  This,  along  with  the  fact  that  the 
high-frequency  perturbation  kinetic  energy  is  roughly  constant  with 
depth  below  2000  dbar,  indicates  that  vertical  wave  number  spectra 
of  the  profiles  below  2000  dbar  could  be  approximately  calculated 
without  use  of  the  above  normal ization  and  stretching  procedures. 

But  to  include  the  thermocline  in  these  calculations  requires  the 
use  of  the  above  procedures. 

Figures  10  and  11  show  an  example  of  a profile  which  has  had 
the  above  procedures  (equations  1 1 - 1 and  1 1 -2 ) applied  to  it.  Fig- 
ures 10  and  11(a)  show  the  velocity  profile  219D  as  it  was  origi- 
nally obtained.  Figures  10  and  11(b)  show  the  same  profile  after 
the  "geostrophic"  shear  flow  has  been  removed.  Finally,  Figures  10 
and  11(c)  show  2190  after  the  above  procedures  for  stretching  and 
normalizing  the  profile  have  been  used.  With  the  exception  of 
values  in  the  top  100-200  stretched  decibars,  these  profiles  pre- 
sent a much  more  uniform  appearance  than  do  the  profiles  of  the 
east  and  north  components  of  219D  before  stretching  and  normalizing. 
The  departure  from  uniformity  in  the  top  200  stretched  decibars 
could  be  explained  either  by  the  fact  that  the  WKB  approximation 
breaks  down  near  the  surface,  where  N is  changing  rapidly,  or  by 
the  fact  that  the  shear  of  the  mean  horizontal  velocity  increases 
near  the  surface.  This  shear  is  not  taken  into  account  by  the  above 


WKB  normal ization. 
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Figure  10.  East  component  of  profile  2190  as  originally  obtained  (a),  after  removal  of  the 
mean  profile  (b),  and  after  stretching  and  normalizing  (c). 
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Figure  11.  North  component  of  profile  21 9D  as  originally  obtained  (a),  after  removal  of  the 
mean  profile  (b),  and  after  stretching  and  normalizing  (c). 


f.)  Spectral  decomposition  of  stretched  profiles 

The  time  series  of  profiles  is  now  in  a form  (the  stretched 
and  normalized  profiles,  such  as  2190  in  Fiqures  10  and  11(c))  in 
which  the  vertical  wave  number  spectrum  of  horizontal  high-frequen- 
cy kinetic  energy  can  be  calculated.  We  now  describe  and  present 
several  different  types  of  spectra  which  have  been  obtained. 

The  first  type,  shown  in  Figure  12,  is  the  simple  autospec- 
trum, or  the  total  energy  of  the  east  and  north  components  as  a 
function  of  stretched  vertical  wave  number.  For  this  spectrum,  a 
subset  of  nine  down  profiles  out  of  the  original  20  was  selected. 

The  only  reason  for  this  selection  was  that  some  profiles  have 
better  data  near  the  surface  than  others.  The  EMVP  quite  often 
does  not  obtain  accurate  velocity  values  for  the  first  50-70  dbar 
below  the  surface.  Since  the  WKB  procedure  described  above  tends 
to  weight  values  near  the  surface  more  heavily  than  values  in 
deeper  water,  it  was  important  to  choose  those  stretched  profiles 
which  have  the  best  data  near  the  surface.  Thus,  the  best  nine 
down  profiles  were  chosen  for  the  calculations.  The  top  50  dbar  in 
the  stretched  profiles  were  also  ignored  when  calculating  the  spec- 
trum. This  further  reduced  the  influence  of  near-surface  measure- 
ment errors.  The  95%  confidence  limits  were  calculated  using  an 
assumed  36  degrees  of  freedom  (Jenkins  and  Watts,  1968).  Each 
final  spectral  estimate  is  formed  by  averaging  the  original  esti- 
mates over  four  adjacent  wave  number  bands  in  each  profile,  and  then 
averaging  over  the  nine  profiles.  This  would  give  72  degrees  of 
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freedom  to  each  estimate,  if  east  and  north  velocities  are  statis- 
tically independent.  However,  Figure  13  indicates  that  the  east 
and  north  velocity  components  are  not  independent.  Therefore,  36 
degrees  of  freedom  appears  to  be  more  realistic. 

The  somewhat  unusual  units  in  the  spectrum  of  Figure  13  arise 
because  of  the  normalization  by  N(z).  Thus  NCM/S  means  "normalized 
cm/sec,"  and  C/SDB  means  "cycles  per  stretched  decibar."  If  m and 
m'  are  vertical  wave  numbers  in  the  unstretched  and  stretched  coor- 
dinates , respectively,  then  the  relation  between  them  is 


Normalized  cm/sec  are  the  units  obtained  when  equation  1 1 - 1 is  ap- 
plied to  an  original  velocity  profile. 

The  most  interesting  characteristic  of  the  autospectrum  is 
that  at  high  wave  numbers  it  has  a slope  of  about  -2.5  on  a loq-loq 
plot.  At  smaller  vertical  wave  numbers  this  slope  decreases.  This 
spectral  shape  will  be  compared  in  Chapter  V with  a theoretical 
vertical  wave  number  spectrum  derived  by  Garrett  and  Munk  (1972, 
1975). 

The  second  type  of  spectrum  which  is  presented  is  a form 
which  treats  the  horizontal  velocity  vector  (u,v)  as  a complex  vec- 
tor, u + iv,  (Gonnella,  1972;  Kooers,  1973).  At  any  vertical  wave 
number,  m',  the  Fourier  transform  of  u and  v separately  gives  sinu- 
soids for  each  component.  These  two  sinusoids,  taken  together, 
will  form  an  ellipse  in  the  u,  v plane,  the  position  of  the  velocity 
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Figure  13.  Clockwise  (CW)  and  counterclockwise  (ACW)  spectra  for  9 

down  profiles.  In  the  text,  CW  is  C(m’)  and  ACW  is  A(m'), 
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vector  on  the  ellipse  being  a function  of  z*,  the  pressure  In 
stretched  decibars  (increasing  downward).  Since  any  ellipse  can 
be  represented  as  the  sum  of  two  complex  vectors  rotating  In  op- 
posite directions,  but  with  the  same  frequency,  the  complex  vector 

u , + 1v  , can  be  represented  by 
m m 
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♦ iv 
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1 1-4 


where  m'  is  always  considered  to  be  positive,  and  u+  and  u are 
complex.  The  energy  spectrum  is  then  given  as  a decomposition  be- 
tween the  clockwise  energy,  C(m'),  and  the  counterclockwise  energy, 
A(m'),  defined  by: 


C(m')  * 1/2  < u_  . u_*>  ; 

II -5 

A(m')  * 1/2  < u+  - u+*>  , 

II -6 

where  the  brackets  denote  an  average  over  realizations  (In  this 
case,  profiles)  and  possibly  over  m',  the  stretched  vertical  wave 
number.  The  total  energy,  T(m'),  is  given  by 

T(m')  * A(rn')  + C(m').  II-7 

T(m')  is  plotted  in  Figure  12,  and  A(m')  and  C(m')  for  the  same 
nine  drops  are  plotted  in  Figure  13. 

Probably  the  most  striking  aspect  of  Figure  13  is  the  fact 
that  clockwise  energy  is  greater  than  counterclockwise  energy  over 
much  of  the  stretched  vertical  wave  number  range.  Only  in  the  high 
wave  number  region  of  the  spectrum,  where  the  energy  is  relatively 
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low,  do  we  find  points  where  counterclockwise  is  greater  than 
clockwise  energy.  This  means  that  the  waves  are  el  1 i pti cal ly 
polarized  in  the  clockwise  sense,  or  in  other  words,  that  there  is 
a distinct  tendency  for  the  horizontal  current  vector  to  rotate  in 
a clockwise  sense  with  depth,  as  seen  by  an  observer  looking  down 
from  above.  This  has  important  implications  for  the  theory  of 
internal  waves,  which  will  be  discussed  in  the  next  chapter. 

The  third  type  of  spectrum  that  we  present  here  is  simply 
that  obtained  by  plotting  the  energy  of  the  east  and  north  compo- 
nents separately.  This  spectrum  is  given  in  Figure  14  for  the 
same  nine  down  drops  as  in  Figure  12.  We  see  that,  in  contrast  to 
the  difference  in  Figure  13  between  the  clockwise  and  counterclock- 
wise spectra.  Figure  14  shows  that  there  is  almost  no  difference 
between  the  energy  of  the  east  and  north  velocity  components.  This 
is  an  indication  of  isotropy;  that  is,  that  in  the  mean  there  is  as 
much  energy  in  the  east-west  as  in  the  north-south  direction. 

Finally,  we  present  in  this  section  results  of  spectral  cal- 
culations performed  on  nine  up  profiles  from  the  same  time  series. 
It  often  happens  that  on  a single  drop  (up  and  down)  the  down  part 
will  have  good  near-surface  data  while  the  up  part  will  not.  For 
this  reason,  the  nine  best  up  profiles  chosen  do  not  coincide  with 
the  nine  best  down  profiles.  These  spectra  are  shown  in  Figures 
15,  16  and  17.  Each  spectral  estimate  represents  an  average  over 
four  adjacent  wave  number  bands  and  over  nine  up  profiles.  The 
major  reason  for  performing  the  same  calculations  for  up  and  down 
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Figure  14.  Spectra  of  east  and  north  velocity 
components  for  9 down  profiles. 


KINETIC  ENERGY  DENSITY 
(NCM/S)  2 /(C/SDB) 


49 


WAVE  LENGTH 
(SOB) 


VERTICAL  WAVE  NUMBER 

(C/SDB) 


Figure  15.  Total  spectrum  (east  and  north  energy)  for  9 up  profiles. 
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Figure  17.  Spectra  of  east  and  north  velocity 
components  for  9 up  profiles. 


profiles  was  to  show  that  there  is  very  little  difference  between 
them  (even  though  all  of  the  same  profiles  were  not  used  in  both 
calculations).  That  the  two  computations  are  very  similar  can  be 
seen  from  comparing  Figures  1 2 and  15,  13  and  16,  and  14  and  17. 

The  95%  confidence  limits  (with  36  degrees  of  freedom)  have 
been  included  for  the  total  vertical  wave  number  spectra  (Figures 
12  and  15)  and  for  the  spectra  of  east  and  north  velocity  compo- 
nents (Figures  14  and  17).  These  confidence  limits  should  be  used 
with  caution,  however,  since  they  strictly  apply  only  to  a white 
noise  spectrum.  The  spectra  presented  in  this  section  are  clearly 
"red."  Fven  though  the  original  profiles  were  pre-whitened  (by 
taking  first  differences)  before  applying  the  Fourier  transform, 
and  the  above  spectra  were  then  produced  by  re-coloring  the  spectra 
obtained  from  the  first-difference  profiles,  it  is  likely  that  the 
pre-whitening  did  not  completely  convert  the  profiles  into  ones 
with  flat  spectra.  This  means  that  the  actual  number  of  degrees 
of  freedom  could  be  less  than  36.  VJe  have  also  not  plotted  confi- 
dence limits  for  the  clockwise  and  counterclockwise  spectra. 

The  clockwise  and  counterclockwise  energy  estimates  depend  not  only 
on  the  probability  distribution  of  the  autospectra  for  the  east  and 
north  velocity  components,  but  also  on  the  probability  distribution 
of  the  quadrature  spectrum  between  east  and  north  velocity  compo- 
nents (Appendix  D).  Instead  of  deriving  a theoretical  estimate  for 
the  probability  distribution  of  the  clockwise  and  counterclockwise 
spectra,  it  is  probably  more  important  to  point  out  that,  as 
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indicated  below,  clockwise  and  counterclockwise  spectra  calculated 
from  four-drop  groups  in  the  time  series  always  showed  that  clock- 
wise energy  was  greater  than  counterclockwise  energy  over  most  of 
the  stretched  vertical  wave  number  range.  Furthermore,  while  av- 
eraging the  spectra  of  two  four-drop  groups  together  gave  clockwise 
and  counterclockwise  spectra  very  similar  to  those  of  Figures  13 
or  16,  averaging  three  four-drop  groups  together  did  not  appreciably 
change  the  spectra.  This  indicates  that  nine  profiles  are  enough 
to  give  stable  estimates  of  the  clockwise  and  counterclockwise 
spectra. 

Before  the  spectrum  of  Figure  12  was  calculated,  five  pre- 
liminary autospectra  were  obtained  by  using  four-drop  groups  in 
the  series  of  20  down  profiles.  That  is,  drops  219D-223D  were  used 
to  obtain  a spectrum,  then  drops  224D-228D,  and  so  on.  All  four- 
drop  groups  showed  a dominance  of  clockwise  over  counterclockwise 
energy.  Each  of  the  five  preliminary  spectra  also  showed  a tend- 
ency for  the  energy  to  peak  at  vertical  wave  lengths  of  100-130 
sdb.  Drops  230D-233D,  in  particular,  showed  a strong  energy  peak 
there.  However,  in  some  cases  (for  example,  219D-223D)  this  peak 
did  not  show  up  as  clearly.  The  fact  that  the  energy  spectra  did 
appear  to  have  some  dependence  on  time  indicates  that  the  spectra 
calculated  from  the  nine  best  profiles  out  of  the  20  will  tend  to 
smooth  out  peaks  in  the  energy  spectra. 

On  the  other  hand,  it  is  clear  that  energy  contributed  at 
vertical  scales  of  100-130  sdb  is  an  important  part  of  the  stretched 


profiles.  From  Figure  3 we  see  that  the  most  apparent  visual  var- 
iability in  this  profile  occurs  in  the  above  range  of  vertical 
scales.  (In  the  original  pressure  coordinate,  the  above  scales 
correspond  to  100-150  dbar  in  the  main  thermocline,  and  reach 
1000  dbar  in  the  deep  water.)  A comparison  of  drops  21 90  and  221 D 
in  the  stretched  vertical  coordinate  also  shows  that  much  of  the 
"mirror-imaging"  that  occurs  between  these  profiles  is  caused  by 
waves  with  the  above  vertical  length  scales. 


g. ) Paired  profiles 

During  the  course  of  the  time  series  experiment,  several 
simultaneous  paired  profiles  were  made,  using  two  EMVP's.  These 
were  done  at  separations  of  from  100  meters  to  about  15  kilometers. 
Figure  18  shows  two  profiles  which  were  made  simultaneously  at  a 
separation  of  4.8  kilometers.  Profile  2351)  was  made  to  the  south 
of  profile  236U.  It  can  be  seen  from  this  figure  that  many  velocity 
features  of  235U  can  be  traced  visually  to  correspondino  features 
in  profile  236U,  although  there  is  quite  often  a depth  chanqe  be- 
tween the  more  energetic  features  of  profiles  235U  and  236U.  This 
indicates  that:  a.)  many  of  the  features  observed  in  the  two 

series  (at  least  these  with  wave  lengths  greater  than  about  100 
dbar)  have  a large  ratio  of  horizontal  to  vertical  length  scales; 
and  b.)  the  energetic  features  in  the  two  profiles  are  not  hori- 
zontal, but  are  slightly  inclined.  We  will  return  later  to  a con- 
sideration of  the  inclination  of  velocity  features  between  2351)  and 


VELOCITY  PROFILES  235U  a 236U 

SIMULTANEOUS  DROPS  48  KM  APART 


(SH313W)  Hid 30 


Figure  i8.  Simultaneous  paired  drops,  235U  and  236U,  separated 
horizontally  by  4.8  km. 
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Chapter  III  Vertical  Polarization  and  the  Vertical  Propagation  of 
Internal -inertial  Waves 

a.)  Summary  of  data 

Before  going  further,  it  will  be  useful  to  summarize  the 
results  of  the  previous  chapter.  Figure  3 indicates  that  high- 
frequency  energy  in  the  profiles  is  dominated  by  waves  near  the 
local  inertial  frequency.  There  is  probably  also  energy  contri- 
buted by  diurnal  and  semidiurnal  tides.  Figures  6 and  7 show  that 
the  velocity  structures  associated  with  these  inertial  waves 
clearly  move  upward  with  time  at  the  time  series  location.  These 
two  figures  also  show  that  the  east  and  north  components  of  hori- 
zontal velocity  are  usually  not  in-phase  with  depth,  but  that  the 
east  component  laqs  behind  the  north  component.  This  indicates 
that  these  near-inertial  waves  are  polarized  with  depth,  a fact 
which  is  born  out  by  the  clockwise  and  counterclockwise  spectra  of 
Figures  13  and  16.  These  spectra  show  that  the  waves  are  ellipti- 
cal 1 y polarized  with  depth,  and  they  show  that  this  polarization 
is  predominantly  clockwise  (the  clockwise  energy  is  greater  than 
the  counterclockwise  energy  over  most  of  the  observed  wave  number 
band).  We  note  that  the  observed  lag  in  Figures  6 and  7 between 
the  east  and  north  components  also  indicates  that  these  waves  are 
polarized  clockwise  with  increasing  depth.  Finally,  Figures  14  and 
17  indicate  that  there  is  very  little  difference  between  the  eneroy 
of  the  east  component  and  that  of  the  north  component.  We  would  now 
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like  to  see  how  the  above  observations  may  be  related  to  internal 
wave  theory. 


b.)  The  vertical  spatial  behavior  of  horizontal  internal  wave 
velocity  components 

Expressions  for  the  east  and  north,  or  u and  v velocity 
components  of  a single  propagating  internal  wave  are  given  in  the 
WKB  approximation  by: 

G = ( + cok  + iflH  x A (N»(z)  - u>2)  V4 
v « ( + ikf  * wi)  J M(ua  . f2)l/2(k2  + ! 2 ) 1 /2 
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where  A is  the  wave  amplitude,  and  the  coordinate  system  is  chosen 
so  that  x is  positive  eastward,  y is  positive  northward  and  z is 
positive  upward.  The  coordinates  2 and  2*  will  be  reserved  for  the 
vertical  pressure  coordinate  (increasing  downward)  in  decibars  and 
stretched  decibars,  respectively.  (The  reader  may  refer  to  Chapter 
V for  the  derivation  of  these  formulas.)  In  order  to  simplify 
these  expressions,  we  rotate  the  coordinate  system  through  an  angle 
a as  shown  in  Figure  19.  In  the  new  coordinate  system: 
u'  = u cosa  + v sina  ; 
v'  = -u  sina  + v cosa  ; 


k ' = k cosa  + 1 sina  ; 

1 ' = -k  sina  + 1 cosa  = 0 , 


1 1 1 -2 


~u  k y. 


•+  m ~ nr~\ 


Figure  20.  Direction  of  rotation  of  the  horizontal  velocity  vector  of  an 
internal  wave  for  positive  and  negative  vertical  wave  numbers 
The  curved  arrow  shows  the  direction  in  which  7 rotates  with 
increasing  z.  ~ 
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where 


cosa  = 


/k2  + l2 


sina  = 


/k2  + l2 


1 1 1 -3 


In  the  new  coordinate  system,  we  have 


u'  = + 1 ^ 

v-  - + a 3 
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We  let 


k'2  1/2  f 2 2 1/2 

J (N2(z)  - ^2)  dz  : m(z ) 


where  m(z),  the  local  vertical  wave  number,  is  given  approximately 


m(z)  = (-r^rr)172  (n2(z)  - u>2) 


V - f2 


1 1 1-6 


Mote  that  m(z)  is  a positive  quantity.  The  choice  of  siqns  in  the 
expressions  for  horizontal  velocity  determines  whether  the  wave  is 
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propagation  upward  or  downward.  Then 
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The  horizontal  component  of  the  wave  number  vector,  k',  is  assumed 

to  be  in  the  direction  of  horizontal  phase  propagation,  so  u > 0. 

1 /2 

We  note  in  passing  that  u'  and  v'  are  proportional  to  N , while 
m(z)  is  proportional  to  N.  This  confirms  the  connection  made  in 
Chapter  II,  without  proof,  between  the  WKB  approximation  and  the 
procedures  used  for  stretching  and  normalizing  the  velocity  pro- 
files. Taking  real  parts  of  u'  and  v': 


Re  ( u ' ) = u ' 
Re  ( v * ) = v’ 


+ cos  (kV  + mz  - wt)  ")  ^ (n2(z)  - U)2)1/4 

+ — • sin  ( k ' x 1 + mz  - wt)  3 (w2  - f2)^ 

II 1-8 


The  upper  signs  are  chosen  if  m is  directed  upward,  and  the  lower 
signs  are  chosen  if  m is  directed  downward.  A is  assumed  to  be 
real.  To  observe  the  behavior  of  u'  and  v'  as  z varies,  we  set 
x'  = t = 0.  Then 


u‘  = + cos  (+  mz)  Y x ^ (N2(z)  - to2)  1/4 

v'  = + — sin  (+  mz)  f (w*  - f2)]/2  . 1 1 1 -9 

(0 


The  two  cases  for  the  choice  of  upper  or  lower  sians  are  shown  in 
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Figure  20.  We  note  that  i'  = v'/(N2(z)  - w2)  . The  horizontal 

vector,  v' , is  plotted  in  Figure  20  instead  of  v'  in  order  to  re- 
move variations  in  the  horizontal  velocity  vector  magnitude  brought 
about  by  changes  in  N(z).  The  ellipses  represent  hodographs  of  the 
horizontal  velocity  vector  as  z varies.  As  w -*  f,  the  wave  becomes 
circularly  polarized.  We  note  that  although  we  have  been  writing 
simply  A for  the  wave  amplitude  in  the  above  expressions,  the  wave 
amplitude  is  more  generally  a function  nr  vertical  wave  number, 
frequency,  and  a,  the  direction  of  horizontal  phase  propagation; 
that  is,  A = A(m,  u,  a).  In  the  + m case,  the  horizontal  velocity 
vector  rotates  counterclockwise  with  increasing  z,  while  in  the  - m 
case,  the  sense  of  rotation  is  the  opposite.  The  result  is  that 
there  is  a unique  relation  between  the  sign  of  m and  the  sense  of 
polarization  in  the  internal  wave  hodographs  of  Fiqure  20. 

We  can  now  relate  the  sign  of  m(z)  to  the  direction  of  the 
vertical  group  velocity  component  for  internal  waves.  Fiqu-e  21 
shows  a propagation  diagram  for  internal  waves.  The  vector  It  is 
the  wave  number  vector  for  the  wave  and  £g  is  the  group  velocity 
vector  (or  the  direction  of  energy  propagation). 

For  simplicity,  only  two  wave  number  vectors  have  been 
shown,  one  pointing  toward  the  sea  surface  (toward  positive  z)  and 
the  other  pointing  away  from  the  sea  surface:  the  dispersion  rela- 

tion gives 


u>2  = N2  sin2  (0)  + f2  cos 2 (0 ) , 
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where  0 is  the  angle  that  the  wave  number  vector,  k,  makes  with 
the  z axis.  For  w close  to  f,  o is  small,  and  £ is  almost  verti- 
cal. For  constant  N,  t is  also  independent  of  depth.  On  the  other 
hand,  if  N varies  with  depth,  then  0 and  the  vertical  component,  m, 
of  k also  change  with  depth.  If  the  wave  propagation  may  be  approx- 
imated by  a WK8  expansion,  however,  equations  1 1 1 - 1 0 and  1 1 1 - 11  re- 
main true  locally,  except  that  m and  N are  now  functions  of  z.  The 
magnitude  of  the  group  velocity  is  given  by: 

Cg  = — L_  -&*1  = (N2  - f2)  sin  0 cos  0 

Ik  I 30  |t|  [N2sin20  + f2cos~2  0]1/2  . II 1-11 

The  group  velocity  vector  is  perpendicular  to  £ and  points  away 
from  the  z axis  if  N > f,  which  it  was  everywhere  in  the  water 
column  during  the  time  series  experiment.  The  vector  v in  Figure 
21  is  the  water  velocity  due  to  the  wave.  Since  the  water  is  as- 
sumed to  he  incompressible,  £ • v = 0 always.  There  is  no  depen- 
dence  on  horizontal  direction  in  the  above  formulas,  so  the  k and 
fg  vectors  can  be  rotated  about  the  z axis  without  chanae  of  fre- 
quency or  magnitude  of  Co.  We  have  also  drawn  two  helices  in 
pinure  21.  These  helices  represent  the  path  that  the  tip  of  the 

velocity  vector,  v,  makes  as  distance  from  the  origin  is  increased 
-► 

along  k. 

Lookina  down  on  this  figure  from  above  the  z axis,  it  is 
readily  seen  that  the  sense  of  polarization  of  the  wave  (that  is, 
the  direction  in  which  v rotates  with  increasing  depth)  changes, 
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dependinq  on  whether  k is  directed  toward  or  away  from  the  sea  sur- 
face.  For  a vector  k directed  toward  the  sea  surface,  the  wave  is 
polarized  clockwise  with  increasing  depth,  while  for  a k vector 
directed  away  from  the  sea  surface,  the  wave  is  polarized  counter- 
clockwise with  increasing  depth.  This  is  simply  a restatement  of 
the  result  expressed  by  equations  1 1 1 -8  and  1 1 1 -9 . 

We  point  out  that  polarization  can  be  inferred  easily  from 
the  behavior  of  internal  waves  in  time.  We  assume  we  are  observing 
at  a fixed  point  within  the  ocean  in  the  Northern  Hemisphere.  Then 
it  is  known  that  a velocity  vector  due  to  an  internal  wave  propa- 
gating through  that  point  will  rotate  in  a clockwise  sense  with 
time,  due  to  the  influence  of  the  Coriolis  aceeleration.  Since 
the  internal  wave  is  assumed  to  be  propagating  along  the  direction 
of  the  wave  number  vector,  ft,  of  Figure  21  (either  toward  or  away 
from  the  sea  surface),  the  waves  have  to  be  polarized  as  shown  in 
the  figure  in  order  to  give  a velocity  vector,  v,  which  rotates 
clockwise  in  time  at  a given  point  in  space. 

Now,  Figure  21  also  shows  that  if  £ is  pointed  toward  the 
*► 

sea  surface,  Cg  is  pointed  away  from  it,  and  vice  versa.  There- 
fore, if  either  the  sense  of  wave  polarization  or  the  direction  of 
the  vertical  phase  propagation  (the  sign  of  m)  can  be  determined, 
then  the  direction  of  the  vertical  component  of  the  group  velocity 
* vector,  and  thus  the  direction  of  vertical  energy  propagation,  can 

be  found  (Leaman  and  Sanford,  1975). 

V- 
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The  data  presented  in  Chapter  II  allows  us  to  determine  both 
the  direction  of  the  vertical  component  of  k (at  least  for  the 
dominant  waves)  and  the  sense  of  wave  polarization  in  the  vertical. 
As  already  pointed  out,  the  dominant  waves  in  Fioures  6 and  7 al- 
most always  appear  to  be  propaaatino  upward  through  the  upper  2.5 
kilometers  of  the  water  column.  This  indicates  that  the  direction 
of  phase  propagation,  or  k,  for  these  waves  is  pointed  toward  the 
sea  surface.  If  this  is  true,  the  above  discussion  shows  that  the 
waves  should  be  polarized  predominantly  clockwise  with  increasinq 
depth.  Both  the  clockwise  and  counterclockwise  spectra  of  Figures 
13  and  16,  and  the  fact  that  the  east  and  north  velocity  components 
are  out-of-phase  in  depth,  with  north  leading,  in  Figures  6 and  7, 
show  that  the  waves  are  in  fact  dominated  by  clockwise  polarization. 

The  Influence  of  the  mean  flow  on  the  observed  waves  has  so 
far  been  ioncred.  The  depth  dependent  part  of  the  mean  flow  pro- 
file has  already  been  presented  (Figures  4 and  5).  Another  depth 
independent  component  of  the  low-frequency  profile  may  also  exist, 
but  its  maqnitude  cannot  be  evaluated  from  EMVP  data  (Sanford, 

1971).  In  order  to  evaluate  a possible  depth  independent  part  of 
the  mean  profile,  we  have  calculated  the  averaoe  horizontal  cur- 
rent at  the  central  mooring  as  measured  by  a current  meter  at 
1513  dbar  (1496  meters).  The  average  horizontal  velocity  compo- 
nents over  the  interval  of  the  profiler  time  series,  as  determined 
from  the  current  meter  record,  are  + 1.6  cm/sec  to  the  west  and 
+ 2.5  cm/sec  to  the  north.  The  magnitudes  of  the  horizontal  current 
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components  at  1513  dbar  in  Figures  4 and  5 are  - 1.5  cm/sec  to  the 

west  and  ~ 1.5  cm/sec  to  the  north.  This  indicates  that  there  is 

very  little  shift  in  the  east  component,  while  ~ + 1 cm/sec  should  be 

added  to  the  north  profile  to  make  the  average  north  flow  observed  by 

the  current  meter  and  the  profiler  agree. 

-► 

Since  the  wave  number  vector,  k,  is  not  exactly  vertical  for 
waves  propagating  energy  in  the  vertical,  we  expect  that  the  advec- 
tion  of  waves  by  the  mean  flow  will  cause  an  apparent  vertical  propa- 
gation. The  amount  and  sign  of  this  apparent  propagation  (upward  or 
downward)  depends  on  the  direction  of  horizontal  wave  propagation 
relative  to  the  mean  flow.  Figures  4 and  5 show  that  the  major  part 
of  the  mean  flow  is  oriented  north-south.  Therefore,  if  the  waves 
were  propagating  northward  (and  upward)  the  advection  would  cause  an 
apparent  downward  phase  propagation  to  be  added  to  the  true  upward 
phase  propagation  of  the  wave  above  1500  dbar.  We  do  not  have  adequate 
information  on  the  directions  of  horizontal  phase  propagation  for  the 
observed  waves.  Some  work  (Hunk  and  Phillips,  1968)  could  be  cited  to 
argue  that  near-inertial  waves  would  preferentially  prooacate  north- 
ward. Their  work  actually  shows  that  standing  waves  would  exist  south 
of  the  critical  latitude  for  these  waves  (the  latitude  at  which  the 
wave  frequency  equals  the  local  inertial  frequency).  However,  if  en- 
ergy is  taken  out  of  the  northward-propagating  waves  as  they  reflect 
off  the  bottom  (as  will  be  argued  later)  standing  waves  will  not  be 
set  up,  and  northward  energy  propagation  would  be  expected  to  be  dom- 
inant south  of  the  critical  latitude.  If  this  is,  in  fact,  true, 
Figures  6 and  7 show  that  k must  be  directed  toward  the  sea  surface: 
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that  is,  the  actual  upward  phase  propagation  of  the  wave s is  enough 
to  overcome  any  apparent  downward  propaoation  caused  by  advection.  A 
strong  indication  that  the  observed  waves  are  propaaatinq  upward  can 
be  found  by  examining  Figures  4 and  5.  It  is  known  from  trackinn 
neutrally  buoyant  floats  that  the  mean,  or  "low-frequency1  flow  is 
small  at  depths  of  around  1500  dbar  (Voorhis,  1974).  Fioures  4 and  5 
also  show  that,  taking  into  account  the  shift  in  the  north  profile  re- 
quired to  make  the  mean  flow  as  observed  by  the  profiler  and  current 
meter  equal  at  1513  dbar,  the  mean  flow  in  the  depth  range  1400-1500 
dbar  is  small  (<  3 cm/sec).  Therefore,  there  should  be  little  advec- 
tion of  waves  at  this  depth.  However,  Figures  6 and  7 show  that  the 
waves  are  moving  upward  at  this  depth.  In  fact,  the  rate  of  upward 
motion  is  greater  at  1500  dbar  than  it  is  at  lesser  depths.  We  note 
that  the  sense  of  polarization  is  not  affected  by  the  advection  of 
waves  Dast  the  observation  point.  Therefore,  for  the  purpose  of 
determining  the  direction  of  vertical  wave  energy  flux,  observation 
of  polarization  appears  to  be  more  useful  than  observation  of  the 
apparent  vertical  phase  propagation. 

The  behavior  of  wave  polarization  and  the  direction  of 
phase  propagation  both  appear  to  indicate  that  the  vertical  energy 
flux  is  directed  downward,  away  from  the  sea  surface.  It  is  nec- 
essary, however,  to  discuss  just  what  the  partition  of  energy  be- 
tween clockwise  and  counterclockwise  means  for  these  waves.  In 
particular,  we  have  limited  the  above  discussion  to  a sinole  wave. 

It  is  more  appropriate  to  relate  the  partition  of  clockwise  and 
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counterclockwise  energy  to  the  energy  carried  by  upward  and  downward 
propagating  waves.  This  will  be  done  in  Chapter  IV. 

It  is  known  from  internal  wave  theory  (see,  for  example, 
Garrett  and  Munk,  1972)  that  internal  waves  with  frequencies  near 
the  inertial  frequency  are  isotropic  in  time;  that  is,  the  magni- 
tude of  the  frequency  spectrum  of  a near-inertial  wave  horizontal 
velocity  component  will  be  independent  of  the  direction  of  that 
velocity  component.  The  same  Isotropy  will  be  true  for  the  verti- 
cal wave  number  spectra  of  near-inertial  waves,  provided  only  that 
east  and  north  spectra  for  a sufficient  number  of  profiles  are 
averaged  together. 

To  see  why  this  is  true,  we  consider  a single  internal  wave 
propagating  energy  downward.  The  u'  and  v'  velocity  components  are 
given  by  equations  III -8 : 


u'  = - cos  (k'x'  + mz  - wt) 
d 

v'  = sin  (k'x'  + mz  - wt) 

d w 


x ft  - “Z)1/4 


(w2  - f2) 


1/2 


III-12 


where  A is  the  wave  amplitude.  We  point  out  that  the  hodograph  in 
z for  this  wave  is  given  by: 


A2 (m,w,a) 


(N2  - w2)1/2 
(w2  - f2) 


1 1 1 - 1 3 


At  any  time,  the  ratio  of  semi-major  to  semi-minor  axis  is  , and 
as  o)  -*•  f,  the  wave  becomes  circularly  polarized.  If  we  assume  the 
presence  of  a second  wave  propagating  upward  with  the  same  vertical 
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wave  number  magnitude,  m,  and  the  same  frequency,  w,  but  with  dif- 
ferent amplitude,  B,  then  the  velocity  components  for  this  wave  are 
given  by: 


u'  = cos  (k 1 x ' - mz  - wt) 
u 

v'  = — sin  (k'x‘  - mz  - wt) 

U U) 


x B 


Aiil 

f2) 


1/4 

1/2 


1 1 1-14 


We  assume  A and  B to  be  real  now,  althouqh  this  is  not  necessary. 
We  also  assume  we  are  observing  at  the  coordinate  origin,  so 
x'  = 0.  The  result  of  adding  the  two  waves  together  gives: 


w 


1/4 


f2) 


1/2 


[(B  - A)  cos(mz)cos(wt)  - (A  + B)  sin(mz)sin(wt)]. 


v' 

T 


4!i 

U2 


1/4 

1/2 


[(A  - B)  cos(mz)sin(wt)  - (A  + B)  sin(mz)cos(wt)]. 

1 1 1 - 1 5 


We  see  that,  for  a fixed  z,  the  magnitude  of  u'T  averaged 
over  a cycle  in  time,  becomes  eoual  to  the  average  of  v'  as  u>  -*•  f 
independently  of  the  magnitudes  of  A and  B.  Thus,  spectra  of  u^. 
and  v^.  in  frequency  may  be  expected  to  show  the  same  energy  when 
the  frequency  is  close  to  f. 

The  same  is  not  true,  however,  if  vertical  wave  number  spec- 
tra of  u^.  and  v^  are  considered,  as  can  be  seen  by  fixina  t in  the 
above  expressions.  If  A or  B is  zero,  the  vertical  wave  number 
spectra  for  u^  and  v^  will  be  approximately  equal  for  w close  to  f. 


71 


However,  if  A = B (which  would  occur  if  standing  waves  existed  in 
the  vertical),  energy  in  the  vertical  wave  number  spectra  of  u^  and 
v;  will  depend  on  t,  the  time  in  the  wave  cycle  at  which  the  spectra 
were  obtained.  For  example,  in  this  case,  at  t = 0,  the  spectrum 
of  V.J.  will  contain  all  the  energy  at  that  particular  vertical  wave 
number,  m,  and  the  spectrum  of  u^.  would,  in  theory,  be  zero. 

However,  if  vertical  wave  number  spectra  of  u^.  and  v_^_  are 
taken  at  enough  different  times,  t,  in  the  wave  cycle,  we  again  ap- 
proach isotropy  for  u>  -*•  f.  This  is  true  because,  at  a given  m,  the 
energy  in  u^.  and  v^  will  be  the  same  if  averaged  over  time.  The 
average  vertical  wave  number  spectra  of  Figures  14  and  17  can  be 
considered  as  the  sum  of  individual  wave  number  spectra  taken  at 
numerous  different  times,  t,  in  a wave  cycle.  Then  the  fact  that 
much  of  the  energy  in  the  profiles  appears  to  be  contributed  by 
waves  with  frequencies  close  to  f explains  why  there  is  an  equi- 
partition  of  energy  between  the  east  and  north  components. 

c.)  Polarization  in  the  deep  water 

As  mentioned  previously,  we  can  approximately  calculate  deep 
water  spectra  (below  2500  dbar)  without  stretching  and  normalizing 
the  velocity  profiles.  This  is  possible  because  N does  not  change 
much  in  the  deep  water.  Calculating  the  clockwise  and  counter- 
clockwise spectra  over  the  depth  range  2500-5500  dbar  for  several 
profiles  Indicates  that  dominance  of  clockwise  energy  over  counter- 
clockwise energy  exists  in  the  bottom  half  of  the  water  column,  at 
least  for  vertical  wave  lengths  greater  than  about  100  dbar.  This 


is  important  because  it  indicates  that,  even  in  deep  water,  the 
net  flux  of  energy  by  these  near-inertial  waves  is  downward. 

Pollard  (1970)  indicated  in  a theoretical  study  that  it  would  be 
difficult  for  near-inertial  waves  to  propagate  down  into  the  deep 
water  with  the  observed  intensity.  The  result  of  the  above  cal- 
culations indicates  that  these  waves  are,  nevertheless,  coming  down 
from  regions  farther  up  in  the  water  column.  The  Pollard  model 
assumes  that  inertial  waves  are  generated  at  the  surface  by  an  im- 
pulsive force  that  acts  over  some  time  and  space  interval.  In 
Chapter  V we  will  examine  the  behavior  of  the  vertical  energy  flux 
of  near-inertial  waves  for  an  assumed  form  of  the  internal  wave 
energy  spectrum  in  frequency  and  vertical  wave  number  space. 

d.)  Comparison  with  meteorological  observations 

As  was  pointed  out  in  the  introduction,  meteorologists  have 
been  conducting  experiments  to  obtain  vertical  profiles  of  atmos- 
pheric winds  for  many  years.  Their  results  may  be  compared  with 
results  obtained  by  the  Ef'VP  in  the  ocean. 

Various  authors  have  pointed  out  that  the  existence  of  in- 
ternal waves  in  the  atmosphere  has  not  been  proved.  For  example, 
Justus  and  Woodrum  (1973)  felt  it  necessary  to  state  that:  'There 

has  been  no  unambiguous  resolution  of  gravity  waves  in  the  upper 
atmosphere  (i.e.  simultaneous  observations  of  amplitudes,  phases, 
and  frequency  sufficient  to  verify  propagation  accordino  to  the 
theoretical  dispersion  equation)."  They  go  on  to  state  that  there 


Is,  however,  strong  circumstantial  evidence  for  the  occurrence  of 
internal  waves.  Some  of  that  evidence  will  now  be  presented. 
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In  1966  radar  tracking  of  radar-reflecti ve  balloons  was  used 
at  Cape  Canaveral  to  obtain  profiles  of  atmospheric  winds  from 
0 kilometers  to  about  16  kilometers  (Endiich,  et.  al.,  1969).  Two 
series  of  profiles  were  made,  each  series  lasting  between  one  and 
two  days.  The  method  allows  one  to  sample  rather  rapidly  in  the 
vertical  (one  point  each  25  meters)  but,  unfortunately,  it  can 
only  determine  the  wind  speed  at  each  point,  not  its  direction. 

The  authors  point  out  that  there  are  strong  velocity  features  in 
the  July  4-5,  1966,  experiment  which  progressively  move  downward 
with  time.  The  profiles  are  complicated,  however,  and  no  definite 
time  or  space  scales  stand  out  (although  it  is  known  that  internal 
waves  can  exist  In  the  atmosphere  at  these  altitudes). 

Although  apparent  high-frequency  features  in  the  atmosphere 
have  been  observed  to  change  altitude,  there  appear  to  have  been 
few  reports  of  any  polarization  or  directional  structure  associated 
with  these  features.  Endiich  could  only  measure  speed  as  a func- 
tion of  altitude  and  therefore  could  not  observe  polarization. 

Hines  (1966)  has  reported  polarization  of  atmospheric  winds  that  is 
apparently  related  to  the  diurnal  atmospheric  tide.  His  data,  ob- 
tained by  tracking  sodium  vapor  trails  at  heights  of  90-130  kilo- 
meters, Indicates  that  winds  with  periods  appropriate  to  the  diurnal 
tide  rotate  with  altitude  in  a given  rocket  profile.  If  an  observer 
is  imagined  to  be  standing  above  the  atmosphere  and  looking  down. 
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then  the  waves  seen  by  Hines  have  horizontal  velocity  vectors  that 
rotate  counterclockwise  with  increasing  "depth,"  or  atmospheric 
pressure.  This  is  the  opposite  of  the  case  observed  in  the  ocean 
by  the  FMVP  over  the  smooth  topoqraphy,  in  which  the  horizontal 
current  vector  rotates  clockwise  with  increasing  depth.  As  Hines 
points  out,  the  atmospheric  observations  indicate  that  the  wave 
energy  of  these  oscillations  is  propagating  from  lower  altitudes  up 
to  the  90-130  kilometers  range.  This  appears  to  be  one  of  the  few 
cases  where  wave  polarization  has  been  observed  in  the  atmosphere 
(C.  0.  Hines,  1975,  personal  communication). 

One  good  reason  for  a lack  of  elliptical  or  circular  polar- 
ization in  the  motions  interpreted  as  being  internal  waves  In  the 
atmosphere  is  the  fact  that  these  motions  quite  often  have  periods 
much  shorter  than  the  local  inertial  period.  Hines  (1960)  points 
out  that  atmospheric  internal  waves  might  be  expected  to  be  most 
energetic  at  periods  of  several  hours.  For  waves  of  this  period, 
the  earth's  rotation  is  not  a major  influence,  and  a single  atmos- 
pheric internal  wave  with  a three  hour  period  would  be  almost 
linearly  polarized  in  the  vertical. 

Hines  also  notes  that  some  observations  in  the  atmosphere 
show  structures  which  proaressi vel.y  move  downward  in  time.  Usinq 
the  same  arauments  as  those  put  forth  in  describing  Figure  21,  he 
points  out  that  because  the  directions  of  vertical  phase  propaga- 
tion and  energy  flux  are  opposite  for  internal  waves,  these  down- 
ward moving  structures  in  the  upper  atmosphere  could  be  interpreted 
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as  internal  waves  carrying  energy  upward,  away  from  the  lower  at- 
mosphere and  the  surface  of  the  earth.  (The  sources  for  atmospheric 
internal  waves  are  thought  by  many  meteorologists  to  be  located  in 
the  lower  atmosphere  or  at  the  earth's  surface.) 

Finally,  there  is  some  evidence  from  fresh  water  lakes  that 
appears  to  indicate  upward  phase  propagation  of  internal  waves. 
Lazier  (1973)  obtained  a series  of  temperature  profiles  in  Lake 
Bala,  Wales,  over  a time  interval  of  about  12  hours.  Profiles  were 
repeated  at  approximately  11  minute  intervals.  In  many  cases, 
features  in  the  temperature  profile  were  observed  to  move  progres- 
sively upward  in  time.  Lazier  attributes  this  upward  motion  of 
temperature  features  to  the  presence  of  propagating  internal  waves. 
If  the  energy  source  for  these  waves  was  located  at  the  lake  sur- 
face, they  would  have  to  propagate  their  energy  downward.  Then, 
as  in  Figure  21,  the  wave  number  vectors  of  the  waves  would  be 
directed  toward  the  lake  surface,  and  the  wave  phase  would  move 
upward  in  time.  Since  Lazier  only  measured  temperature  profiles, 
however,  no  data  on  wave  velocity  polarization  were  obtained. 
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Chapter  IV  The  Use  of  the  Ratio  of  Clockwise  to  Counterclockwise 
Energy  as  a Reflection  Coeffici ent.  _Re f lection  of  Internal  Waves 
From  a Smooth  Bottom 

a.)  Calculation  of  a reflection  coefficient  for  the  observed  waves 
We  pointed  out  in  Chapter  III  that  the  dominance  of  clock- 
wise over  counterclockwise  vertical  wave  number  spectra  for  hori- 
zontal kinetic  energy  appears  to  indicate  that  the  net  vertical 
energy  flux  for  these  waves  is  downward.  However,  a further  step 
must  be  taken  if  we  are  to  relate  the  partition  between  clockwise 
and  counterclockwise  energy  to  the  energy  in  upward  and  downward 
propaaating  waves.  This  step  is  to  assume  that  the  observed  waves 
have  frequencies  close  enough  to  f so  that  any  single  wave  should 
be  almost  circularly  polarized  in  the  vertical  (see  equation  1 1 1 - 1 3 ) . 
The  sense  of  polarization  is  determined  by  whether  the  wave  energy 
flux  is  downward  (clockwise  polarization)  or  upward  (counterclock- 
wise polarization).  With  this  assumption,  the  horizontal  kinetic 
energy  of  the  downward  propagating  waves  is  approximately  equal  to 
C(m'),  the  clockwise  energy  in  the  vertical  wave  number  spectrum, 
while  the  horizontal  kinetic  energy  of  the  upward  propagating  waves 
is  approximately  equal  to  A(m'),  the  counterclockwise  energy  in  the 
same  spectrum.  (Strictly  speaking,  A(m')  and  C(m')  should  be  multi- 


, plied  by  the  Brunt-Vaisala  frequency  to  get  the  true  energy  level 

(see  Chapter  V).  However,  since  we  are  finally  interested  only  in 

f 


the  ratio  of  clockwise  to  counterclockwise  energy,  this  step  has 
been  ianored.) 
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We  have  assumed  in  the  previous  paragraph  that  inertial 
waves  are  the  primary  cause  of  high-frequency  energy  in  the  pro- 
files. However,  frequency  spectra  of  horizontal  velocity  obtained 
at  the  central  mooring  (M.  Briscoe,  1975,  personal  communication) 
indicate  that  the  diurnal  and  semidiurnal  tides  also  contribute 
to  the  internal  wave  frequency  range  (see  Figure  29).  Because  of 
the  relatively  short  length  of  the  time  series,  it  is  difficult  to 
resolve  semidiurnal  tidal  and  inertial  oscillations  in  frequency, 
and  it  is  impossible  to  resolve  diurnal  tidal  and  inertial  oscilla- 
tions. On  the  other  hand,  we  can  make  inferences  concerning  the 
probable  influence  of  the  tides  on  the  vertical  wave  number  spectra 
computed  in  Chapter  II.  First,  that  part  of  the  tiual  energy  which 
is  caused  by  the  surface  or  barotropic  tide,  since  it  is  depth- 
independent,  is  not  sensed  by  the  EMVP.  Second,  work  of  other  in- 
vestigators (Hendry,  1974)  appears  to  show  that  most  of  the  tidal 
energy,  at  least  for  the  semidiurnal  tide,  is  concentrated  in  the 
first  two  or  three  barocllnic  modes  (along  with  the  surface  tide). 
In  terms  of  the  vertical  wave  number  spectra  of  Chapter  II,  this 
means  that  any  tidal  contributions  to  the  spectra  are  located  in 
the  lowest  and,  possibly,  the  second-lowest  wave  number  estimates. 
We  may  then  infer  that  above  the  first  one  or  two  vertical  wave 
number  estimates,  the  primary  contributions  to  the  spectra  come 
from  inertial -internal  waves.  There  is  also  direct  evidence  from 
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the  profiles  to  indicate  that  most  of  the  semidiurnal  tidal  energy 
is  contributed  by  waves  with  small  vertical  wave  numbers.  This 
evidence  will  be  presented  in  Chapter  V. 

We  now  assume  that  the  observed  waves  are  generated  at  or 
near  the  surface  and  propagate  their  energy  downward.  Upon  reach- 
ing the  bottom,  some  fraction  of  their  energy  is  reflected  upward 
toward  the  surface,  while  the  remainder  of  the  energy  is  lost  in  a 
frictional  boundary  layer  at  the  bottom.  There  may  also  be  some 
energy  loss  as  the  waves  propagate  through  the  interior.  We  also 
assume  that  the  smooth  bottom  does  not  act  as  an  energy  source  for 
these  waves  (that  is,  a boundary  layer  at  the  bottom  will  remove 
energy  from  the  incoming  waves,  but  other  processes  do  not  act  to 
put  energy  into  waves  that  are  reflected  from  the  bottom).  The 
theoretical  model  presented  later  in  this  chapter  will  hopefully 
indicate  how  reasonable  the  above  assumptions  are. 

If  we  accept  these  assumptions,  then  the  ratio  A(m * )/C (m ' ) 
can  be  thought  of  as  a reflection  coefficient  for  these  waves. 

This  ratio  is  shown  in  Figure  22,  using  the  spectra  of  Figures  13 
and  16.  Since  the  bottom,  by  assumption,  cannot  act  to  put  energy 
into  the  waves,  a good  test  of  the  above  model  is  that  A(m')/C(m') 
should  always  be  less  than  one.  An  examination  of  Fioure  22  shows 
that  this  is  substantially  true,  although  at  high  wave  numbers  (or 
low  energy),  the  ratio  A(m')/C(m')  behaves  erratically.  This  er- 
ratic behavior  at  low  energy  levels  may  be  an  indication  that  the 
signal-to-noise  ratio  of  the  FMVP  is  becoming  small.  A low 
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polarized  results  in  an  elliptically  polarized  wave).  This  argument 
indicates  that  A(m')/C(m')  is  probably  an  over-estimate  of  the  re- 
flection coefficient  for  the  near-inertial  waves.  The  reader  may 
refer  to  Appendix  D for  a more  complete  discussion  of  the  influence 
of  hiqher  frequency  waves  on  the  clockwise  and  counterclockwise 
spectra. 

b.)  A model  for  the  reflection  of  internal  waves  from  a smooth 
bottom 

A theoretical  model  for  internal  wave  reflection  can  be 
developed  that  supports  some  of  the  assumptions  made  in  the  pre- 
vious section.  This  model  is  a variation  of  one  studied  by  Phil- 
lips (1963)  that  analyzes  the  reflection  of  purely  inertial  waves 
from  a flat  boundary.  The  present  model  includes  buoyancy  effects 
but  otherwise  gives  much  the  same  results  as  Phillips'  model. 

We  start  with  the  usual  equations  for  internal  waves,  with 

buoyancy  and  rotation: 

XL  - fv  = - 3E.  + vV2u  ■ 

at  ax  > 
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The  reader  can  refer  to  Appendix  B for  definitions  of  the 
quantities  in  the  above  equations.  The  only  additional  terms  in 
the  above  equations  are  the  frictional  terms  vV2u,  vV2v  and  vV2w, 
where  v is  the  kinematic  viscosity  coefficient.  The  Brunt-Vaisala 
frequency,  N,  is  assumed  to  be  constant.  This  is  not  a serious 
assumption,  since  near  the  bottom  where  we  expect  the  model  to  ap- 
ply, N is  fairly  constant  with  depth. 

If  all  variables  are  independent  of  the  y coordinate,  a 
stream  function,  ip,  can  be  introduced  that  satisfies  continuity: 

u = -^  and  w = - M 

3Z  . IV-2 


If  \p  is  substituted  into  equations  IV-1,  these  equations  can  be 
reduced  to  two  equations  in  ip  and  v: 

[_L  v2  (_3  . vVz)  + H2  -it  ] ,i,  - f &L-  = 0; 

Lit-  'at-  ' »v2  J ' i + i-r 


f M + [-2.  - vV2l  v = 0 
3z  L3t 


IV-3 


The  coordinate  origin  is  assumed  to  be  located  at  the  bot- 
tom, so  that  z = 0 there,  and  the  coordinates  form  a right  hand 

system.  Then  equations  IV-3  must  satisfy  the  boundary  conditions 


We  look  for  solutions  of  the  form: 


kyz  i(kx-wt)- 
tp  = A e e > 


kyz  i(kx-uat) 
v = V e e 
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If  these  expressions  for  ^ and  v are  substituted  in  equations  IV-3 
there  results  a sixth  order  equation  in  y: 


[R2]  y6  + [2iw'R  - 3R2]  y" 


+ [1  - oj’2  - 4io)'R  + 3R2  + i — R]  y2 

fw 


+ [U1 


2 - N 


+ 2ioj'R  - R2  - i i_  R]  = 0 , 
f2  fu) 
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where  <*)'=  — and  R = . If  we  assume  that  R = 0,  the  charac- 

f f 

teristic  equation  for  y becomes 

N2  - w2 


y 


2 - _ 


W2  - f2 


. /N2  - ,,2v  1/2 

or  y = + i(-2 

1 «2  w2  - f2 


I V- 7 


These  two  roots  give  the  z-dependence  that  would  be  expected  from 
two  internal  waves  propagating  without  friction. 

If  we  now  assume  R « 1 , we  can  factor  out  the  above  two 
roots  to  get  approximate  expressions  for  the  remaining  four: 
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, V 
,2  = - + 


(l/2)[-  _1-  li  (Ni-  V - 2 

P2  R uf 


Ni-Ud)]’72 

u>2  - f2 


If  we  assume  that  the  - JL  term  dominates  under  the  radical,  then 

R2 

the  expression  for  the  four  remaining  roots  becomes 

u2  = -1  (V  + 1)  . IV-9 


Noting  that  u is  assumed  positive,  we  find  that  the  two  boundary 
layer  roots  which  give  solutions  that  decay  into  the  interior  are 


* J-'-l  (a)'  * 1)1/2 

».*•  / 2R 
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The  four  roots  of  jj  that  have  been  obtained  for  the  two  propagating 
internal  waves  and  the  two  boundary  layer  solutions  are  then: 

u - i (^-^)1/2; 

1 a)2  - f2 

_ , /N2  - w2\l/2. 

VJ  - - i v ) * 

2 O)2  - f2 

V = LjlI  (w‘-l)1/2; 

3 S2R 


= i_l_i  (w'+l) 


IV-11 
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Ue  note  that  the  thickness  of  the  boundary  layer  is  determined  by 

1 '2 

the  quantity  (to1  - 1)  ' //^R,  for  w1  positive  and  close  to  one.  As 
u) 1 -*-1,  or  as  w ■»  f,  the  boundary  layer  penetrates  more  and  more 
deeply  into  the  interior. 

The  expressions  for  v and  ip  are  now 


ip  = [A^e 
v = [Vte 


kv^z 


kuiz 


A e 
2 


V2e 


ku2z 


kp2z 


A e 

3 


Vje 


ku3z 


ku3z 


ky,,z,  i(kx  - wt) 
+ A e ] e 

4 


kp,z  i (kx  - wt) 
+ V e * ] e 

4 J 
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where  vij  through  p4  are  given  above.  The  A’s  and  V‘s  must  be  deter- 
mined from  the  boundary  conditions  at  the  bottom  (equations  IV-4). 
The  terms  with  V3  and  in  the  expressions  for  \p  and  v give  the 
amplitude  of  the  incoming  wave,  while  the  terms  involving  V2  and  A2 
oive  the  amplitude  of  the  reflected  wave.  The  average  horizontal 
kinetic  energy  for  the  incoming  wave  is  given  by 

E.  = IV  |2  (1  + ^)  . IV-13 

11  f2 

The  average  horizontal  kinetic  energy  for  the  reflected  wave  is 


E 


o 


|vj2  (i 


I V - 1 4 


where 


A = u>'a 


— i-£fi _ i-^R(w'-l)  f 1 1 

(i-l)L'-l)1/2  2(1-1)  Vtco’-l  )V2  (oo'+l  )1/2 

IV-15 
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and 


ui  a 


i A* 

(i-l)U'-l) 


1/2 


+ i ^R(u/-U 

2(i-D 


(w'-l)1/2  (ui'+l)172 


IV-16 


and 


a 


(m2  - 
N2  - a)2 


IV-17 


Using  and  EQ,  we  can  form  a reflection  coefficient  for  the  hori- 
zontal kinetic  energy: 

h.  - , - M <Ni  ~ “,)l/2  f + 0(R). 

E u>(u,l2-l)1/2  ( (W-i)1/2  (u>'+l)1/2j 

1 I V - 1 8 

This  reflection  coefficient  has  several  characteristics  that 
are  relevant  to  the  present  data.  First,  although  for  dimensionless 
frequencies  u>’  much  greater  than  one,  EQ/E^  ~ 1 , as  u>'  -*•  1 this  re- 
flection coefficient  decreases  rapidly  (EQ/Ei  does  not  limit  to  zero 
as  iii  + f due  to  the  approximations  made  in  the  above  derivation). 
This  shows  that  energy  in  waves  near  the  inertial  frequency  may  be 
greatly  reduced  as  these  waves  reflect  off  the  bottom.  Most  of  the 
energy  loss  suffered  by  these  waves  takes  place  in  the  bottom  fric- 
tional boundary  layer  (a  boundary  layer,  we  recall,  that  grows  in 
thickness  as  u>  ■*  f).  For  frequencies  close  to  f,  EQ/Ei  is  very 


sensitive  to  w.  Since  we  cannot  resolve  ui  as  closely  as  would  be 
required  to  test  the  above  theory,  we  can  only  say  that  the  ob- 
served "reflection  coefficients"  (Figure  22)  and  the  possible 
boundary  layer  thickness  (described  below)  can  be  accounted  for  if 
the  frequencies  of  the  observed  waves  are  close  enough  to  f. 

Figure  8 appears  to  indicate  that,  from  2000  dbar  to  about 
5300  dbar,  at  least,  the  average  perturbation  horizontal  kinetic 
energy  follows  closely  the  average  profile  of  N.  This  would  seem 
to  imply  that  the  waves  are  not  undergoing  strong  frictional  losses 
as  they  propagate  through  the  water  column,  and  therefore  that  most 
of  the  energy  loss  from  these  waves  takes  place  as  they  reflect  off 
the  bottom. 

There  is  also  a suggestion  in  Figure  8 that  the  mean  per- 
turbation kinetic  energy  decreases  within  about  100-150  dbar  above 
the  bottom.  This  may  indicate  the  presence  of  a frictional  boundary 
layer  at  the  bottom.  In  theory,  the  horizontal  kinetic  energy  must 
go  to  zero  at  the  bottom.  However,  the  EMVP  only  measures  horizon- 
tal velocity  to  about  10  dbar  above  the  bottom.  Thus,  it  probably 
would  not  see  that  part  of  the  frictional  boundary  layer  nearest 
the  bottom,  where  the  greatest  energy  loss  takes  place.  The  error 
estimate  (Chapter  II)  calculated  for  the  mean  perturbation  kinetic 
energy  profile  of  Figure  8 indicates  that  the  decrease  in  energy 
near  the  bottom  is  significant.  We  again  must  emphasize,  however, 
that  this  error  estimate  only  takes  into  account  the  influence  of 


instrumental  noise. 


fioure  8 also  shows  an  apparent  bulge  in  the  mean  pertur- 
bation horizontal  kinetic  energy  at  about  200  dbar  above  the  bot- 


88 


l 


L 


tom.  This  bulge  remains  even  if  the  kinetic  energy  profile  is  nor- 
malized by  the  mean  Brunt-Vaisala  profile  (N  also  increases  as  the 
bottom  is  approached).  The  theory  presented  above  predicts  only 
that  the  mean  perturbation  horizontal  kinetic  energy  should  increase 
smoothly  from  zero  at  the  bottom  to  values  appropriate  to  the  in- 
coming and  outgoing  waves  beyond  the  frictional  boundary  layer. 
Therefore,  this  apparent  bulge  in  the  kinetic  energy  near  the  bot- 
tom is  not  explained  by  the  above  theory.  The  possible  existence 
of  small  bottom  slopes  may  be  crucial  here.  The  above  theory  as- 
sumes a flat,  horizontal  bottom.  Wunsch  (19G9)  has  shown  that  if 
the  group  velocity  vector  for  the  reflected  wave  has  an  inclination 
angle  to  the  horizontal  which  is  close  to  the  bottom  slope,  an  in- 
tensification of  the  internal  wave  energy  near  the  bottom  can  take 
place.  Since  the  group  velocity  vectors  for  waves  near  the  iner- 
tial freguency  are  only  slightly  inclined  to  the  horizontal,  the 
near-bottom  behavior  of  these  waves  may  depend  stronnly  or.  small 
bottom  slopes. 
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Chapter  V Comparison  of  Vertical  Wave  Number  Spectra  with  Theoret- 
ical Models.  Computation  of  Dropped,  Time-lagged,  Rotary  Coherence 
(DLRC).  The  Vertical  Energy  Flux  of  Internal-inertial  Waves. 

a.)  Vertical  wave  number  spectra 

Recently,  several  models  have  been  derived  which  describe 
the  amount  of  energy  contributed  to  a region  of  the  ocean  by  in- 
ternal waves  of  various  length  and  time  scales.  Perhaps  the  pri- 
mary model  of  this  type  Is  that  of  Garrett  and  Munk  (1972,  1975). 
These  authors  have  constructed  an  internal  wave  energy  density 
spectrum,  E(a,  to),  as  a function  of  horizontal  wave  number,  a,  and 
frequency,  to.  The  energy  spectrum  is  assumed  to  be  isotropic; 
that  is,  the  magnitude  of  E is  independent  of  the  direction  of  a. 
The  general  form  of  E(a,  to)  is  (Garrett  and  Munk,  1972,  eq.  6.3): 

E(ot,to)  = Cp'1A(A)f2(to).  V-l 

This  is  a dimensionless  equation.  Time  has  been  made  dimensionless 
by  a reference  Brunt-Vaisala  frequency, 

N = 3 cph,  V-2 

o 

and  distance  by  a reference  cyclical  wave  number, 

M = 0.122  cpkm.  V- 3 

o 


C is  a constant,  w is  called  the  "a-bandwidth  of  the  spectrum"  as 
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a function  of  frequency.  The  parameter  A is  qiven  by 

X = a/u  . V-4 


The  form  chosen  for  fi(w)  was 
C>t  \ "P+^S  / 2 

Q((jd)  = 0)  (u>  - UK  ) , 


V-5 


where  u^-  is  the  (non-dimensional)  local  inertial  frequency.  The 
constants  p and  s were  chosen  to  be  2 and  1/2,  respectively.  Then 


J2(u>) 


1 


U)(w2 


1/2  . 


V- 6 


In  the  1972  version  of  the  model,  A(A)  was  chosen  to  have  the  form 
of  a "top  hat": 

A(A)  = 1 , for  0 £ X <_  1 

and  A(A)  = 0,  for  X > 1.  V-7 

Finally,  u(go)  was  chosen  to  be 

u(oj)  = j^r(u)2  - w.2)1/2,  V-8 

where  is  an  integer.  In  the  1972  model,  j was  called  the 
"equivalent  number  of  modes  at  the  inertial  frequency"  by  the 
authors.  (It  should  be  emphasized  that  their  model  does  not  really 
depend  on  the  existence  of  modes  in  the  vertical.  The  parameter  j* 
can  be  thought  of  as  characterizing  the  maximum  vertical  wave  number 
found  to  contain  significant  energy  in  the  internal  wave  energy 
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density  spectrum).  We  would  now  like  to  see  how  this  theoretical 
model  compares  to  vertical  wave  number  spectra  observed  by  the  EMVP. 

In  order  to  proceed  further  we  must  introduce  the  dimension- 
al forms  of  the  above  equations.  In  what  follows,  a hat  (~)  above 
a quantity  will  indicate  that  it  is  dimensional.  A prime  (')  will 
indicate  that  vertical  wave  numbers  in  the  stretched  coordinate 
are  being  used. 

The  dispersion  relation  for  internal  waves  (under  a WKB 
approximation)  gives 

a = m(z)  (S.2  for  ui  « f)(z) 

A(z)  v-9 


where  m is  the  local  vertical  wave  number.  Also, 

. s , , , J/2  _ i (w2  - P)V2 

p " J*TT  I U)  - - J*TT  -1 

0 


V-10 


Then 


X = a/y  = 


& "(2)  “o 


UM0  (j*u)  M0f)(z) 


V-ll 


However,  since  in  the  stretched  vertical  coordinate 


m'  = m(z) (NQ/N(z) ) , 


we  have 


V-12 


(j**)  Mn  . 

o 


V - 1 3 


Then  the  dimensional  form  of  fi(w)  becomes 
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n(u)  = 


(£2  - fz)1/2  ’ 


/V  / /V 

CO 


and 


E(i,m')  = C A(x) 


(j*ir)  £(£*  - f2)  . 


V - 1 4 


V-15 


The  total  energy  in  horizontal  velocity  is  given  by 


Et0Ul  ' f S'  E(“’“) 


IF  dado) 


a o) 


where 


U»  = Nlzl 
No 


V-16 


(Garrett  and  Munk,  1972).  V-17 

However,  since  velocity  in  the  profiles  has  been  normalized  by 
N(z)/Nq,  we  must  use 


■ 2 , _ 


tt 


(A(z)/no) 


V -1 8 


instead  of  U in  equation  V-16.  Converting  a and  <d  to  their 
dimensional  counterparts , where 


a = ot/ M_  and  u>  = w/N  , 
o o 


V - 1 9 


we  have 


’total  M 


— ) \ U2 ' E(a,o))  da  du> 


V-20 


oo 


a w 


I 

L ... 
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Instead  of  Integrating  over  S,  we  wish  to  carry  out  the  integration 


over  m'.  Since 

da  * (ui2  - ?2-)-1/2  dA  V-21 

fi(z) 

and 

dm'  = (N0/fl(z))  dm,  • V-22 


the  integral  becomes 
E 


total 


JL  J C IP  ' E(ffi ' ,£j)  tu2  _i2) 

N0 


1/2 


dm'  dw 


oo  m a) 


Mn  1 / a) 


£*+i^  (u2  - f2) 


1/2 


/v 

dm'  du 


o o 


m'uj 


■ H 

(j*ir,M0  m'i 


(a2  + f2J  - 
£3(£2  - ?2)1/2 


A ( 


m * 

m 


j^M0 


/V  A 

dio  dm 1 . 
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Etotai  must  **e  multiplied  by  pN2M0‘3  to  obtain  the  dimensional 
energy  spectrum,  where  p is  sea  water  density  (assumed  equal  to 
1 gm/cm3).  Removing  the  integral  over  m'  leads  to  the  following 
dimensional  form  of  the  energy  density  spectrum  in  stretched  verti- 
cal wave  number,  m': 


Ea.(fi')  * 


m 


ecN0 

j*7TMr 


A [-^r) 

J*nM0 


X 


U) 


+ f2 


dco 


u> 


'(w2  - f2) 


1/2 


3Ep_ 

2j*TT 


A ( 


_JL_) 

J>Mo 


V-24 


( f.  Garrett,  1974,  personal  communication). 
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E„  is  a dimensionless  constant  which  is  related  to  C by 
o 

Efl  s C ^ fi(uO  dw  . 

“1 


V-25 


A A 

To  obtain  the  actual  energy  levels  at  a given  depth,  zQ,  E~,(m') 
is  multiplied  by  (N(z)/NQ).  At  that  level 

m'  = m(z0)  (N0/N(zq)).  V- 26 

A 

The  quantity  E~,(m'),  as  obtained  by  observations,  has  been 

given  already  in  Figures  12  and  15.  Clearly,  the  form  of  the 

stretched  vertical  wave  number  spectrum  depends  on  A(x).  In  the 

1972  model,  A(X)  was  chosen  to  be  a top  hat.  Comparing  equations 

V-7  and  V-13,  we  see  that  A(X)  * 1 for  m'  <_  and  A(X)  = 0 for 

in'  > j itM  . (Again  we  see  that  in  the  top  hat  model  j.  controls 
w o * 

the  maximum  vertical  wave  number  that  contributes  energy  to  the 

A 

spectrum).  This  means  that  the  theoretical  E~,  is  independent  of 
m'  for  X < 1,  and  is  zero  for  X > 1.  It  is  clear  from  Figures  12 
or  15  that  £jjj,  is  not  Independent  of  m',  but  in  fact  decreases  with 
increasing  in' . 

In  a revised  form  of  their  model  (Garrett  and  Hunk,  1975), 
these  investigators  chose  a new  form  for  the  function  A(X): 

A(x)  - (t  - 1)  (1  + X)  _t  V-27 

where  t is  chosen  to  fit  observations.  Garrett  and  Hunk  picked 

t = 2.5.  V-28 
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(The  choice  of  t * 2.5  appears  to  agree  reasonably  well  with  verti- 
cal displacement  spectra  obtained  by  Hayes,  et.  al .(1975)).  For 
NQ  = 3 cph,  MQ  * 0.122  cpkm,  E0  = 2it  x 10”^  and  j*  = 6 (in  the 
1972  model,  j*  was  20),  the  vertical  stretched  wave  number  spectrum 
becomes 


and 


E~((m*)  = 19.1  A(X)  Lem/ seel2 
m cpkm 

X = 0.435  m*  . 


V-29 

V-30 


We  have  repeated  the  stretched  vertical  wave  number  spectrum 
of  Figure  12  in  Figure  23.  Figure  23  also  shows  the  theoretical 

a A 

stretched  vertical  wave  number  spectrum,  Ejjj.(m'),  using  the  form  of 
A(X)  given  by  equation  V-27,  for  three  choices  of  j*:  j*  = 6, 

J*  = 15  and  = 20.  (The  cyclical  frequency  has  been  converted 
from  cycles  per  kilometer  to  cycles  per  sdb).  The  value  of  the 
Constant  for  the  spectral  energy  level,  E(j*),  (corresponding  to 
tQ  in  the  above  derivation)  has  been  adjusted  to  make  the  three 
theoretical  curves  have  the  same  energy  levels  at  large  stretched 
Vertical  wave  numbers.  If  E(j*)  is  the  actual  value  of  the  con- 
stant used  for  a given  j*,  and  EQ  = 2tt  x 10"^,  we  have  the  follow- 
ing values  for  E(jJ:  E(j*  * 6)  = 2 x Efl;  E(j*  * 15)  = 0.6  x E0; 
E(j*  ■ 20)  * 0.4  x Eq.  The  theoretical  curve  for  a of  15 
appears  to  best  represent  the  observed  stretched  vertical  wave  num- 
ber spectrum,  both  in  the  general  shape  of  the  spectrum  and  in  the 
total  energy  level. 


VERTICAL  WAVE  NUMBER 
(C/SOB) 

Figure  23.  Comparison  of  theoretical  and  observed  total  energy  spectra 
in  stretched  vertical  wave  number  (see  text). 
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A value  of  15  for  j*  is  larger  than  the  value  of  6 assumed 
in  the  1975  model.  If  we  return  for  the  moment  to  a consideration 
of  Figures  10  and  11,  the  reason  why  j*  = 6 Is  apparently  too  small 
becomes  clearer.  Visually,  Figures  10  and  11  show  that  waves  con- 
taining relatively  high  energy  have  wave  lengths  as  short  as 
100  sdb.  The  total  depth  in  the  stretched  vertical  coordinate  is 
1860  sdb.  If  we  now  assume  that  j*  is  in  some  sense  representative 
of  the  greatest  number  of  stretched  vertical  wave  lengths  that 
"fit"  into  1860  sdb  (and  contain  significant  energy),  we  find  that, 
for  199-130  sdb  vertical  wave  lengths,  should  be  at  least  15. 

It  is  also  interesting  to  note  that  the  value  of  E required 
to  make  the  theoretical  (j*  * 15)  and  observed  stretched  vertical 
wave  number  spectra  contain  (roughly)  the  same  total  energy  is 
about  60%  of  the  constant,  Eq,  chosen  by  Garrett  and  Munk.  The 
fact  that  high-frequency  Internal  waves  (those  with  frequencies 
near  N)  are  probably  not  energetic  enough  to  be  sensed  by  the  EMVP 
may  partially  explain  this.  However,  it  is  also  likely  that  the 
total  energy  of  the  Internal  wave  field  at  the  MODE  central  moor- 
ing Is  less  than  the  total  energy  at  other  locations  (at  Site  0 
(39°  20'  N,  70°  00'  W),  for  example).  Since  observations  at  Site  D 
were  used  (at  least  In  part)  to  determine  E0,  the  fact  that 
E(j*  s 15)  Is  less  than  E^  may  also  reflect  an  actual  decrease  in 
Internal  wave  energy  at  the  central  mooring,  relative  to  other 
locations  in  the  ocean  where  measurements  have  been  made. 
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b.)  Dropped,  time-lagged,  rotary  coherence  of  profiles 

The  following  calculations  were  carried  out  In  order  to 
compare  profiler  horizontal  velocity  data  with  some  theoretical  re- 
sults based  on  the  Garrett  and  Munk  model  presented  In  section  a. 

(C.  Garrett,  1974,  personal  communication) . Since  we  have  a time 
series  of  profiles.  It  is  possible  to  calculate  the  velocity  coher- 
ence between  two  profiles  separated  by  a time  Interval,  T,  for  any 
stretched  vertical  wave  number.  This  Is  called  the  dropped,  lagged, 
coherence  (or  DLC)  by  Garrett  and  Munk.  However,  we  know  that  much 
of  the  variability  In  the  profiles  Is  associated  with  waves  having 
periods  near  the  local  Inertial  period.  For  this  reason.  Instead 
of  the  DLC  we  have  calculated  the  coherence  between  two  profiles 
when  the  horizontal  velocities  of  the  second  profile  have  been 
rotated  through  an  angle  y: 

y * uT, 

U)  = 2*f,  V-31 

where  f is  the  Inertial  frequency  and  T is  the  time  interval  be- 
tween the  two  profiles.  This  Is  called  the  dropped,  lagged,  rotary 
coherence  (or  DLRC).  If  the  first  velocity  profile  Is 

W,(z*)  * U,(z*)  + 1Vj(z*);  V- 32 

and  the  second  profile  Is  given  by 

W2(z*)  = U2(z*)  + 1V2(z*),  V- 33 

then  the  coherence  Is  calculated  between  W,  and 

W,  * W2eiY  « U,(z*)  + IV, (D  V-34 
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We  then  obtain  the  following  four  quantities: 


PU1U3  = 


PU1U3 


[PU1U1]1/2[PU3U3]1/2 


pviv3  = Em3 

[PV1V1]1/2[PV3V3]1/2 


QU1U3 


[PU1U1 ]1/2[PU3U3]1/2 


QV1V3  = 


[PV1V1]1/2[PV3V3]1/2  * 


where  PU1U3  is  the  normalized  cospectrum  between  velocity  compo- 

A 

nents  U,  and  U3,  QU1U3  is  the  normalized  quadrature  spectrum  between 
the  same  components,  and  PV1V3  and  QV1V3  are  the  same  quantities 
for  the  velocity  components  V,  and  V3.  PU1U1,  PU3U3,  PV1V1  and 
PV3V3  are  the  autospectra  for  Uj,  U3,  V,  and  V3,  respectively. 

All  of  the  above  quantities  are  functions  of  stretched  vertical 
wave  number  and  the  time  lag  between  profiles. 

We  now  describe  the  actual  method  of  computation.  Time  lags  , 
T,  were  calculated  for  all  pairs  of  drops  in  the  time  series.  These 
pairs  were  then  sorted  in  such  a manner  that  all  pairs  with  lags 
within  + 1/2  hour  of  a given  whole  hour  were  assigned  a lag  equal 
to  that  whole  hour.  For  example,  all  pairs  with  time  lags  of  from 
2 1/2  to  3 1/2  hours  were  given  a lag  equal  to  3 hours  for  the  cal- 
culation. Calculations  were  done  only  for  those  lags  which  had 
four  or  more  pairs  of  drops.  The  results  of  these  calculations  are 
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shown  in  Figures  24,  25,  26  and  27  for  lags  up  to  almost  two  days. 
For  any  lag  at  which  the  calculation  was  made,  there  are  three 
points  shown  for  each  of  the  four  normalized  spectral  quantities 
given  above.  Each  of  these  points  represents  an  average  over  all 
pairs  available  at  that  time  lag.  In  addition,  estimates  have 
been  averaged  over  several  stretched  vertical  wave  number  bands. 

The  number  n represents  the  number  of  the  stretched  vertical  wave 
number  estimate;  that  is,  the  vertical  wave  length  is  given  by 
1800/n  sdb,  where  1800  sdb  is  the  total  length  of  the  series. 

Thus,  the  three  points  (at  a given  lag  and  for  a given  quantity) 
are  averages  over  stretched  vertical  wave  number  estimates  n = 1-5, 
n * 6-10,  and  n 3 11-15. 

It  is  somewhat  difficult  to  put  error  bars  on  PU1U3,  PV1V3, 

/N  A 

QU1U3  and  QV1V3  since  each  point  may  be  the  result  of  averaging 
over  a different  number  of  pairs  of  drops.  In  an  attempt  to  esti- 
mate the  error  we  have  calculated  Puil)3,  0V1V3,  QU1U3  and  QV1V3 
separately  for  each  pair  of  profiles  in  a case  when  four  pairs  of 
profiles  were  used  to  give  the  estimates  shown  in  Figures  24-27. 

(We  recall  that  the  smallest  number  of  pairs  used  was  four.)  The 
resulting  values  of  PU1U3,  PV1V3,  QU1U3  and  QV1V3  have  an  average 
standard  deviation  about  their  mean  values  of  about  +0-2.  Pre- 
sumably, points  which  are  obtained  by  averaging  over  more  than  four 
pairs  of  drops  are  even  more  stable;  that  is,  the  scatter  of  points 
about  their  mean  values  is  less  than  t 0-2  standard  deviation. 


T (HOURS) 


Figure  24.  Dropped,  lagged,  rotary  cospectrum  (normalized)  between  velocity  components 
U,  and  Uj  (see  text  for  a further  description  of  Figures  24  to  27). 
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Figure  25.  Dropped,  lagged,  rotary  cospectrum  (normalized)  between  velocity 
components  V,  and  V3. 


T (HOURS) 


T (HOURS) 


Figure  27.  Dropped,  lagged,  rotary  quadrature  spectrum  (normalized)  between 
velocity  components  and  V3. 
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Figures  24  and  25  also  show  the  curve  of  PU1U3  or  PV1V3  ob- 
tained theoretically  (C.  Garrett,  1974,  personal  communication)  from 
the  Garrett-Munk  model  under  the  assumption  of  an  isotropic  spec- 
trum. The  curve  for  PU1U3  or  PV1V3  does  not  depend  on  wave  propa- 
gation in  the  vertical.  In  contrast,  Qtll U3  and  §V1V3  do  depend  on 
vertical  propagation.  If  only  vertical  modes  exist,  QU1U3  and  QV1V3 
are  theoretically  always  equal  to  zero.  On  the  other  hand,  if  all 
waves  are  propagating  energy  downward,  the  theoretical  curve  for 
QU1U3  or  QV1V3  starts  out  at  zero  for  zero  lag  and  slowly  increases 
to  about  + 0.25  for  a lag  of  50  hours.  If  we  accept  the  above 
"confidence  limit"  of  + 0.2,  it  is  difficult  to  say  whether  the  ob- 
served values  of  QU1U3  and  QV1V3  are  significantly  different  from 
zero.  Therefore,  the  following  paragraphs  will  be  limited  to  a dis- 
cussion of  the  observed  behavior  of  PU1U3  and  PV1V3.  These  quan- 
tities appear  to  contain  the  most  interesting  information  on  the 
behavior  of  the  DLRC  with  varying  time  lag  and  stretched  vertical 
wave  number. 

Using  the  above  error  estimate,  it  is  possible  to  say  sev- 
eral things  about  PU1U3  and  PV1V3.  First,  the  behavior  of  PU1U3 
or  PVIV3  for  estimates  containing  the  largest  vertical  wave 
lengths  (n  = 1-5)  is  significantly  different  from  the  behavior  of 
these  quantities  for  shorter  vertical  wave  lengths  (n  = 6-10  and 
11-15).  The  estimates  for  n = 1-5  decrease  much  more  rapidly  with 
increasing  lag  than  the  estimates  for  higher  n.  In  fact,  the  ob- 
served PU1U3  and  PV1V3  curves  for  n = 1-5  have  an  oscillatory 
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character  with  minima  (greatest  negative  values)  for  lags  of  about 
12  hours,  and  again  for  lags  of  about  30-40  hours.  This  Is  a 
strong  indication  that  the  longest  vertical  wave  lengths  have  sig- 
nificant energy  contributions  from  waves  with  periods  near  12  hours. 
The  most  obvious  candidate  for  this  energy  contribution  is,  of 
course,  the  barotlinic  semidiurnal  tide.  In  contrast,  the  esti- 
mates for  higher  n show  little  evidence  of  energy  contributions  at 
approximately  12  hour  periods.  This  implies  that  most  of  the 
semidiurnal  tidal  energy  is  in  the  longest  vertical  wave  lengths. 
Other  work  (Hendry,  1974)  has  shown  that  the  baroclinic,  semi- 
diurnal tide  has  most  of  its  energy  in  the  lowest  three  vertical 
modes.  The  above  result  appears  to  confirm  this. 

Second,  PU1U3  and  PV1V3  for  higher  n (n  = 6-10,  11-15) 
follow  the  theoretical  curve  more  closely.  Given  the  above  "error" 
estimate,  PU1U3  and  PV1V3  curves  for  n = 6-10  and  n = 11-15  cannot, 
in  general,  be  distinguished.  The  theoretical  carve  does  not  de- 
pend on  stretched  vertical  wave  number,  but  the  error  is  probably 
large  enough  so  that  dependence  on  vertical  wave  number  cannot  be 
tested  adequately.  It  is  interesting  to  note  that  PU1U3  and  PV1V3 
for  n = 6-10,  11-15  are  almost  always  greater  than  the  theoretical 
curve.  A possible  explanation  for  this  is  that  while  the  theoret- 
ical curve  takes  into  account  all  internal  wave  frequencies  (from 
inertial  to  Brunt-V'disSlS) , the  sensitivity  of  the  profiles  is  such 
that  higher  frequency  internal  waves  are  "filtered  out."  These 
higher  frequencies  tend  to  make  £lllU3  or  Pvi V3  decrease  more 


rapidly  with  increasing  laq  than  would  be  true  if  only  near-inertial 
waves  were  beinq  observed.  The  results  of  this  section  may  be  sum- 
marized as  follows:  1.)  The  behavior  of  PU1IJ3  and  PV1V3  for  wave 

lengths  outside  of  the  reoion  in  which  semidiurnal  tidal  energy  is 
sinnificant  does  not  appear  to  refute  the  theoretical  calculations, 
taking  into  account  the  possible  sources  of  error  described  above. 
2.)  The  result  of  other  calculations,  showinn  that  most  of  the 
semidiurnal  tidal  energy  is  confined  to  the  largest  vertical  wave 
lengths,  is  also  obtained  from  the  present  calculations. 

c.)  The  vertical  energy  flux  of  near-inertial  waves 

We  have  pointed  out  that,  because  of  the  relatively  short 
time  interval  over  which  the  time  series  was  obtained,  it  is  dif- 
ficult to  obtain  good  frequency  resolution  for  the  observed  waves. 
However,  there  are  several  independent  characteristics  of  the  data 
which  lead  us  to  believe  that  the  waves  which  we  have  been  calling 
"inertial"  actually  have  periods  somewhat  shorter  than  the  inertial 
period. 

First,  a frequency  spectrum  of  horizontal  velocity  obtained 
near  the  central  mooring  (Figure  29)  appears  to  show  that  what 
we  mioht  call  the  "inertial  peak"  is  actually  displaced  toward 
frequencies  somewhat  greater  than  the  local  inertial  frequency. 
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Second,  consideration  of  Figure  18  shows  that  many  of  the 
more  energetic  velocity  features  which  can  be  Identified  between 
drops  235U  and  236U  (these  drops  being  separated  by  4.8  kilometers) 
do  not  occur  at  the  same  depth  In  both  profiles,  but  rather  are 
offset  slightly.  This  depth  offset  is  too  large  to  be  accounted 
for  by  errors  in  the  EMVP  pressure  calibration.  The  Implication  Is 
that  the  phase  planes  of  these  waves  are  tilted  slightly  away  from 
horizontal.  Figure  21  then  Indicates  that  the  wave  number  vector, 
k,  for  these  waves  must  be  tilted  at  a small  angle  to  the  vertical, 
and  equation  III-10  shows  that  the  wave  frequency,  has  to  be 
slightly  greater  than  f. 

The  third  Indication  that  the  wave  frequencies  are  somewhat 
greater  than  f comes  from  attempts  to  least-squares  fit  sinusoidal 
waves  to  the  data.  At  any  pressure,  z,  it  is  possible  to  fit  sinu- 

ft 

soids  of  the  form 

ufU  (Zp.t^  = A(zn)  sin(<t>(zn)  + wt^  , 

vfU  (zn,ti)  = B(zn)  s1n(\p(zn)  + uti ) V-36 


over  a group  of  profiles,  t.  Is  the  time  of  a given  profile,  u and 
v are  east  and  north  velocity  components.  Then  the  quantities 


and 


U(2„)  = £ <um  (*„.»,)  - Uobs  (2„,t())’ 
v<*n>  = 5 ' “obs  (Tn'tt))! 


V-37 


V-38 
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give  an  estimate  of  the  residual  error  In  the  fit  at  depth  zn. 

These  fits  can  then  be  calculated  over  a range  of  depths,  zn.  The 
quantity 

D = t U(z  ) + V(zn)  V- 39 

z n 
n 

represents  a measure  of  the  total  residyal  error  of  the  fit  over 
the  time  and  depth  Interval  studied.  Attempts  to  do  such  fits  in 
the  main  thermocline,  where  the  strongest  signal  is  found,  indi- 
cate that  the  total  residual,  D,  is  minimized  if  the  period,  T, 
used  in  the  above  formulas  (w  = 2ir/T)  is  slightly  less  than  the 
inertial  period.  In  fact,  the  best  fit  occurs  if  we  choose  w so 
that  T a 22-23  hours.  The  inertial  period  is  25.56  hours  at  the 
latitude  of  the  profiler  time  series. 

We  have  argued  in  previous  chapters  that  the  near-inertial 
waves  are  propagating  their  energy  downward  from  the  surface.  Since 

A 

waves  with  frequencies  close  to  f have  group  velocity  vectors  that 
are  only  slightly  Inclined  to  the  horizontal,  it  is  at  first  glance 
surprising  that  these  waves  appear  to  make  significant  contributions 
to  the  energy  even  in  the  deep  water.  It  could  be  argued  that  those 
waves  which  have  the  maximum  (downward)  vertical  energy  flux  will 
be  able  to  penetrate  from  the  surface  to  the  deep  ocean  most  effi- 
ciently. This  vertical  energy  flux  is  given  by 

Fz  = Cg2  • F V-40 

where  Cgz  is  the  vertical  component  of  the  group  velocity  vector, 
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and  E Is  the  total  (kinetic  and  potential)  wave  energy.  It  Is  of 
interest  to  obtain  the  local  Cg2,  as  a function  of  z,  and  use  It, 
along  with  the  theoretical  energy  spectrum  of  the  previous  section, 
to  calculate  as  a function  of  frequency  and  vertical  wave  number. 

Because  the  derivation  of  the  Garrett-Munk  spectrum  assumes 
that  a WKB-type  description  of  the  wave  behavior  is  adequate,  we 
will  assume  In  what  follows  that  this  is  again  true.  The  equations 
to  be  solved  are  the  same  as  equations  B-6,  except  that  now  f)  Is  a 
function  of  depth.  The  equation  for  the  vertical  velocity,  w,  is 


♦ s’  <8’(Z) . . o 

dz2 


V-41 


with 


A.  A 

. . k2  + l2 
6 " aS2  - ?2 
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To  obtain  the  WKB  solution,  we  let 


w = e 


ejx(z)dz 


A„U) 


n=0  B 

A 

The  solutions  for  n = 0 and  n = 1 give,  respectively,  x and  AQ: 
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Then 


X - t 1 [NZ(z)  - S2]1/2 


A = 1 

o - ~ 1/4 


[N2  - oo2] 
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[82(z)  - w2] 
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where  the  region  of  Integration  Is  over  an  Interval  of  z for  which 

A ^ 

N2(z)  > u*.  From  continuity  ,we  get  u and  v: 


u » (+  wk  ? Ifl)^ 

v - (♦  IK?  * u)T)/~ri 

• Ml  fi\*  - 


- (o' 


f x(z)dz 


Z{Z2  - f*)1/2(k*  + ll)1/2 


e zo 
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The  buoyancy  Is 


+ B \ X dz 


b 3 6 

1w(f)2  - w2)1/4 

The  total  energy,  F,  is  given  by 


E » 1/2  ( < |u| 2 + | v | 2 + |w| 2 + 4±2  ) 


[5*  - ^)'^2 


+ f1 


to)2  + N2 


w2  - f2  w2(N2  - Z2)'1^  , 


where  the  angle  brackets  denote  an  average  over  time.  The  ratio 
of  vertical  kinetic  energy  plus  potential  energy  to  horizontal 
kinetic  energy  is 


vert.  K.E.  + P.E.  _ (N2  4 Z2)  [Z2  - f2) 
horlz.  K.E.  = (N2  - Z2)  (Z2  + f2) 


For  frequencies  close  to  ?,  almost  all  of  the  energy  appears  as 
horizontal  kinetic  energy.  The  vertical  component  of  the  group 


velocity  vector  at  depth  z is 

Cg,(m,u,z)  * & I ■ - & 


2 - f2)  (N2  - Z2) 


m(z)  Z (A2  - ?2)  . 


The  theoretical  spectrum  of  Garrett  and  Munk  has  already  been  given 
(equation  V-24).  This  spectrum  contains  only  horizontal  kinetic 
energy.  However,  equation  V-50  shows  that  for  frequencies  close  to 
t,  most  of  the  wave  energy  is  horizontal  kinetic  energy.  From 
equation  V-24,  the  energy  density  spectrum  is 

- 6c  <-£->  «L*J1 

*0  VMo  w3(w2  - ?2) ]/Z  * V-52 

The  vertical  energy  flux  for  waves  in  a small  Interval  of  frequen- 
cy-stretched vertical  wave  number  space  Am' Aw  is 

F = Cg  Hill  E(m',w)  Am'Aw  , V-53 

z z n0 

A A A A 

where  m'  and  w are  taken  to  be  in  the  middle  of  the  element  Am'Aw. 

A A 

The  vertical  group  velocity  can  be  expressed  in  terms  of  w and  m* 


by  using  equation  V-12 


(h  _ lu 2 _ r.2' 


* , — . » "o  (wz  - fz)  (N2  - w 

C,z  ( .*)  - — Vl  - V — 


F , . ac  A I*2  * f1)  2 - “2>  - f2 

M03  VMo  (N2  - f2) 


>V2  . 

! Am' Aw 


/V  A /\ 

We  note  that  for  w close  to  f,  Fz  is  almost  independent  of  z. 
From  the  WKB  approximation,  it  can  be  shown  that 


3F,  a 

sT  ’ 31  <C9Z  • E)  - o 


However,  since  we  are  only  considering  horizontal  kinetic  energy, 

8F z/3z  is  not  equal  to  zero  but  approaches  zero  as  u>  approaches  f. 

We  are  particularly  interested  in  examining  ?2  as  a function 
of  a),  for  fixed  m'.  The  frequency  dependence  of  Fz  is  given  by 

go) « ♦ n (y  - - S-) 

i>(H!  - ?!) 

c*  * f>)  (i!  - ?«)V20  - ^ > 

to"  , V- 57 

A A A /\ 

since  N » f.  If  w and  N are  made  non-dimensional  by 


then 

(i2  + 1)  (£2  - 1)1/2(1  - &L) 

f(w)  * fg(u)  = 

““  * V-59 

f(u>)  is  plotted  in  Figure  28,  with  N = 5f.  It  can  be  seen  that 

f(u>)  increases  rapidly  with  increasing  i until  w * 1.2,  and  then 

A ^ 

decreases.  This  implies  that,  for  fixed  m',  Fz  reaches  a maximum 
for  uj  * 1.2f.  For  an  inertial  period  of  25.56  hours,  the  vertical 
flux  of  horizontal  kinetic  energy  reaches  a maximum  (at  a given  m ' ) 
for  waves  with  periods  of  around  21  hours.  This  period  appears  to 
be  somewhat  shorter  than  that  which  might  have  been  expected.  (We 
recall  that  the  least-square  fit  indicated  that  the  dominant  waves 
in  the  profiles  seemed  to  have  periods  of  22-24  hours  in  the  main 
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thermocl ine. ) However,  It  is  significant  that  the  downward  energy 
flux  appears  to  have  its  maximum  value  close  to  the  local  inertial 
frequency  in  this  model. 

A certain  amount  of  inconsistency  is  inherent  in  using  a 
spectrum  such  as  the  Garrett-Munk  model  to  calculate  energy  flux. 

The  first  point  of  inconsistency,  of  course,  is  that  the  true  spec- 
trum will  not  have  an  infinity  at  the  inertial  frequency.  This 
problem  could  be  removed  by  assuming,  for  example,  that  there  should 
really  be  an  energy  cut-off  near  the  inertial  frequency;  that  is, 
instead  of  the  infinity  at  the  inertial  frequency,  the  spectrum 
could  be  altered  locally  (at  and  near  f)  so  that  a simple  peak  oc- 
curs at  the  inertial  frequency.  This  alteration  would  not  change 
the  flux  curve  of  Figure  28  appreciably.  In  fact,  it  would  force 
the  observed  flux  maximum  to  be  concentrated  even  more  about  per- 
iods of  21-22  hours.  The  second  problem  is  more  subtle,  and  in- 
volves the  question  of  why  one  would  expect  the  spectral  peak  to 
be  exactly  at  the  inertial  frequency.  In  a case  (such  as  the 
Garrett-Munk  model)  where  vertical  modes  have  been  assumed  to  have 
been  established  a priori,  it  seems  that  a peak  at  the  inertial 
frequency  would  not  be  inconsistent,  since  the  Waves  are  not  re- 
quired to  be  supported  by  energy  sources  located,  for  example,  at 
the  surface  or  bottom  boundaries.  However,  once  we  require  that 
the  waves  be  supported  by  a flux  of  energy  into  the  ocean  from  the 
surface  (or  bottom)  boundary  it  is  not  consistent  to  expect  that 
the  peak  in  the  observed  energy  spectrum  at  some  depth  will  be 
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found  at  that  frequency  (f)  where  the  vertical  energy  flux  becomes 
zero.  Both  of  the  above  points  argue  that  the  "Inertial  peak" 
should  actually  appear  at  frequencies  somewhat  greater  than  f. 

Figure  29  shows  a frequency  spectrum  of  total  horizontal  kinet- 
ic energy  calculated  from  data  obtained  by  a current  meter  at  1500 
meters  near  the  central  mooring  from  April  18,  1974  to  July  27,  1974 
(D.  Porter,  1975,  personal  communication).  (Note  that  this  is  one 
year  after  the  profiler  time  series,  but  is  at  approximately  the  same 
location  and  is  in  the  same  season  as  the  profiler  series.)  Figure 
29  indicates  that  the  center  of  the  inertial  "peak"  is  displaced 
toward  frequencies  greater  than  f (of  course,  this  peak  also  contains 
diurnal  tidal  energy).  The  center  of  this  wide  peak  is  in  the  period 
range  22-23  hours,  in  agreement  with  the  above  discussion. 

It  Is  possible  to  use  the  above  expression  for  the  vertical 
group  speed  (equation  V-54)  along  with  the  observed  stretched  verti- 
cal wave  number  spectrum  (Figure  16)  to  calculate  a rough  estimate 
of  the  vertical  energy  flux  (downward  and  net),  integrated  over  all 
vertical  wave  numbers.  In  the  absence  of  accurate  frequency  infor- 
mation, we  have  assumed  that  the  clockwise  and  counterclockwise 
spectra  are  dominated  by  waves  with  periods  of  the  calculated  max- 
imum flux  (or  of  the  center  of  the  near-inertial  "peak"  in  Figure 
29);  that  is,  periods  of  22-23  hours.  Using  this  period  range  in 
the  expression  for  the  vertical  group  speed,  and  assuming  that 
roughly  half  the  energy  in  the  clockwise  and  counterclockwise  energy 
estimate  for  the  smallest  two  or  three  stretched  vertical  wave 
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Figure  29.  Freauency  spectrum  of  horizontal  velocity  (east  and  north 
components)  at  1500  m near  the  central  mooring.  The  data 
on  which  the  computations  were  performed  occupied  the 
time  interval  April  18,  1974  to  July  27,  1974.  Estimates 
have  been  averaged  over  4 adjacent  frequency  bands,  f is 
the  local  inertial  frequency.  1 T.U.  = 1 hour. 
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number  estimates  Is  contributed  by  the  semidiurnal  tide  (see  Chap- 
ter V-b),  we  obtain  estimates  of  the  total  downward  wave  energy 
flux  of  0. 6-1.0  erg/cm2/sec,  and  of  the  net  downward  wave  energy 
flux  of  0.2  - 0.3  erg/cm2/sec.  The  estimate  of  the  total  down- 
ward flux  is  consistent  with  other  estimates  of  vertical  energy 
flux  by  internal  waves  (Bell,  1975).  We  have  suggested  in  pre- 
vious chapters  that  a bottom  boundary  layer  could  be  a significant 
energy  sink  for  near-inertial  waves.  An  order-of-magnitude  cal- 
culation indicates  that  a net  downward  energy  flux  of  0.2  - 0.3 
erg/cm2/sec  could  be  matched  by  energy  loss  in  a boundary  layer. 

A measure  of  this  possible  energy  loss  is 

L = 4Kvp((az)2  + AZ  er9/cn,2/sec  v"60 

(Yih,  1969) 

where  Kv  is  the  coefficient  of  vertical  eddy  viscosity,  p is  den- 
sity, u and  v are  horizontal  velocity  components,  and  Az  is  a scale 
depth  of  the  boundary  layer.  Choosing  values  for  Au2  + Av2  and 
Az  from  Figure  8 (Au2  + Av2  - 4 cm2/sec2,  Az  ~ 5 x 103  cm),  we  get 
L - (0.3  x 10-2)  Ky  erg/cm2/sec.  For  Ky  * 100  cra2/sec,  L - 0.3 
erg/cm2/sec.  Bell  (1975)  quotes  a value  for  Ky  of  ~ 10  cm2/sec  in 
the  interior  but  notes  that  Ky  could  be  at  least  an  order  of  mag- 
nitude larger  in  a turbulent  bottom  boundary  layer.  Thus,  the 
value  of  Ky  needed  to  make  L comparable  to  the  estimated  net  down- 


Chapter  VI  Some  Observations  over  Rough  Topography 

In  previous  chapters  we  have  presented  an  analysis  of  data 
obtained  over  smooth  topoqraphy  (at  the  MODE  central  mooring). 

Some  of  the  calculations  presented,  such  as  reflection  coefficients, 
depend  on  the  assumption  that  the  bottom  is  not  acting  as  a wave 
source.  Over  the  smooth  topography  this  does  not  appear  to  be  an 
unreasonable  assumption.  However,  it  is  well  known  from  theoretical 
studies  that  a rough  bottom  may  interact  with  either  a mean  flow, 
producing  lee  waves,  or  with  an  oscillatory  flow  such  as  a tide. 

Bell  (1975)  has  discussed  the  possible  generation  of  internal  waves 
in  the  open  ocean  by  rough  topography.  In  an  earlier  work, 
Hendershott  (1964)  analyzed  a case  where  inertial  waves  could  be 
generated  by  the  oscillatory  flow  of  a diurnal -tide  component  over 
a rough  bottom.  His  study  indicated  that  strong  inertial  wave 
currents  might  be  generated  at  latitudes  slightly  to  the  south  of 
the  critical  latitude  (the  latitude  at  which  the  tidal  and  .nertial 
periods  are  equal).  The  critical  latitude  for  the  24-hour  diurnal 
tide  would  be  30°  N,  for  example.  His  theory  actually  predicts 
that  there  should  be  a somewhat  complicated  band  structure  of  the 
inertial  energy  south  of  the  critical  latitude,  with  a primary  peak 
in  the  inertial  energy  just  south  of  the  critical  latitude,  and 
secondary  maxima  farther  south. 
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It  Is  Important  to  realize  that  there  are  numerous  tidal 


lines  that  are  called  diurnal.  Although  the  critical  latitude  for 


the  24-hour  tide  is  30°  N,  the  critical  latitude  for  the  24-hour 


50-minute  tide  moves  down  to  about  28.9°  N.  This  means  that  the 


actual  energy  levels  of  the  inertial  waves  (generated  by  topography) 


might  not  depend  as  strongly  on  latitude  as  is  predicted  by  the 


theory.  In  any  case,  the  fact  that  the  major  diurnal  tidal  compo- 


nents at  28°  N have  frequencies  within  the  internal  wave  frequency 


range,  but  near  the  Inertial  frequency.  Indicates  that  it  would  be 


worthwhile  to  look  for  near-inertial  energy  propagating  away  from 


the  bottom  in  the  rough  topography  region. 


During  the  MODE  experiment,  a set  of  profiles  was  obtained 


around  a ridge  in  the  rough  topography  region  east  of  the  central 


mooring.  These  drops  are  indicated  by  a dagger  (t)  in  Appendix  A. 


The  center  of  the  ridge  is  located  at  approximately  28°  04'  N, 


68°  31'  W.  The  main  part  of  the  ridge  is  oriented  along  the  axis 


NNE-SSW  and  rises  about  700-800  meters  above  the  adjacent,  smoother 


part  of  the  ocean  floor.  Because  this  set  was  obtained  in  order  to 


make  a survey  of  the  velocity  structure  around  the  ridge  it  is  not 


a time  series  at  one  location.  In  fact,  the  set  is  made  up  of  a 


series  of  pairs  of  profiles,  except  in  a few  cases  where  only  single 


profiles  were  made.  Each  pair  was  obtained  at  one  location  around 


the  ridge,  and  the  second  profile  was  usually  taken  at  about  one 


half  of  an  inertial  period  after  the  first  profile.  Because  a time 


series  of  profiles  is  not  available  over  the  rough  topography,  it 
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is  difficult  to  accurately  remove  the  low-frequency  shear  profile. 
For  this  reason,  the  results  of  calculating  clockwise  and  counter- 
clockwise vertical  wave  number  spectra  are  more  difficult  to  inter- 
pret than  are  the  corresponding  results  over  the  smooth  topography. 
However,  some  significant  features  in  these  rough  topography  pro- 
files can  be  seen  visually;  that  is,  without  recourse  to  spectrum 
calculations. 

Figures  30  and  31,  for  example,  show  hodograph  plots  (plots 
of  the  horizontal  velocity  vector  with  depth)  for  profiles  190D  and 
192D  between  3000  dbar  and  4500  dbar.  These  two  profiles  were  made 
at  a time  laq  of  one  half  of  an  Inertial  period.  The  original  pro- 
files have  been  smoothed  by  applying  a running  150  dbar  triangular 
filter  to  each  velocity  component.  This  makes  it  easier  to  vis- 
ually follow  the  behavior  of  the  longer  wave  length  (higher  energy) 
waves  in  the  profiles.  Smoothed  points  are  plotted  every  10  dbar, 
and  depth  is  marked  every  100  dbar.  The  most  striking  feature  of 
these  hodographs  is  that  they  show  wave-like  structures  for  which 
the  horizontal  velocity  vector  is  clearly  rotating  counterclockwise 
with  increasing  depth.  The  structure  is  clearly  time-dependent,  as 
can  be  seen  by  comparing  the  hodographs  for  profiles  1900  and  192D. 
Comparison  of  the  hodographs  appears  to  indicate  that  these  waves 
have  periods  near  the  inertial  period.  We  have  already  pointed  out 
in  Chapter  III  that  such  counterclockwise  polarization  indicates 
that  these  waves  are  propagating  energy  upward,  away  from  the  ridge. 


Figure  30.  Smoothed  hodograph  of  profile  190D  between  3000  dbar 
and  4500  dbar.  Depths  are  in  hundreds  of  decibars. 
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Although  spectral  estimates  obtained  from  only  two  profiles 
are  probably  not  very  stable,  it  should  be  mentioned  that  vertical 
wave  number  spectra  (without  stretching  the  vertical  coordinate) 
obtained  in  the  deep  water  (2000  dbar  and  deeper)  for  drops  190D, 
191D,  192D  and  1930. show  a consistent  dominance  of  counterclockwise 
energy  over  clockwise  energy.  These  four  drops  were  made  at  the 
southwest  end  of  the  ridge. 

In  contrast,  deeper  water  (greater  than  2000  dbar)  spectra 
and  hodogr^ihs  for  a drop  pair,  1950  and  197D,  that  was  made  north- 
east of  the  ridge,  do  not  show  any  strong  tendency  for  counter- 
clockwise rotation  with  depth.  If  anything,  these  drops  show  a 
tendency  for  the  velocity  vector  to  rotate  clockwise  with  depth  in 
the  same  manner  as  was  observed  over  the  smooth  topography.  Swallow 
floats  located  near  the  ridge  when  these  profiles  were  made  indi- 
cated that  the  mean  flow  (at  4000  dbar)  was  to  the  southwest. 

Since  the  pair  195D  and  197D  was  located  to  the  northeast  of  the 
ridge  while  1900,  1 9 1 D , 192D  and  193D  were  to  the  southwest,  we 
might  speculate  that  if  the  ridge  was  generating  near  Inertial  waves, 
these  waves  would  be  advected  to  the  southwest  and  therefore  would 
be  seen  primarily  by  profiles  to  the  southwest  of  the  ridge.  Al- 
though it  is  difficult  to  draw  general  conclusions  from  a small  set 
of  profiles,  the  above  observations  in  the  deep  water  around  the 
ridge  can  be  summarized  as  follows:  profiles  made  to  the  southwest 

of  the  ridge  show  clear  evidence  of  counterclockwise  polarization 
with  depth;  profiles  made  near  the  ridge  crest  (such  as  189D,  for 


example)  do  not  appear  to  show  a dominance  of  either  polarization; 
two  profiles  (195D  and  197D)  made  over  a smoother  region  to  the 
northeast  of  the  ridge  appear  to  have  reverted  to  the  clockwise 
polarization  with  depth  characteristic  of  the  smooth  topography 
profiles,  at  least  for  the  longer  vertical  wave  lengths,  which  are 
most  energetic. 

Bell  (1975)  has  calculated  a value  for  the  possible  upward 

energy  flux  of  internal  waves  generated  by  rough  topography  of 

~ 1 erg/cm2/sec.  This  value  is  comparable  to  the  total  downward 

flux  in  the  smooth  topography  region,  as  estimated  in  Chapter  V. 

If  the  average  vertical  energy  flux  values  for  bottom  topography 

and  sea  surface  generation  are  comparable,  it  might  be  expected 

that  energy  input  at  the  bottom  could  sometimes  exceed  energy  input 

« 

at  the  surface  (especially  since  botn  input  processes  are  probably 
intermittent  in  time).  This  would  result  in  waves  polarized 
counterclockwise  with  increasing  depth,  at  least  in  the  deep  water 
over  the  rough  topography.  (We  have  also  examined  hodographs  of 
the  rough  topography  profiles  in  the  upper  half  of  the  water 
column.  In  contrast  to  the  plots  in  the  lower  half,  these  hodo- 
graphs do  not  show  a clear  dominance  of  one  polarization  over  the 
other,  indicating  that  in  the  upper  half  of  the  water  column  a 
complicated  addition  of  the  top  and  bottom  energy  sources  may  be 
taking  place.) 


Chapter  VII  Conclusions,  a Further  Discussion  of  Some  of  the  Data 
and  Recommendations  for  Future  Experiments 

The  previous  chapters  have  been  devoted  to  an  analysis  of 
internal  wave  data  obtained  by  vertical  profiling.  It  Is  now 
worthwhile  to  summarize  the  significant  results  of  this  analysis 
and  place  these  results  in  the  general  framework  of  internal  wave 
theory . 

One  of  the  more  fruitful  results  of  this  analysis  has  been 
the  ability  to  relate  the  observed  characteristics  of  wave  polar- 
ization and  vertical  phase  propagation  to  the  locations  of  energy 
sources  for  these  waves;  that  is,  to  the  possible  existence  of 
internal  wave  energy  sources  at  the  top  and  bottom  boundaries  of 
the  ocean.  As  Muller  and  Olbers  (1975)  have  pointed  out,  a two- 
sided,  symmetric,  vertical  wave  number  spectrum  (equal  energy  in 
positive  and  negative  vertical  wave  numbers)  is  not  consistent  with 
the  existence  of  net  boundary  sources  at  the  surface  and  bottom. 

A major  cause  of  difficulty  with  internal  wave  observations  until 
now  has  been  the  inability  to  observe  asymmetry  of  internal  wave 
vertical  wave  number  spectra  (or,  in  fact,  vertical  wave  number 
spectra)  of  horizontal  velocity  with  the  techniques  that  have  been 
employed.  Time  series  of  horizontal  velocity  obtained  from  current 
meters  at  fixed  locations  in  space  cannot  distinguish  between  an 
internal  wave  having  an  upward  or  a downward  phase  propagation, 
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since  the  horizontal  velocity  vector  rotates  In  the  same  sense  for 
either  wave  (clockwise  in  the  northern  hemisphere).  The  tech- 
nique, employed  In  Chapter  II,  of  decomposing  vertical  profiles 
into  clockwise  and  counterclockwise  energy  contributions  is  a 
useful  method  for  observing  asymmetry  in  the  internal  wave  vertical 
wave  number  spectrum  and  for  inferring  the  direction  of  net  verti- 
cal energy  flux  from  this  asymmetry.  A similar  Instance  where 
vertical  wave  polarization  has  been  used  to  infer  the  direction  of 
vertical  energy  propagation  In  the  atmosphere  has  already  been 
pointed  out  (Hines,  1966).  A case  where  vertical  polarization  of 
tidal  currents  was  apparently  observed  In  an  analysis  of  current 
records  obtained  from  current  meters  on  moorings  has  been  reported 
by  Titov  and  Fomin  (1971).  We  have  mentioned  in  Chapter  III,  how- 
ever, that  the  ability  to  identify  clockwise  and  counterclockwise 
energy  with  waves  having  upward  and  downward  wave  number  vector 
components  (that  Is,  with  the  two-sided  vertical  wave  number  spec- 
trum) requires  that  the  observed  energy  be  primarily  contributed 
by  waves  with  frequencies  near  the  local  inertial  frequency.  This 
is  true  because  individual  waves  are  nearly  circularly  polarized  If 
their  frequencies  are  near  the  local  inertial  frequency.  If  a 
single  wave  with  a frequency  much  greater  than  inertial  were  pres- 
ent, for  example,  this  wave  would  provide  energy  contributions  to 
both  the  clockwise  and  counterclockwise  vertical  wave  number  spec- 
tra, since  the  wave  would  be  elliptically  polarized  in  the  vertical. 
The  clockwise  energy  would  still  be  greater  than  the  counterclockwise 
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energy  If  the  w a/e  were  propagating  energy  downward.  But  the 
true,  two-sided  vertical  wave  number  spectrum  of  this  wave  would 
contain  energy  only  on  one  side  of  the  spectrum  (at  the  positive 
or  negative  vertical  wave  number  of  the  wave,  depending  on  which 
way  the  phase  of  the  wave  is  propagating  in  the  vertical). 

Fortunately,  the  data  presented  and  analyzed  in  previous 
chapters  indicates  that  the  two  dominant  sources  of  energy  in  the 
profiles  are  near-inertial  waves  (at  all  vertical  wave  numbers), 
and  the  semidiurnal  tide  (at  the  smallest  vertical  wave  numbers). 
Inertial  and  diurnal  tidal  contributions  to  the  profiles  cannot, 
however,  be  resolved  because  of  the  short  length  of  time  over  which 
the  profile  time  series  was  obtained.  The  frequency  spectrum  at 
1500  meters  (Figure  29)  shows  that  a large  part  of  the  energy  ob- 
served at  this  depth  in  the  internal  wave  band  comes  from  frequen- 
cies near  the  inertial -diurnal  tidal  frequency,  and  from  the  semi- 
diurnal tide.  The  observed  "mirror-imaging"  of  profiles  (Figure 
3)  and  the  calculations  of  dropped,  lagged,  rotary  coherence  (Fig- 
ures 24  to  27),  support  this. 

The  observation  that  much  of  the  variability  in  the  profiles 
was  caused  by  the  near-inertial  waves  (and  the  semidiurnal  tide  at 
the  largest  vertical  wave  lengths)  led  to  efforts  to  evaluate  quan- 
titatively the  total  and  net  downward  energy  flux  over  the  smooth 
topography.  Making  assumptions  about  the  amount  of  energy  contri- 
buted by  the  semidiurnal  tide  at  the  smallest  vertical  wave  numbers, 
we  obtained  values  for  the  total  vertical  energy  flux  downward  of 
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- .6-1.0  erg/cm2/sec,  and  for  the  net  vertical  energy  flux  down- 
ward of  -.2  - .3  erg/cm2/sec.  Although  both  estimates  are  subject 
to  uncertainty  because  of  the  presence  of  the  semidiurnal  tide, 
neither  estimate  seems  unreasonable  (see  Chapter  V-c).  In  partic- 
ular, It  appears  that  the  net  downward  energy  flux  can  be  balanced 
by  a rough  estimate  of  the  rate  of  energy  loss  per  cm2  in  the  bot- 
tom boundary  layer  (if  we  assume,  of  course,  that  the  observed  de= * 
crease  of  mean  perturbation  horizontal  kinetic  energy  near  the 
bottom  is  actually  caused  by  the  presence  of  a boundary  layer). 

In  relation  to  the  contour  plots  (Figures  6 and  7),  a fur- 
ther statement  can  be  made  which,  although  somewhat  obvious,  is  no 
less  important.  The  observed  relations  among  wave  period,  direc- 
tion of  vertical  phase  propagation,  and  vertical  wave  polarization 
In  Figures  6 and  7 constitute  at  least  a partial  verification  of 
the  dispersion  relation  for  Internal  waves  (those  near  the  inertial 
frequency).  Other  observations  and  analyses  of  data  (from  current 
meter  records,  for  example)  have  assumed  that  the  observed  oscil- 
lations between  the  inertial  and  Brunt-Vaisala  frequencies  obey 
the  dispersion  relation  for  internal  waves.  In  the  absence  of  clear 
evidence  of  the  propagation  of  oscillations  in  this  frequency 
range,  other  investigators  (Webster,  1969,  for  example)  have  ar- 
gued that  some  other  process  (such  as  turbulence)  could  explain 
' the  observed  oscillations.  Although  turbulence  is  difficult  to 

*.  define,  one  of  its  characteristics  Is  that  It  does  not  follow  a 

" well-defined  dispersion  relation.  In  this  sense,  then,  the  results 

r 
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of  Figures  6 and  7 clearly  arque  for  internal  waves  as  being  the 
cause  of  the  observed  variability  in  the  profiles,  at  least  for 
the  near-inertial  waves. 

Another  significant  result  brought  out  by  this  analysis  is 
the  fact  that  both  profiler  data  and  other,  independent  observa- 
tions (Figure  29)  have  shown  that  the  energy  concentration  in  what 
one  might  call  the  "inertial  peak"  is,  first,  not  concentrated  in 
a well-defined  peak  (as  is  the  semidiurnal  tide,  for  example)  and 
second,  is  not  concentrated  at  exactly  the  inertial  peak  but  is 
rather  centered  at  a frequency  about  20*  greater  than  the  local 
inertial  frequency.  Some  previous  observations  have  also  shown  that 
the  Inertial  "peak"  is  quite  often  displaced  toward  frequencies 
somewhat  greater  than  the  local  inertial  frequency.  Webster  (1968), 
for  example,  shows  a frequency  spectrum  of  current  at  38°  28.8'  N, 
70°  00.5'  W,  in  3300  meter  water,  with  the  observations  made  about 
30  meters  above  the  bottom.  This  spectrum  has  a clear  concentra- 
tion of  energy  near  the  local  inertial  freouency,  but  the  center 
of  this  peak  occurs  at  a frequency  about  5%  higher  than  the  local 
inertial  frequency.  White  (1972)  considered  temperature  data  ob- 
tained in  the  upper  several  hundred  meters  a;  a site  in  the  Pacific. 
His  temperature  spectra  always  show  an  inertial  "peak"  at  frequen- 
cies up  to  30%  higher  than  the  local  inertial  frequency. 

Various  attempts  have  been  made  to  explain  this  shift  theor- 
etically. Munk  and  Phillips  (1968),  for  instance,  argued  that 
internal  waves  generated  to  the  south  of  a given  latitude  of 


observation  (in  the  northern  hemisphere)  would  propagate  northward, 
and  that  as  they  propagated  northward  their  freguencies  would  ap- 
proach the  local  inertial  frequency.  A test  of  this  idea  by  the 
use  of  profile  data  would  require  better  information  on  the  direc- 
tion of  horizontal  phase  propagation  for  the  observed  waves.  We 
will  return  to  this  question  later,  when  the  possible  influence  of 
the  mean  vertical  shear  on  the  propagation  of  the  waves  seen  in 
the  profiles  is  discussed. 

White  (1972)  pointed  out  that  the  observed  frequency  shift 
of  the  inertial  waves  could  be  explained  as  the  result  of  a 
Doppler  shift  due  to  the  mean  flow.  Because  he  always  observed  a 
shift  toward  higher  frequencies,  he  was  forced  to  conclude  that 
the  near-inertial  waves  he  observed  were  always  propagating  in  the 
direction  of  the  mean  flow.  The  data  presented  in  previous  chap- 
ters, however,  show  that  it  is  unlikely  that  Doppler  shifting  could 
be  the  cause  of  the  observed  frequency  shift.  The  frequency  spec- 
trum at  1500  meters  (Figure  29)  shows  that  there  is  an  apprecia- 
ble (20%)  shift  of  the  "inertial  peak"  toward  higher  frequencies. 
But  it  is  known  from  a calculation  of  the  average  flow  observed 
by  a current  meter  at  1500  meters  near  the  location  of  the  five- 
day  time  series  (Chapter  II)  that  the  mean  flow  at  that  depth, 
when  the  time  series  was  made,  was  small  (<  3 cm/sec).  It  does  not 
seem  that  such  a small  mean  flow  could  account  for  the  observed 
frequency  shift.  This  does  not  necessarily  negate  White's  conclu- 
sion, since  he  was  observing  at  a point  where  the  mean  flow  was 
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about  five  times  that  observed  at  1 bOO  meters  in  the  present  exper- 
iment. It  would  seem  fortuitous,  however,  that  the  near-inertial 
waves  would  always  propagate  in  the  direction  of  the  mean  flow, 
especially  if,  as  the  present  data  indicate,  they  are  geing  gen- 
erated by  atmospheric  forcing  at  the  ocean  surface.  In  fact,  there 
is  a suggestion  in  the  data  (discussed  below)  which  indicates  that 
at  least  the  more  energetic  feature*  in  the  profiles  have  a compo- 
nent of  propagation  to  the  north.  Since  the  mean  flow,  at  least  to 
about  1500  meters,  was  to  the  south  during  the  present  experiment, 
this  would  imply  that,  if  anything,  the  inertial  peak  would  be 
shifted  to  lower  frequencies. 

The  results  of  the  profiler  observations  now  allow  us  to 
propose  another  mechanism  to  explain  why  the  observed  "inertial 
peak"  frequency  shift  almost  always  appears  to  be  toward  higher 
frequencies.  This  explanation  is  that  the  observed  waves  are 
clearly  propagating  in  the  vertical.  The  contour  plots  (Figures 
6 and  7)  show  that  the  dominant  structures  associated  with  the 
near-inertial  waves  are  moving  upward  in  time,  corresponding  to 
downward  energy  propagation.  (Of  course,  these  cor  tour  plots  are 
influenced  by  the  mean  flow.  However,  we  have  already  pointed  out 
that  around  1500  meters  the  mean  flow  is  small,  and  yet  the  great- 
est speed  of  upward  phase  propagation  in  these  contour  plots  is 
found  in  the  1500  dbar  to  2500  dbar  ranae.)  In  order  to  propagate 
vertically,  these  waves  must  have  . requencies  greater  than  the 
local  inertial  frequency.  The  calculators  of  Chapter  V-c  showed 
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that,  using  an  assumed  form  of  the  internal  wave  energy  spectrum, 
the  maximum  vertical  flux  of  horizontal  kinetic  energy  occurs  at 
frequencies  about  20%  higher  than  the  inertial  frequency.  We 
pointed  out  that,  once  the  energy  sources  are  placed  at  the  surface 
(or  bottom)  boundary,  it  is  inconsistent  to  have  the  peak  in  the 
energy  spectrum  appear  exactly  at  the  inertial  frequency  when  waves 
of  this  frequency  cannot  propagate  into  the  interior.  It  is  con- 
sistent, however,  to  have  an  "inertial  peak"  that  is  shifted  to 
higher  frequencies.  It  seems  likely  that  the  waves  which,  when 
generated  at  the  surface,  have  the  greatest  downward  energy  flux 
would  be  most  able  to  penetrate  through  the  total  water  column. 
(They  would  presumably  be  able  to  travel  a greater  distance  in  the 
vertical  before  losing  an  appreciable  fraction  of  their  energy  due 
to  friction,  non-linear  interactions,  shear  interactions,  or  some 
other  cause.) 

We  have  argued  that  a significant  part  of  the  energy  put  in 
by  the  atmosphere  is  lost  in  a bottom  boundary  layer,  at  least  for 
the  near  inertial  waves.  Calculations  of  the  net  downward  energy 
flux  and  the  rate  of  energy  loss  in  what  we  are  interpreting  in 
Figure  8 to  be  a bottom  boundary  layer  show  that  (in  terms  of  a 
rough  calculation)  the  energy  flux  into  the  boundary  layer  could  be 
accounted  for  by  the  energy  loss  in  the  boundary  layer.  It  is  cer- 
tainly true  that  other  factors  can  act  to  remove  energy  from  the 
wave  field.  For  example,  it  is  clear  that  the  observed  waves  tend 


to  increase  in  amplitude  and  decrease  in  vertical  wave  length  in 


the  thermocline,  where  the  Brunt-Vaisala  frequency  is  large.  For 
waves  shorter  than  some  critical  value,  the  shear  caused  by  these 
waves  will  be  large  enough  so  that  the  Richardson  number  would  be 
less  than  1/4,  and  the  flow  caused  by  the  wave  might  become  unstable. 
This  would  extract  energy  from  the  waves  (possibly  putting  that 
energy  into  turbulence).  Such  a process  would  be  most  active  at 
the  high-vertical-wave  number  end  of  the  vertical  wave  number  spec- 
trum. It  is  also  reasonable  to  expect  that  frictional  losses  will 
be  important  throughout  the  water  column  (not  just  in  a bottom 
boundary  layer)  for  waves  with  short  vertical  wave  lengths.  As  the 
vertical  wave  length  decreases,  the  frictional  term,  v 32v/dz2,  in 
the  horizontal  momentum  equation  (Chapter  IV)  for  near-inertial 
waves  grows.  In  terms  of  the  boundary  layer  discussed  in  Chapter 
IV,  we  can  say  that  (assuming  the  frequency  is  fixed  and  close  to 
inertial)  as  the  vertical  wave  length  decreases,  the  horizontal 
wave  length  will  also  decrease.  This  will  cause  the  parameter 
R = vkz/f  of  Chapter  IV  to  increase.  When  this  occurs,  frictional 
losses  are  no  longer  confined  to  a thin  bottom  layer,  but  rather 
become  important  in  the  interior  of  the  water  column  as  well.  This 
indicates  that  the  argument  that  near-inertial  waves  are  losing  most 
of  their  energy  in  a bottom  boundary  layer  is  probably  only  appli- 
cable for  near-inertial  waves  with  small  vertical  wave  numbers. 

This  process  of  energy  loss  due  to  friction  in  the  interior  of  the 
water  column  would  also  be  expected  to  be  most  active  at  the  high- 
vertical  -wave  number  end  of  the  vertical  wave  number  spectrum. 


Another  cause  of  energy  loss  by  these  waves  would  be  by 
interaction  of  the  waves  with  the  mean  flow.  In  processing  the 
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profiler  data  we  chose  to  normalize  the  profiles  by  the  mean  Brunt- 
Vaisala  profile.  This  type  of  normalization  was  used  because  it 
made  the  calculations  easier  to  compare  with  existing  theories. 
However,  it  is  well  known  (Phillips,  1966)  that  a vertical  shear 
in  the  mean  horizontal  velocity  profile  can  also  cause  changes  in 
the  amplitude  and  wave  length  of  internal  waves  as  these  waves 
propagate  through  the  shear  profile.  Perhaps  the  greatest  diffi- 
culty in  evaluating  the  influence  of  the  shear  on  the  profiles  is 
lack  of  data  on  the  direction  of  horizontal  phase  propagation  of 
the  waves  relative  to  the  mean  shear  flow.  We  have  made  several 
attempts  to  determine  the  direction  of  horizontal  phase  propaga- 
tion. One  attempt  involved  calculating  the  coherence  between  temp- 
erature and  horizontal  velocity  components  as  obtained  by  the  pro- 
filer in  individual  drops.  In  theory,  observed  phase  differences 
between  temperature  and  east  velocity  component  and  temperature  and 
north  velocity  component,  in  the  case  of  a single  dominant  wave, 
would  allow  one  to  determine  the  direction  of  horizontal  phase 
propagation.  In  practice,  the  observed  levels  of  coherence  between 
temperature  and  the  two  velocity  components  were  too  low  to  infer 
anything  about  horizontal  phase  propagation  direction.  Using  the 
observed  energy  levels  of  the  near-inertial  waves  in  the  main  ther- 
mocline  and  the  mean  vertical  temperature  gradient  observed  there, 
a rough  calculation  indicates  that  a near-inertial  wave  with  a 
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frequency  20%  greater  than  the  local  Inertial  frequency  would  pro- 
duce a temperature  signal  that  would  be  measurable  in  the  profiles, 
given  the  resolution  level  attainable  by  the  EMVP  temperature  meas- 
urements. The  fact  that  the  coherence  between  temperature  and 
velocity  measurements  is  so  low  may  be  due  to  the  presence  of 
temperature  variations  caused  by  waves  with  frequencies  much  greater 
than  the  inertial  frequency  (the  semidiurnal  tide,  for  example). 

In  order  to  use  the  temperature  profiles  to  determine  horizontal 
propagation  direction,  better  frequency  resolution  of  the  profile 
data  must  be  obtained.  Comparison  of  drops  235U  and  2361)  (Figure 
18),  both  visually  and  by  doing  lagged  correlations  between  the 
east  components  and  north  components  of  the  two  profiles,  suggests 
that  strong  features  which  can  be  identified  between  the  two  pro- 
files occur  at  a slightly  greater  depth  in  profile  236u  (4.8  kilo- 
meters to  the  north  of  235U).  if  this  is  true  it  implies  that 
these  waves  have  a component  of  propagation  toward  the  north. 

Taking  into  account  the  fact  that  the  mean  shear  flow  is  primarily 
in  the  north-south  direction  (Figure  7),  it  would  be  expected 
(Frankignoul  ,1972)  that  a wave  group  with  a component  of  propaga- 
tion to  the  north  and  propagating  energy  downward  would  be  acted 
upon  by  the  shear  in  such  a way  that  the  vertical  wave  length  would 
decrease  with  increasing  depth.  In  fact,  if  the  wave  is  thought  to 
be  propagating  directly  north  (against  the  major  part  of  the  mean 
shear),  an  analysis  of  the  WKB  expressions  for  the  change  in  verti- 
cal wave  length  brought  about  by  changing  N and  the  presence  of  the 
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shear  suggests  that  the  two  Influences  could  he  of  comparable  mag- 
nitude. However,  Fiqure  3 shows  that  the  most  energetic  waves  have 
wave  lengths  in  the  deep  water  that  are  longer  relative  to  observed 
wave  lengths  In  the  thermocline.  We  have  also  tried  to  use  the 
stretched,  normalized  profiles  to  see  if  there  is  a shift  of  ener- 
getic spectral  peaks  in  the  upper  versus  lower  parts  of  the 
(stretched)  water  column.  We  chose  a set  of  stretched  profiles 
(230D  through  2330)  which  contained  a distinct  peak  at  about 
130  sdb.  Vertical  wave  number  spectra  for  the  upper  and  lower 
halves  of  these  profiles  were  then  calculated  separately.  On  the 
assumption  that  the  stretching  of  these  profiles  has  removed 
changes  in  vertical  wave  number  and  wave  amplitude  due  to  variations 
of  N,  any  change  in  the  position  of  the  peak  in  the  upper  half  rela- 
tive to  the  lower  half  might  reflect  the  influence  of  shear.  The 
result  showed  that  there  is  a slight  shift  toward  higher  stretched 
vertical  wave  number  in  the  lower  half  relative  to  the  upper  half, 
but  it  is  very  small  ( < 5 sdb).  One  might  expect  a shift  in  this 
direction  if  the  waves  have  a component  of  horizontal  propagation 
to  the  north.  The  fact  that  the  shift  is  so  small  might  indicate 
that  the  waves  are  propagating  primarily  east-west,  with  a small 
propagation  component  to  the  north. 

The  major  conclusion  of  the  above  paragraphs  is  that  a good 
evaluation  of  the  Influence  of  the  mean  shear  on  wave  propagation 
awaits  better  estimates  of  horizontal  phase  propagation  directions. 
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A final  Important  result  of  the  present  work  Is  that  wave 
structures  have  been  observed  over  the  rough  topography  which  have 
polarization  appropriate  to  upward  energy  propagation.  As  pointed 
out  in  Chapter  VI,  there  is  a theoretical  basis  for  believing  that 
rough  bottom  topography  at  the  latitude  of  this  experiment  ( ~ 28°N) 
could  generate  near-inertial  waves  that  must  propagate  their  energy 
upward,  away  from  the  bottom.  Bell  (1975)  has  calculated  that 
rough  topography  could  support  an  upward  energy  flux  of  ~ 1 erg/cm2/ 
sec.  If  we  accept  our  estimate  of  .6-1.0  erg/cm2/sec  downward 
for  the  energy  flux  to  internal  waves  input  by  some  (as  yet  only 
poorly  understood)  atmospheric  forcing,  then  we  see  that  rough 
topography  generation  of  upward  energy  flux  can  be  comparable  to 
atmospheric  generation  of  downward  energy  flux.  Some  data  from 
current  meters  at  4000  meters  in  the  MODE  region  (D.  Porter,  1975, 
personal  communication)  indicates  that,  in  the  band  of  periods  from 
37  hours  to  14  hours  (encompassing  the  inertial  period),  the  total 
energy  in  the  frequency  spectra  obtained  at  4000  meters  depth  over 
the  smooth  topography  is  significantly  less  (at  least  by  a factor 
of  2)  than  the  total  energy  in  the  above  frequency  band  in  spectra 
obtained  at  the  same  depth  over  the  rough  topoqraphy.  The  above 
observations  suggest  that,  at  least  in  regions  of  rough  topography, 
one  might  sometimes  expect  to  see  a convergence  of  energy  toward 
the  thermocline,  due  to  the  combined  effects  of  atmospheric  and  bot- 
tom forcing.  According  to  Frankignoul  and  Strait  (1972),  such  a 
convergence  in  the  near-inertial  frequency  band  can  be  seen  in 
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records  obtained  at  Site  0 (39°  20'  N,  70°  00'  W).  We  again  empha- 
size that  we  have  not  seen  any  evidence  of  dominant  upward  energy 
propagation  (i.e.  counterclockwise  polarization  with  depth)  in  the 
smooth  topography  time  series. 

The  above  discussion  makes  clear  some  of  the  problems  asso- 
ciated with  interpreting  the  present  analysis.  These  problems,  in 
turn,  suggest  some  future  experiments  which  could  be  performed  In 
order  to  clarify  the  influence  of  various  processes  on  internal 
wave  propagation.  We  would  now  like  to  propose  several  relevant 
experiments.  First,  a time  series,  similar  to  the  one  obtained 
over  smooth  topography,  should  be  performed  over  a rough  topography 
region.  One  of  the  primary  difficulties  with  interpreting  vertical 
wave  number  spectra  of  deep  water  profiles  over  the  rough  topography 
is  that  the  low-frequency  shear  profile  cannot  be  accurately  removed 
if  only  two  profiles  have  been  obtained  at  a given  location  (at  a 
time  interval  of  one  half  of  an  inertial  period).  Any  remaining 
low-frequency  components  to  the  profiles  will  cause  errors  in  the 
calculated  magnitudes  of  the  clockwise  and  counterclockwise  spectra. 
Although  the  smoothed  hodographs  of  Figures  30  and  31  are  clearly 
rotating  counterclockwise  with  depth,  and  the  main  features  in  the 
hodographs  appear  to  reverse  direction  over  the  one-half  inertial 
period  time  interval  between  the  two  profiles,  it  is  difficult  to 
put  any  confidence  In  the  interpretation  of  spectra  calculated  from 
single  or  paired  profiles.  A time  series  over  rough  topography 
could  result  in  other  interesting  observations  as  well.  For 
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example,  one  might  speculate  that  the  observed  semidiurnal  tidal 
contribution  to  the  profiles,  which  resides  in  roughly  the  lowest 
five  stretched  vertical  wave  number  estimates  at  the  smooth  topog- 
raphy observation  side,  might  show  a change  of  character  over  the 
rough  topography.  Interaction  between  the  rough  topography  and 
"low-mode"  semidiurnal  tidal  waves  might  lead  to  the  scattering  of 
the  "low-mode"  energy  into  higher  modes.  It  would  be  interesting 
in  this  regard  to  compare  the  DLRC  calculations  (Figures  24  to  27) 
done  for  the  smooth  topography  time  series  with  the  same  calcula- 
tion done  for  a time  series  over  the  rough  topography.  If,  as 
seems  likely,  the  bottom- induced  generation  of  near-inertial  waves 
is  Intermittent  in  time  or  space,  then  the  "quick-look"  capability 
associated  with  the  EMVP  data  processing  system  would  be  an  in- 
valuable aid  in  determining  when  and  where  to  perform  the  rough 
topography  time  series. 

A second  experiment  involves  using  a pair  of  EMVP's  to  per- 
form repeated  paired  profiles,  such  as  those  shown  in  Figure  18. 
The  single  pair  in  Figure  18,  as  already  pointed  out,  has  shown 
that  there  is  an  apparent  tilt  to  energetic  features  that  seem  to 
correspond  to  each  other  in  the  two  profiles.  Repeated  paired 
profiles  would  allow  one  to  make  a much  more  certain  estimate  of 
this  tilt  and  would  also  (if  series  of  paired  profiles  were  made 
along  two  orthogonal  directions)  allow  one  to  make  an  estimate  of 
the  horizontal  wave  propagation  direction.  These  paired  profiles 


could  also  provide  information  on  whether  the  field  of  internal 
waves  is  horizontally  isotropic.  This  information  on  isotropy,  added 
to  observations  of  asymmetry  of  the  spectrum  of  internal  waves  in 
vertical  wave  number  (provided  by  the  present  work),  could  be  used 
to  test  models  that  attempt  to  calculate  momentum  fluxes  (in  the 
vertical,  for  example)  caused  by  internal  waves  (Muller,  1974). 

A final  possible  experiment  would  involve  the  combined  data 
of  an  EMVP  time  series  performed  at  some  location  (at  the  central 
moorinn,  in  smooth  topography,  for  example)  along  with  a current 
meter  mooring  at  the  same  location  which  is  densely  instrumented 
for  the  first  several  hundred  meters  off  the  bottom.  This  exper- 
iment v/ould  both  help  to  verify  the  indication  in  the  present  EMVP 
data  that  there  is  a layer  near  the  bottom  in  which  the  average, 
high-frequency,  horizontal  kinetic  eneroy  appears  to  decrease 
toward  the  bottom  and  would  also  provide  freouency  information 
which,  as  indicated  in  Chapter  IV,  is  crucial  in  determining  the 
rate  of  energy  loss  undergone  by  near-inertial  waves  as  they 
reflect  off  the  (smooth)  bottom. 

Perhaps  the  usefulness  of  a set  of  experimental  observations 
can  best  be  measured  by  whether  analysis  of  that  set  of  data  clearly 
illuminates  the  way  along  which  more  comprehensive  theoretical  and 
experimental  work  may  be  carried  out.  The  present  analysis  has 
generated  both  significant  results  and  significant  questions  which 
must  be  answered  by  further  work,  and,  therefore,  under  the  above 
criterion  this  set  of  observational  data  has  been  extremely  useful. 
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APPENDIX  A 


Profile  Locations  and 

Times 

Drop 

Date 

1973 

Drop  Time  (GMT) 

Position 

177 

April 

10 

1627 

32°21 '.ON 

64°34'.5W 

178 

II 

10 

1738 

32°21 '.ON 

64°34'.5W 

179 

II 

10 

1858 

32°21 '.ON 

64°34'.5W 

180 

li 

12 

1539 

32°21 ' .ON 

64° 34 1 .5W 

181 

II 

12 

1658 

32°21 ' .ON 

64°34'.5W 

182 

II 

12 

1833 

32°21 ' .ON 

64°34'.5W 

183 

May 

14 

1458 

29°24 ' .ON 

69°39'.8W 

184  t 

ll 

15 

1423 

28°01'.5N 

68°30'.0W 

185  t 

II 

15 

2045 

28°01'.9N 

68°32'.8W 

186  + 

ll 

15 

2238 

28°02'.7N 

68°33'.8W 

187  t 

II 

16 

0430 

28°01'.4N 

68° 30 ' . 3W 

188  t 

ll 

16 

1115 

28°03’.2N 

68° 34' .OW 

00 

VO 

ll 

16 

1750 

28°02'.9N 

68°30'.9W 

190  t 

II 

17 

0029 

27°56'.1N 

68°44 ' . 5W 

191  t 

li 

17 

0547 

27°57'.9N 

68° 39' . 8W 

192  + 

II 

17 

1305 

27°56 ' .IN 

68°44'.3W 

193  t 

II 

17 

1842 

27°58'.0N 

68° 39' .OW 

194  t 

ll 

18 

0132 

28°08'.2N 

68°28'.3W 

195  * 

II 

18 

0801 

28"08'.0N 

68°18'.0W 

196  t 

ll 

18 

1418 

28°08'.5N 

68°28'.5W 

197  t 

ll 

18 

2044 

28°08'.6N 

68°18'.3W 
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Drop 

Date 

1973 

Drop  Time  (GMT) 

Position 

198  t 

May 

19 

0231 

28°02'.5N 

68°33'.6W 

199 

II 

22 

0045 

27°06 ' ,6N 

69°40'.12W 

200 

H 

22 

1612 

28°00*.0N 

69°19'.OW 

201 

li 

23 

1512 

27°43\2N 

70°0T.8W 

202 

II 

24 

0737 

27°59'.3N 

69° 39' .OW 

203 

II 

24 

1454 

27°57 1 .ON 

69°36'.2W 

204 

II 

24 

2034 

27°59'.2N 

69°38'.0W 

205 

II 

25 

0310 

27°59'.1N 

69°38'.4W 

206 

II 

25 

0915 

27°58‘ .ON 

69° 38' .OW 

207 

II 

25 

2126 

29°21'.7N 

70°42'.2W 

208  t 

June 

1 

0204 

28°02'.0N 

68°30'.8W 

209  t 

II 

1 

0559 

28°02 1 .ON 

68°30'.9W 

210 

II 

8 

2257 

27°23'.5N 

70°01'.2W 

211 

II 

8 

2257 

27°23'.5N 

70°01'.2W 

212 

II 

9 

0530 

27°29'.8N 

70°00'.4W 

213 

li 

9 

1217 

27°23'.5N 

70°0T.2W 

214 

li 

9 

1217 

27°23'.5N 

70°01’.2W 

215 

li 

9 

1829 

27°29'.8N 

70°00’.5W 

216 

II 

9 

1829 

27°29'.8N 

70e00'.5W 

217 

II 

10 

1220 

27° 39' . 5N 

69°31'.0W 

218 

II 

10 

1220 

27°39'.5N 

69°31 1 .OW 

219  * 

II 

11 

0523 

27°59'.8N 

69°39'.1W 

220  * 

II 

11 

1152 

27°59'.7N 

69°38’.9W 

♦ 

CNJ 

CNJ 

II 

11 

1815 

27°59'.8N 

69°39'.2W 
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Drop 

Date 

1973 

Drop  Time  (GMT) 

Position 

in 

June 

11 

1815 

27°59‘.9N 

69° 39 ' .0W 

223  * 

M 

12 

0043 

27°59\7N 

69°39'.3W 

224  * 

II 

12 

0714 

27°59 1 ,8N 

69°39'.2W 

225 

II 

12 

1301 

27°52 1 .ON 

69° 39' .2W 

226  * 

II 

12 

1945 

27°59'.8N 

69°39'.0W 

227  * 

II 

13 

0219 

27°59'.8N 

69°39 1 . 1W 

228  * 

II 

13 

0759 

27°59'.8N 

69°39'.2W 

229 

H 

13 

1220 

27°51'.6N 

69°42'.9W 

230  * 

II 

13 

1800 

27°59'.7N 

69°39'.0W 

231  * 

II 

13 

2117 

27°59'.7N 

69°38'.6W 

232  * 

II 

14 

0102 

27°59'.8N 

69°38'.8W 

233  * 

II 

14 

0451 

27°59'.8N 

69°39'.1W 

234  * 

II 

14 

0828 

27°59'.8N 

69°38'.9W 

235 

It 

14 

1210 

27°57'.6N 

69°38 ' . 7W 

236  * 

II 

14 

1206 

27°59'.8N 

69°39' .1W 

237  * 

ft 

14 

1630 

27°59'.8N 

69°39 ' . 1W 

238 

II 

14 

2019 

28°03'.6N 

69°42'.0W 

239  * 

II 

14 

2016 

27°59'.9N 

69° 38' .8W 

240  * 

II 

15 

0034 

27°59'.4N 

69°39’.1W 

241 

li 

15 

0438 

28°05'.7N 

69°44 1 . 7W 

242  * 

li 

15 

0438 

27°59' .8N 

69°39'.1W 

243  * 

li 

15 

0800 

27°59'.7N 

69°39'.1W 

244 

II 

15 

1219 

28°00'.6N 

69°39'.8W 

245  * 

li 

15 

1219 

27°59'.9N 

69°39‘.0W 

146 


Drop 

Date  1973 

Drop  Time  (GMT) 

Position 

246 

June  15 

1935 

27°52'.1N  69°39'.0W 

247 

" 15 

2057 

27°52\2N  69°38'.8W 

(t)  rough  topography  drops 
(*)  time  series  drops 
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APPENDIX  B 


Use  of  Electro-magnetic  Velocity  Profilers 


to  Measure  Internal  Waves 


As  was  pointed  out  in  Chapter  I,  a conversion  equation  must 
be  used  to  obtain  horizontal  water  velocity  from  values  of  horizon- 
tal electrical  current  density  recorded  by  the  EMVP.  This  conver- 
sion equation  involves  not  only  the  above  variables,  but  also  the 
vertical  water  velocity  at  the  point  of  measurement,  as  well  as  the 
horizontal  gradient  of  any  electrical  potential  that  exists  in  the 
water.  Vertical  water  velocities  and  horizontal  potential  gra- 
dients, therefore,  introduce  errors  into  the  calculation  of  hori- 
zontal water  velocity.  We  will  now  solve  a rather  simple  problem 
that  allows  us  to  obtain  an  estimate  of  this  error. 

The  model  chosen  is  that  of  a rotating  ocean  with  constant 
Brunt-VSisSlS  frequency  and  constant  depth,  H,  and  with  a rigid 
surface  and  bottom.  We  assume  that  internal  wave  vertical  modes 
exist  in  the  water  column.  This  assumption  is  not  crucial  in  esti- 
mating the  error,  as  will  be  seen  later. 

The  equation  for  the  electrical  potential,  4>,  and  the  elec- 

A ^ 

trical  current  density,  J = (Jx,  Jy,  Jz),  generated  by  an  ocean 
moving  with  velocity  $(x,y,z,t)  = (u,v,w)  is 

A 

* 2 + 

V<J>  = v x F - J/o  , 


B-l 


r 
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where  o Is  the  conductivity  of  the  sea  water,  which  is  assumed  to 
be  constant,  and  t = ( 0 , Fy , Fz ) Is  a vector  representing  the  terres- 
trial magnetic  field.  The  coordinates  x,  y,  and  z are,  respectively, 
magnetic  east,  magnetic  north,  and  vertical.  The  coordinate  origin 
is  located  at  the  sea  surface,  and  z is  positive  upward. 

Equation  B-l  will  be  called  the  full  conversion  equation, 

/V  A 

„ ■+  ■* 

since  it  contains  all  terms  that  relate  #,  v and  J.  The  horizontal 
components  of  equation  B-l  are: 


J 

-X 

a 


-uF  - M 
z 3y 


B-2 


The  underlined  terms  make  up  the  "applied"  conversion  equation,  or 
the  equation  that  is  actually  used  to  convert  EMVP  electrical  cur- 
rent readings  to  horizontal  water  velocity.  The  remaining  terms 
cause  the  error.  (The  reader  should  refer  to  Sanford,  et.  al. 
(1974)  for  a more  comprehensive  discussion  of  EMVP  measurements.) 

A 

Since  V • 3 = 0 and  V x F = 0,  the  equation  to  be  solved 

becomes 

V2<j>  = V • (v  x F)  = F • (v  x v)  B-3 

with  the  boundary  conditions: 
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A 

w = Jz  = 0 at  z = 0; 

w = Jz  = 0 at  z = -H  . B-4 

A 4 

Evaluating  7 x F,  these  boundary  conditions  may  be  rewritten: 
w(z  = 0)  = w (z  = -H)  = 0 ; 


(z  = 0)  = F • u(z  - 0)  ; 

dZ  J 

A 

& (z  = -H)  = F • G(z  = -H)  . B-5 

9z  y 

A A 

Given  $(x,y,z,t)  we  can,  at  least  in  principle,  solve  for  4>  and 

A 

then  3. 

The  internal  wave  equations  are: 


A 

3u 

_ 
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II 

< > 
4- 

at 

ax 

A 

a£ 
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fG . - ^ 

at 

»y 

aw 

- i 
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32 

A 

ab 

+ 
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at 

/A 

au 

+ 

av  + aw  _ 0 

ax 

ay  32  " 

where  u,  v and  w are,  respectively,  the  east  magnetic,  north  mag- 
netic and  vertical  velocity  components,  f is  the  Coriolis  parameter. 
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and  p is  pressure  divided  by  a reference  density,  pq.  We  have  used 
the  Cousslnesq  approximation,  in  which  the  density,  p , is  assumed 
to  be  given  by  three  terms: 

p = pQ  + p(z)  + p(x,y,z,t) , B-7 

and  the  water  is  assumed  to  be  incompressible.  The  water  pressure 

A 

is  made  up  of  two  parts,  p0p  and  P0p.  The  hydrostatic  part  of  the 
pressure  is  given  by 

p <Mzl=  _g(-  +£(z)),  B-8 

0 dz  0 

where  g is  the  acceleration  due  to  gravity.  Then  p0p  is  the  devia- 
tion of  the  pressure  from  hydrostatic,  the  buoyancy,  b,  is  -gp/pQ 
and  the  Brunt-Vaisala  frequency,  N,  is 

N=(--^)1/2  B-9 

p dz 

By  Fourier  transforming  equations  B-6,  where 

~i  i .l  ■>  \ ^kx  + ■ wt) 

u(x,y,z,t)  ■ u(z;k,l,u))  e , B-10 


and  similarly  for  the  other  variables,  we  obtain  the  following 
equation  for  w: 

d2w  + (kJ  + l2)  .N.2._-.  .oi.2,,  w = o.  B-l  1 

dz2  w2  - f2 

The  form  of  w that  satisfies  the  boundary  conditions  on  vertical 
velocity  at  the  surface  and  bottom  is 
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nirz 

~ _ A sin  H 
w = 

M 

H 

with  n an  integer,  and 


2-2 


n“TT 


(k2  + l2) 


2 _ , .2 


2 x 1 2 \ N - 0) 


U>2  - f2 


Equations  for  u and  v are: 


*c: 


mrz 


u = - 


v = - 


A 

N2  - 

w2 

COS  H 

nir 

iu)(ui2 

- f2) 
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H 
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mrz 

A 

N2 

- a)2 

cos  H 

(ilf  + k«*>); 


nir  iu)(w2  - f2)  n^ 
H H 


(lw  - ikf) 


Then 


- (k2  + l2)  $ = F • (V  x v)  = 


dz 


= F 


nrcz 
sin  H 

y n^ 

H 


fA 

\ iu>(u)2  - fJ 


(ilf  + ku>)  - ikA 


where 


+ F,  -C0S 


mrz 

H 


2 mr 
H 


Bf (k2  + l2)  L 


B = - ^--N2.  - u)2 


Djl  ia)(u)2  - f2) 
H 


B - 1 2 


B-13 


B-l  4 


B-l  5 


B-16 


B-l  7 


V 


The  boundary  conditions  on  $ are  then: 


and 


■B  _ * k"> 

dz  mr 

H 


at  z * 0 


...  F B(-l)n(11f  + ku) 
a<t  _ y 

dz  “ n^ 

H 


at  z 


-H 


B-18 


B-19 


The  system  of  equations  B-16,  B-18  and  B-19  can  now  be 
solved  to  give: 


$= 


UL 


N2  - u>2 


^rTT2(2^-+  k2  + l2)  \ 

H 


- 7717  - ik“>  • ik 


FyA(N2  -a)2)  (If  - ikoi) 
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mrz 
sin  H 


01 
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FzAf(Nz  - a,2)  (k2  ♦ l2)  cos' 


3(a)2  - f2) 


mr._ 

H 


E-l 


The  full  conversion  equations  have  been  given  above  (equa- 
tions B-2).  The  applied  conversion  equations  (those  actually  used 

* /v 

to  convert  from  3 to  $)  are 


J * 

= (v  - v ) Fz  . 
o 

JJL  ’ - <Q  - i*)  Fz  , 
a £ 


where  ¥*  and  v*  are  depth-independent  terms  (Sanford,  1971). 

A 

Since  we  are  only  interested  in  the  part  of  1 that  varies  with 


depth,  we  set 


— * 

¥ = v = 0. 


If  we  now  let  u and  v be  the  velocity  components  as 
true  true 

given  by  equations  B-14  and  B-15  and  uQbs  and  v ^ be  the  compo- 
nents calculated  from  the  applied  conversion  equations  B-21  and 


B-22 , then: 


Z - U =1M 

true  obs  fz  3y  } 


’true  ‘ * obs  ' V ( If  * "V  • B-24 

We  see  that  the  expressions  for  $ and  w can  be  used  to  estimate  the 
error  in  the  conversion  equations.  At  a given  depth,  2,  the  energy 
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In  the  u error  signal  is  given  by 


Eu  (z)  = 1/2  < (G  - u . ) (u  - u ) > 
error  true  obs'  'true  obs* 


2F 


]T<  l^|2> 


3y 


B-25 


where  the  brackets  < > indicate  an  average  over  time.  Similarly, 
the  v error  signal  is  given  by 


= ^ < ) (V*  - V ) > 

error  true  obs  true  obs 


3$ 

3x 


+ v?  F 

y 


> 


We  note  that 


and 


Eu^  (z)  = 1/2  < u • u*  > 
true  true  true 


Ev  (z)  * 1/2  < v 
true  true 


* 

v > 
true 


B-26 


B-27 


B-28 


Returning  momentarily  to  the  expression  for  $ (equation 
B-20)  we  see  that  the  first  two  terms  are  proportional  to 
cosh  *'K+1*  z and  (sinh  {A2  + l2  z)/  A2  + l2).  We  now  make  the 
assumption  that  the  horizontal  component  of  the  wave  number  vector, 
A2  + l7,  is  small  enough  so  that  the  cosh  term  is  independent  of  z 
and  the  sinh  term  becomes  a linear  function  of  z.  Since  we  are  only 
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Interested  In  that  part  of  <p  that  is  depth-dependent,  we  ignore 
the  hyperbolic  cosine  term.  The  term  (sinh  A2  + l2  z)/(A2  + 1*)  ■* 
z as  A2  + l2  0.  If  we  examine  the  coefficient  multiplying 
(sinh  A2  + l2  z)/(^2  + l2)  we  see  that  for  k2  + l2  « (n2n2)/(H2) 
the  last  term  In  the  coefficient  cancels  the  first  term,  and  we  are 
left  with  (-i  kFyA)/(^Mj— ) as  the  coefficient  of  the  term  propor- 
tional to  z.  For  k = 0 this  term  disappears.  Even  if  k f 0 this 
term  decreases  rapidly  as  n increases;  that  is,  for  higher  modes  or 
for  frequencies  close  to  f.  Since  we  are  only  interested  in  those 
terms  that  are  sinusoidal  in  z,  we  will  ignore  the  terms  in  $ that 
are  either  independent  of  z or  are  only  weakly  dependent  on  z. 


I I IN*  I I 

* = ’ k*  ♦ v - f!> 0 


f - ikw)  - ik 


mu 

'l  sin  H 
j mr 


F Af (N2  - a)2 ) (k2  + l2) 


ttbU*  - f‘> 

H 


cos  H 


(The  reason  for  this  assumption  has  been  explained  in  Chapter  II. 
Any  depth-independent  contributions  to  the  profiles  examined  in 
Chapter  II  have  been  removed  prior  to  analysis.  Thus,  the  (approx- 

A 

imately)  depth-independent  contributions  of  <j>  should  not  appear  in 
the  profiles.) 

We  are  particularly  interested  in  the  relative  error  in 
u_h(.  and  v „Ke,  given  by  Eu /Eu.  ^ and  Ev /Ev  , 


error  true 


error  true 


respectively.  Tor  simplicity  , we  have  averaged  the  quantities 

[u  , Ev  , Eu  , and  Ev  over  depth,  z,  from  C to  -H. 
error  error  true  true 


Then 


Eu  = — ( < u..  • u.  > dz, 

true  H J true  true 


B-30 


and  similarly  for  the  other  three  quantities.  The  final  result  is 

that  Eu  /Eu  and  Ev  /Ev  are  given  by: 
error  true  error  true 
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B-31 


B-32 


Considering  first  the  u-component  error,  it  can  be  shown 
that  Tu  is  a maximum,  for  a given  <o  and  n,  if  the  wave  is  propa- 
gatinq  in  the  north-south  direction  (k  = 0).  Further,  for  waves 
traveling  in  the  east-west  direction  1 = 0,  and  Tu  = 0.  For 

l2  «(°^)2,  we  find  that  max  {Tu}  — - — , if  the  frequency,  to, 

(“)2 
1 H' 

is  close  enough  to  f so  that  (u>2  - f 2 ) / (N2  - w2)  is  small.  Since 


l2  _ to2  - f2 
(™)2  N2  - u)2 


from  equation  B-13, 
max  (Tu) 

N - w 


B-33 


B-34 


Equation  B-33  or  B-34  shows  that  the  relative  error  in  the  u-compo- 
nent energy  depends  only  on  frequency,  u>,  or  (what  amounts  to  the 
same  thing)  on  the  north-south  "aspect  ratio"  of  the  wave,  l/(mr/M). 

Similar  statements  may  be  made  about  relative  error  in  the 
v-component  energy.  In  this  case,  Tv  is  a maximum,  for  a given  to 
and  n,  if  the  wave  is  propagating  in  the  east-west  direction  (1  = 0). 
However,  Tv  does  not  go  to  zero  as  k goes  to  zero,  because  of  the 
presence  of  the  term  involving  vertical  velocity  in  the  full  conver- 
sion equation  for  v.  If  (to2  - f 2 )/ ( N2  - to2)  is  small  we  find  that 

k2 


max  {Tv}o( 


(nir/H)2  . 


B-35 


As  was  the  case  for  Tu,  max  {Tv}  for  small  (id2  - f 2 )/ (N2  - id2) 
depends  only  on  the  frequency,  or  al ternati vely,  on  the  east-west 
aspect  ratio  k/(mr/H). 
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Using  values  of  the  Brunt-VSisSla  and  inertial  frequencies 
appropriate  to  the  MODE  center  at  a depth  of  4500  dbar,  we  find 
that  the  relative  error  terms  max  {Tu}  and  max  {Tv}  are  both  less 
than  5%  for  periods  down  to  about  one-half  day.  This  range  of 
periods  from  inertial  down  to  one-half  day  contains  roughly  three 
fourths  of  the  energy  observed  in  internal  wave  spectra  from  the 
MODE  area. 


We  have  ignored  vertical  variations  of  N and  the  fact  that 
the  observed  waves  are  propagating  in  the  vertical.  Since  the 
error  terms  are  proportional  to  (id2  - f 2 )/ (N2  - <d2),  and  this 
would  not  change  if  waves  were  propagating  vertically,  the  presence 
of  vertically  propagating  waves  would  not  change  the  error  estimate 
significantly. 

The  Influence  of  varying  Brunt-Vaisala  frequency  is  somewhat 
more  difficult  to  determine.  If  N varies  weakly  with  z,  this  means 
the  forcing  terms  on  the  right  hand  side  of  equation  B-17  have 
coefficients  that  are  weak  functions  of  z.  If  we  assume  that  N 
varies  an  a scale  that  is  appreciably  longer  than  the  scale  of  the 
waves  and  try  to  solve  the  forced  problem  using  two  depth  scales, 
we  would  find  that  the  lowest  order  estimate  of  the  relative  error 
terms  would  again  be  proportional  to  (<d2  - f 2 )/ ( N2  - id2)  but  that 
now  N is  a weak  function  of  z.  The  value  of  N used  above  to  get 
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relative  error  estimates  was  taken  to  be  that  at  4500  dbar  because 
this  value  is  close  to  the  minimum  average  value  in  the  N profile 
and  thus  presumably  maximizes  the  relative  error. 
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APPENDIX  r 


Accuracy  and  Precision  of  Observations 


It  is  important  to  distinguish  between  accuracy  and  preci- 
sion when  interpreting  velocity  mearurements  obtained  by  the  EMVP. 
This  is  true  because  of  certain  characteristics  of  the  EMVP  it- 
self. A detailed  description  of  the  arrangement  of  sensors  on  the 
EMVP  tan  be  found  elsewhere  (Sanford,  et.  al.,  1974).  It  has  been 
pointed  out  (Sanford,  10/1)  that  the  electrodes  that  sense  hori- 
zontal electric  current  density  only  measure  the  depth-varying  part 
of  the  current.  That  is,  there  is  a depth- independent  part  of  the 
current  that  is  not  sen'ed  by  the  EMVP.  Therefore,  the  instrument 
measures  horizontal  current  in  much  the  same  way  that  geostrophic 
shear  calculations  do,  when  a "level  of  no  motion"  must  be  assumed. 
This  is  clearly  a loss  of  accuracy  in  the  instrument. 

A second  cause  of  loss  of  accuracy  in  the  EMVP  involves  the 
actual  location  of  the  electrodes  on  the  instrument.  As  the  pro- 
filer falls.  It  rotates.  This  rotation  causes  an  oscillatory  cur- 
rent to  flow  in  a coil  within  the  profiler,  since  the  coil  is 
rotating  in  the  terrestrial  magnetic  field.  This  oscillating  cur- 
rent is  then  used  to  determine  the  point  in  a cycle  of  rotation 
when  the  coil  is  pointing  toward  magnetic  north.  The  orientations 
of  the  two  independent  sets  of  electrodes  on  the  profiler  are  de- 
fined relative  to  the  direction  of  this  "compass  coil."  In  the 
present  configuration  of  the  instrument,  one  set  of  electrodes  is 
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oriented  parallel  to  the  compass  coil,  and  one  set  is  oriented 
perpendicular  to  It.  However,  these  "perpendicular"  and  "parallel" 
orientations  relative  to  the  compass  coil  can  only  be  measured  to 
about  a degree.  The  result  is  that  errors  in  defining  the  orien- 
tation of  the  electrodes  on  the  instrument  cause  errors  in  the 
determination  of  the  magnetic  bearing  of  the  horizontal  electric 
current  density  vector,  and  thus  errors  in  the  measured  water  ve- 
locity. This  is  a systematic  error,  which  is  another  cause  of  loss 
of  accuracy  in  the  measurements.  If,  in  a small  section  of  a pro- 
file, velocity  readings  from  the  two  independent  electrodes  are 
compared,  it  is  usually  found  that  although  the  two  velocity  read- 
ings track  each  other  closely  with  changing  depth,  there  can  be  a 
constant  offset  between  the  two  of  up  to  1 cm/sec.  If  this  offset 
is  removed,  the  actual  RMS  difference  between  the  two  readings  is 
usually  less  than  0.5  cm/sec.  This  RMS  error  can  be  taken  to  be 
an  estimate  of  the  precision  of  tne  velocity  measurements. 

A more  detailed  discussion  of  EMVP  velocity  measurement 
errors  can  be  found  in  Sanford  et.  al.  (1974).  The  above  sources 
of  error  should  be  kept  in  mind  when  interpreting  velocity  profiles. 

As  mentioned  in  the  introduction,  the  EMVP  also  carries 
pressure,  temperature  and  conductivity  sensors.  We  will  consider 
the  pressure  sensor  error  first.  Because  data  in  the  EMVP  is  re- 
corded in  digital  format,  there  is  a quantizing  error  in  any  re- 
corded variable.  In  the  case  of  pressure,  this  quantizing  error 
amounts  to  about  1.5  dbar.  Bottom  pressure  values  as  obtained  by 
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the  EMVP  for  each  drop  in  the  time  series  have  been  compared  in 
order  to  estimate  the  precision  of  pressure  measurements.  If  we 
assume  that  the  bottom  at  the  time  series  site  is  flat,  so  that  the 
EMVP  hit  bottom  at  the  same  depth  in  each  drop,  then  a comparison 
of  bottom  pressures  as  recorded  by  the  EMVP  indicates  a root-mean- 
square  error  of  about  + 1 dbar  at  about  5500  dbar.  This  is  a meas- 
ure of  the  precision,  or  repeatability,  of  the  pressure  measurements. 
On  the  other  hand,  accuracy  is  limited  by  the  fact  that  values  of 
pressure  computed  from  PGR  (Precision  Graphic  Recorder)  readings 
could  be  determined  to  only  about  5 dbar. 

Temperature  and  conductivity  measurements  were  notably  less 
successful  than  pressure  measurements.  The  quantizing  errors  for 
temperature  and  conductivity  are  0.017°  C and  0.02  mmho/cm.  The 
quantizing  error  in  temperature  is  twice  what  it  should  have  been, 
since  the  least  significant  bit  in  the  digitizer  was  usually  stuck 
in  the  "1"  mode.  In  addition,  due  to  a failure  of  the  conductivity 
head  on  the  profiler  (the  quartz  liner  in  the  head  cracked  under 
high  pressure  and  fell  out),  the  conductivity  measurements  of  the 
head  were  seriously  influenced  by  temperature  and  pressure.  This 
made  the  computation  of  salinity  difficult.  Although  a curve  was 
derived  which  corrects  conductivity  for  the  above  influences  of 
temperature  and  pressure  on  the  head,  the  resulting  computed  salin- 
ity was  only  accurate  to  about  0.035  °/oo.  This  was  not  judged  to 
be  sufficient  for  density  computations,  and  so  these  measurements 
have  not  been  used  in  the  present  work. 


Finally,  a comment  should  be  made  on  navigational  accuracy 
during  the  time  series.  Ship  positions  were  determined  by  Loran  C. 
Although  nominally  performed  at  one  location,  a comparison  of  re- 
corded ship  positions  for  all  drops  during  the  time  series  indicates 
that  the  ship  positions  have  a standard  deviation  of  about  + 230 
meters  (east-west)  and  + 270  meters  (north-south)  around  a mean 
position  of  27°  59.75'  N,  69°  39.04'  W.  It  should  be  noted  that 
the  above  numbers  are  based  on  ship  position  at  the  beginning  of 
each  drop.  The  EMVP  usually  spent  5 to  10  minutes  on  the  surface 
before  it  began  its  descent.  During  this  time  the  instrument  was 
drifting  free  of  the  ship  and  was  presumably  carried  to  the  south 
by  the  mean  flow  (see  Figure  5).  This  additional  drift  of  the 
instrument  to  the  south  would  result  in  a shift  of  the  mean  posi- 
tion of  the  time  series  to  the  south  (probably  by  about  100  to  200 
meters)  but  presumably  would  not  increase  the  standard  deviation  of 
the  drop  positions  appreciably. 
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APPENDIX  D 

A More  General  Treatment  of  the  Partition  of  Energy 
Between  Clockwise  and  Counterclockwise  Spectra 


We  let  m be  the  vertical  wave  number  in  the  stretched  verti- 
cal coordinate,  f be  the  inertial  frequency,  and  w be  the  wave  fre- 
quency. In  what  follows,  we  will  use  the  vertical  coordinate  s, 
where  s = -z  . Thus  s is  a stretched  vertical  coordinate  increas- 
ing upward.  The  horizontal  velocity  components  of  a single  internal 
wave  may  be  represented  by: 


LI  (s,t;uj,m) 

L 


e 


2Tti(ms  - wt) 


f 2iri(ms  - wt) 
U (s,t,w,m)  = - e 

T W 


where  UL  and  Ux  are  the  horizontal  velocity  components  parallel 
and  perpendicular  to  the  direction  of  horizontal  phase  propagation, 
respectively.  Then  (Garrett  and  Munk,  1972),  the  east  and  north 
(u  and  v)  velocity  components  are  given  by 


and 


u = i cos0  - Ux  s i n<p 
v = i Ul  sin<))  + UT  cos0  , 


D-2 


where  0 is  the  horizontal  phase  propagation  direction.  We  adopt 
the  convention  that  m is  always  positive.  A wave  with  upward  phase 
propagation  then  has  a positive  frequency,  w,  while  a wave  with 
downward  phase  propagation  has  a negative  frequency.  All 
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contributions  to  the  profiles  from  upward  propaqatino  waves  may 

then  be  represented  by 
N 00 

o . ( f ./  2Tri(ms  - ait) 

L ( s , t ; <j> ) = J ) A(w,m,<j>)  e 


U 

and 


dm  dw 


f 0 
N °° 


lT  (s,t;$)  • j f l A(oj,m;<j))  e 


2iri(ms  - wt) 


dm  dm 


D-3 


f 0 


The  amplitude  function,  A,  is  considered  to  be  non-zero  only  for 
f < uj  < N and  m > 0.  Therefore,  the  double  integrals  can  be  con- 
sidered to  be  Integrals  from  -«  to  *».  We  similarly  represent  the 
contributions  to  the  profiles  from  downward  propagating  waves  as 


Ul  (s,t l<p)  = 


2iri  (ms  - wt) 

B(w,m;4>)  e dm  dio 


and 

00  OO 

o ( ( f , , 2iri(ms  - ut)  . . 

U°  (s,t; ♦)  - J ) ^B(u>,m;*)e  dm  do)  . D_4 

-00  ”00 

In  this  case,  B is  non-zero  only  for  -f  < ui  < -N  and  m > 0. 

The  total  horizontal  velocity  components  are  then  given  by 
U = U°  + U°  and  II  * U + U . D"^ 

i_  l v.  T T T 

The  cross  correlation  between  u and  v and  the  autocorrela- 
tion for  the  u and  v components  are 


j 

] 
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p-uv  (s»t;4>)6(T)  = « |u(s,t)v*(s  + s,t  + T) 

+ u*(s,t)v(s  + s.t  + T)^  6 ( t ) 6 ( t - T)  >t>s 
and 


D-6 


[p  ( s , T ; d>)  + p (s,T;<|>)]  6 ( T) 

>w  V v 

= «^u(s,t)u*(s  + s\t  + T)  + u*(s,t)u(s  + F,t  + T) 

+ v(s,t)v  (s  + S\t  + T)  + v (s,t)v(s  + ?,t  + T)]  6(t)6(t  - T)  > > 

J t s 

D-7 

where  the  angle  brackets  denote  averages  over  s and  t.  The  6-func- 
tions arise  because  we  are  making  a linear  section  of  the  random 
wave  field  along  the  s axis  (Konyaev,  1973).  (The  contributions 
to  u and  v at  a given  frequency  and  vertical  wave  number  are 
assumed  to  be  a sum  of  waves  with  random  phases.) 

We  assume  that  the  above  integral  limits  may  be  taken  as  - °° 
to  + °°.  This  implies  an  infinitely  deep  water  column.  Strictly 
speaking,  the  limitation  to  a finite  depth  will  cause  the  vertical 
wave  number  spectra  calculated  below  to  be  smoothed  by  a window 
which  is  a function  of  (stretched)  vertical  wave  number.  The  finite 
depth  limitation  has  been  ignored  in  what  follows.  Using  the  above 
expressions  for  the  correlation  functions,  we  can  calculate  the 
cospectrum  (P),  quadrature  spectrum  (Q),  and  the  total  energy 
spectrum  (TE)  for  ze  ) time  lag.  These  quantities  are  given  by 
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P(m;4>)  + 1 0 (m ;<t>)  = j j v (s  ,T;4>)S(T)e 


-2ni  (ms  - utT) 

ds  dT 


( ( _ -2iTi(ms  - wT)  _ 

TE (m ;4> ) = \ (p  + p ) (s ,T;$)<5(T)e  ds  dT 


Carrying  out  the  required  calculations,  we  find 

N 


P(m;<J>)  + iQ(m;4>)  = 1/4  jcos4>sin<J>  [ (1  - | A(u),m;<|>)  1 2dw 

^ f 

-N 

+ j (1  - ) I B(io,m;4>)  1 2 dw 


n -n 

-ill  - | A(a),m;4>)  1 2 <fc>  + I - | B(to,m;d>)  | 2 du) 

L.  ' u (i) 


TE(m;4>)  = 1/i  j ) |A(m,ur,<|>)  | 2 (1  + ) 


I'l 

j"  |B(w,m;<t>)  I2  + ^ ^ du>  ^ 


We  can  simplify  these  expressions  if  we  assume: 


J 
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|A(m,u>;<j>)  | 2 = (1  - a(m,w' ) )E(u>’ )A(m)R(4>)  ; 

| B(m,oo;<j>)  | 2 - afm.io1  )E(u>'  )A(m)R(<>) , D-12 

where  u>’  > 0.  The  value  of  a(m,u)')  is  between  zero  and  one.  Use 
of  this  coefficient  allows  us  to  indicate  how  much  energy  is  con- 
tributed to  a one-sided  frequency  spectrum  by  waves  with  upward 
and  downward  phase  propagation.  R(<t>)  is  assumed  to  be  normalized 

A 

so  that  its  integral  around  2ir  is  one.  E(gj')  and  A(m)  could  be 
identified  with  the  frequency  and  vertical  wave  number  dependence, 
respectively,  of  the  Garrett-Munk  spectrum  given  in  Chapter  V. 
Integrating  around  <p,  we  get 
N 

TE(m)  = 1/4  f E(u')A(m)  (— ' d«'  D-13 

/ 0) 

f 

and 

N 

Q(m)  * 1/4  I E(<*)'  )A(m)  (^)  (1  - 2a)  do'  . D-14 

' oo 

f 

The  clockwise  and  counterclockwise  spectra  are  given  (to 
within  a constant)  by: 

C(m)  = TE(m)  + 2Q(m)  ; 

A(m)  =»  TE(m)  - 2Q(m)  (Gondlla,  1972).  D-15 

We  point  out  that  neither  TE  nor  Q contains  multiplicative  factors 
involving  cos<f>  or  s 1 nd>  in  equations  D-10  and  D-ll.  These  quantities 
(and  therefore  C(m)  and  A(m)  do  not  depend  on  the  distribution  of 


R( 4> ) around  a circle.  This  Is  another  way  of  saying  that  C(m) 
and  A(m)  are  Invariant  with  respect  to  a rotation  of  the  hori- 
zontal coordinates.  If  we  now  assume  that  E(w’)  has  a sharp  peak 
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at  a frequency  co ' near  the  inertial  frequency,  so  that  w'/f  ~ 1, 
o o 

then 

A(m)  _ 1 - g(m,  u>p) 

C(m)  a(m,  u/)  ’ D_16 

which  is  the  ratio  of  downward  to  upward  energy  propagation.  For 
a small  a (which  means  that  most  of  the  energy  is  propagating  down- 
ward), A(m)  is  much  larger  than  C(m).  It  is  important  to  remember 
that  s is  positive  up.  If  we  were  to  use!*,  the  stretched  verti- 
cal pressure  coordinate  which  increases  downward  (as  for  the  spec- 
trum in  Figure  13,  for  example),  then  C(m),  the  clockwise  spectrum, 
would  dominate.  (Changing  the  vertical  coordinate  from  s to  7*  is 
equivalent  to  inter-chanqinq  A and  C.)  Several  other  points  can 
be  made.  First,  if  a = 1/2  everywhere,  the  quadrature  spectrum  will 
be  zero,  and  clockwise  and  counterclockwise  energy  will  be  equal. 

The  choice  of  a = 1/2  corresponds  to  the  case  when,  on  the  average, 
equal  energy  is  contained  by  waves  with  upward  and  downward  energy 
propagation.  Second,  if  an  energetic  wave  is  present  with  a fre- 
quency much  greater  than  inertial  (such  as  the  semidiurnal  tide), 
we  might  expect  that  the  addition  of  this  wave  will  cause  the  ratio 
A(m)/C(m)  to  approach  one.  It  is  interesting  to  note  in  this  con- 
text  that  the  estimates  of  A(m)/C(m)  (using  z as  the  vertical 


I 


coordinate)  at  the  smallest  vertical  wave  numbers  in  Figure  22 
(where  we  know  that  the  semidiurnal  tide  is  important)  do  seem  to 
increase,  relative  to  estimates  at  neighboring,  higher  stretched 


vertical  wave  numbers. 
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